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Preface 

The  road  to  the  end  was  a  bit  tortuous.  I  set  out  to 
compare  the  rate  of  convergence  to  solution  of  five  or  six 
iterative  methods  of  solving  the  Implicit  finite  difference 
representation  of  the  equation  of  transient  heat  transfer, 
to  do  an  analysis  on  the  effect  on  the  rate  of  convergence 
of  a  spiral  grid  soan,  and  a  detailed  error  analysis.  To¬ 
ward  the  last  third  of  the  allotted  time,  the  relationship 
between  the  adapted  Wegsteln  technique  and  successive  over¬ 
relaxation  appeared  and  the  deoision  was  made  to  follow  this 
line  of  investigation.  To  this  end,  this  relationship  be¬ 
tween  the  two  methods  is  now  rather  obvious,  though,  in  my 
earlier  work  on  this  thesis  it  was  not  at  all  obvious  to  me. 
Time  did  not  permit  a  thorough  investigation  designed  to  use 
this  relationship  for  possible  Improvement  of  the  adapted 
Wegstein  technique. 

I  wish  to  express  my  appreciation  to  Dr.  Bernard  Kaplan 
of  the  Physics  Department  of  the  Air  Force  Institute  of 
Technology  for  his  assistance  and  guidance.  Host  of  all,  my 
gratitude  (and  more)  goes  to  my  wife  who  still  insists  that 
having  me  constantly  off  in  my  litter-strewn  corner  of  the 
house  is  better  than  having  me  away  on  te^orary  duty  some¬ 
where. 


Robert  Poppe 
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Absti^ot 

Tills  study  investigu tes  tae  practical  application  of 
two  convergence  techniques  designed  tc  increase  tae  rate 
of  convergence  of  the  method  of  successive  displacements 
( Gbuss-wledel )  for  tae  implicit  numerical  solution  of  the 
diffusion  equation  of  transient  neat  transfer,  a  sample 
problem  of  determining  tae  temperature  distribution  In  a 
cube  with  a  constant  internal  neat  source  and  fixed  bounda¬ 
ry  temperatures  Is  solved  to  provide  the  necessary  data. 

The  results  provide  a  theoreticel  basis  ror  tne  adapted 
•ve ostein  tecnnioue  that  was  not  previously  availacle.  ihis 
tneoretical  basi3  brings  to  light  the  feet  that  successive 
overrelaxation  and  the  adapted  .vegstein  technique  are  cased 
on  the  same  theoretical  o&cjcgrcund. 

A  procedure  based  on  estimating  the  maximum  eigenvalue 
of  tne  metnoa  of  successive  displacements  is  used  to  make 
an  approximation  of  the  relaxation  factor  for  successive 
overrelaxation.  This  procedure  is  shown  to  be  a  practical 
metnou  of  finding  the  proper  relaxation  factor  to  estimate 
the  diffi cult- to-determine  optimum  factor.  The  savings 
using  tnis  procedure  was  about  5C>  of  the  iterations  re¬ 
quired  to  obtain  tne  same  solution  by  successive  displace¬ 
ments. 

A  comparison  of  the  two  accelerating  techniques  is 
made.  Items  of  comparison  are:  the  number  of  iterations 


viii 


GA/rhys/63-8 


Abstract 

required  to  obtain  a  solution  of  the  sample  problem,  the 
progress  of  the  solution  with  increasing  iterations,  the 
characteristics  of  the  accelerating  factor  determined  while 
computing  in  the  normal  successive  displacements  mode,  and 
the  error  associated  with  the  solution  of  the  sample  problem. 

The  principal  results  of  this  comparison  are: 

1.  A  65*  reduction  in  iterations  over  successive 
displacements  when  the  optimum  relaxation  factor  is  used 
and  about  50*  and  30*  reductions  for  successiver  overrelax¬ 
ation  with  an  estimated  relaxation  factor  and  the  adapted 
rfegstein  technique,  respectively. 

2.  Successive  overrelaxation  produces  a  smooth 
convergence  to  solution,  whereas,  the  adapted  ftegsteln  tech¬ 
nique  is  ragged. 

3.  Both  the  single  maximum  eigenvalue  used  to 
compute  the  relaxation  factor  of  the  method  of  successive 
overrelaxation  and  the  #egsteln  slopes  used  to  oompute  an 
accelerating  factor  for  each  individual  node  oonverge  to  the 
same  value,  but  at  different  rates. 

4.  The  error  associated  wltn  the  accelerated 
solutions  is  less  than  that  encountered  in  the  method  of 
successive  displacements. 

An  additional  investigation  is  made  of  the  effects  on 
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Abstract 

the  rate  of  convergence  of  successive  displacements  of 
scanning  the  finite  difference  grid  in  a  spiral  mode  from 
tne  fixed  boundary  conditions,  ±ne  rate  of  convergence  is 
compared  to  that  obtained  using  a  conventional  serial  scan¬ 
ning  procedure.  The  results  obtained  snow  a  small  decrease 
in  the  number  of  iterations  required  for  the  spiral  mode 
over  the  serial  scan,  but  tne  savings  are  not  significant. 

A  tentative  conclusion  is  reached  that  no  change  ii.  scan¬ 
ning  procedure  from  serial  scan  will  produce  a  significant 
decrease  in  the  iterations  required  to  solve  the  transient 
heat  transfer  problem  by  successive  displacements. 
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I.  Introduction 


Background 

The  analytical  solution  to  partial  differential  equa¬ 
tions  of  engineering  and  applied  physios  can  not  in  general 
be  obtained  exoept  in  very  apeoial  cases.  Thus,  approximate 
methods  of  solution  have  been  developed,  and  the  most  popu¬ 
lar  of  these  procedures  is  the  method  of  finite  differences. 
Prior  to  the  development  of  large  scale  digital  computers, 
practical  numerical  solution  of  the  finite  difference  equa¬ 
tions  was  almost  entirely  performed  by  pencil  and  paper  re¬ 
laxation  methods.  When  computers  beoame  available,  it  was 
found  that  the  original  relaxation  methods  were  not  as  suc¬ 
cessful  as  some  systematic  computation  of  points  in  a  con- 
vienent  eyello  order  (fief  4:848).  Is  a  result,  the  develop¬ 
ment  of  new  and  the  adaptation  of  old  solution  procedures 
for  use  with  digital  computers  has  been  the  subjeot  of  in¬ 
tensive  study  in  the  past  ten  to  fifteen  years. 

The  state-of-the-art  at  the  present  time  provides  a 
number  of  digital  oomputer  routines  for  solution  of  the 
finite  difference  representation  of  the  oommon  types  of 
partial  differential  equations  encountered  in  engineering. 
For  instance,  the  finite  differenoe  representation  of 
Poisson's  equation,  V*u=$(ul)  for-  U.=  $(*>%)  in  the  rect¬ 
angular  region  fi  with  U  described  on  the  boundary  a—  the 
equation  encountered  in  suoh  diverse  engineering  routines 
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as  elasticity,  heat  transfer,  electromagnetism  and  fluid 
dynamics— has  at  least  ten  possible  representations  for  use 
in  the  solution,  (rtef  4:282).  The  theoretical  bests  for  these 
diverse  methods  arc  **11  developed  for  the  simple,  two  dimen¬ 
sional,  rectangular  geometry  involved.  The  theoretical  basis 
for  extending  these  routines  into  cases  of  irregular  geometry 
and  three  dimensions  is  still  rather  fragmentary.  The  state 
of  development  does,  however,  permit  an  experimental  approach 
to  the  investigation  of  possible  extensions  to  the  existing 
theory. 

The  specific  case  of  the  solution  of  the  equation  of  non- 
steady  heat  transfer  represents  an  important  computer  routine 
in  a  given  engineering  facility,  in  three  dimensions,  this 
equation  is  of  the  form 

-£■  V-  KVU(x^,i,jL)  *  U(x,yftt)  +  ftx,<0,2,J 1) 

The  presence  of  tne  time  dependence  permits  formation  of  two 
popular  forms  of  the  finite  difference  representation,  de¬ 
pending  on  whether  the  approximations  of  the  space  derivatives 
ere  set  to  provide  the  forward  difference  (or  explicit)  equa¬ 
tion  or  the  baokward  difference  (or  implicit)  equation.  The 
implicit  form  leads  to  the  requirement  to  solve  a  large  system 
of  simultaneous  linear  equations  (realistically,  500  or  more). 
The  explicit  form,  lacking  the  mathematical  complexities  of 
the  implicit  form,  can  be  solved  by  a  step-by-step  solution  of 
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explicit  relationships,  starting  with  given  initial  values 
of  u  in  the  region  it  for  a  given  time  t. 

The  explicit  form  has  one  basic  disadvantage,  however. 

In  moving  forward  in  time  from  the  time  of  given  initial 
conditions,  stability  considerations  require  that  a  relation¬ 
ship  between  grid  size  and  time  increments  be  maintained  such 
that  decreasing  the  size  of  the  grid  requires  a  corresponding 
decrease  in  time  step.  In  computations  over  long  periods  of 
time  where  data  requirements  for  specific  times  indicate 
fairly  large  increments  oould  be  used,  the  size  of  the  time 
step  is  restricted  by  the  useable  grid  size.  This  restric¬ 
tion  becomes  a  frustration  that  creates  a  vast  and  often 
prohibitive  amount  of  computation  to  obtain  a  solution.  The 
speed  and  capacity  of  modern  computers  has  permitted  wider 
application  of  the  explicit  method  to  problems  involving  a 
relatively  large  number  of  computational  steps,  but  there  is 
still  a  practical  limit  to  the  number  of  time  steps  to  be 
oomputed.  Sven  so,  the  solution  of  the  explicit  equation 
remains  a  prime  computer  method  of  obtaining  approximate  solu¬ 
tions  to  parabollo  partial  differential  equations  of  interest. 

The  solution  of  systems  of  linear  equations  has  been  a 
subject  of  intensive  study  by  mathematicians  for  many  years. 

As  a  result,  oertaln  systematic  iteration  procedures  suitable 
for  hand  calculation  have  been  available  to  solve  the  linear 
systems  of  the  implicit  method.  Yet,  the  additional  com¬ 
plexities  of  these  procedures  and  the  need  for  iteration  to 
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solution  made  the  trade-off  point  between  i;:.y.icit  and  ex¬ 
plicit  methods  heavily  in  favor  of  the  explicit  form  for 
the  class  of  problems  that  could  be  attacked  by  hand  numer¬ 
ical  calculation.  After  large  scale  digital  computers 
offered  potential  for  raptdily  solving  large  numbers  of 
simultaneous  equations,  interest  in  tne  implicit  form  was 
renewed.  UrunK  and  Nicholson  seem  to  nave  been  the  first 
to  use  implicit  methods.  In  1947,  tnese  men  demonstrated 
that  the  stability  restrictions  of  tne  explicit  method  do 
not  apply  to  the  implicit  form  (Hef  4:10£).  Ibis  promise 
of  freedom  of  choice  of  space  and  time  increments,  subject 
only  to  truncation  error  and  convergence  conei derations , 
offers  attractive  advantages  that  offset  the  disadvantage 
of  complexity  of  the  implicit  form. 

As  previously  mentioned,  a  number  of  routines  for  solu¬ 
tion  of  implicit  aquations  by  digital  computers  are  In  cur¬ 
rent  use.  Since  the  system  of  equations  is  generally  large, 
the  more  direct  method  of  matrix  inversion  is  not  practical; 
thus  most  routines  use  an  Iterative  procedure.  Of  these 
useful  routines,  the  ones  that  are  in  common  use  are  based 
on  the  metnod  of  successive  displacements— often  referred  to 
as  the  Oauss-Siedel  method.  Successive  displacements  is 
probably  tbe  simpliest  iterative  procedure,  is  easily  adapted 
to  computer  programs,  end  is  therefore,  widely  used  (tof 
10:374).  Routines  have  been  developed  to  speed  tne  rate  of 
convergence  of  the  method  of  successive  displacements. 
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In  his  doctoral  dissertation  in  1950,  0.  Young,  working  in 
two  dimensions,  extended  earlier  work  by  Prankel  and  pro¬ 
vided  a  theoretical  basis  for  a  process  that  increases  the 
rate  of  convergence  by  systematically  overrelaxing  the  Gauss- 
biedel  solution  of  boundary  value  problems  where  u  is  speci¬ 
fied  on  the  boundary.  Though  the  theoretical  basis  is 
somewhat  restrictive,  the  practical  success  of  overrelaxation 
has  led  to  the  application  of  the  method  to  other  problems 
with  a  significant  increase  in  the  rate  of  convergence  being 
realized  (hef  4:242).  specifically,  overrelaxation  has  been 
applied  to  the  problem  of  transient  heat  transfer.  In  1956, 

J.  tfegsteln  suggested  a  means  of  accelerating  the  conver¬ 
gence  of  iterative  solutions  of  problems  of  the  form 
(Hef  9).  This  procedure  was  empirically  adapted  to  the 
solution  of  the  implicit  representation  of  the  transient  heat 
transfer  equation  by  B.  Kaplan  and  N.  Clark  in  1956  (Hef  5). 

In  summary,  analytical  solutions  to  partial  differential 
equations  are  generally  not  available,  so  finite  difference 
methods  are  used  to  obtain  approximate  solutions.  Two  pop¬ 
ular  forms  of  the  finite  difference  representation  are  avail¬ 
able—  the  explloit  and  the  implicit  form.  The  expliolt  form 
may  be  solved  by  a  simple  step-by-step  procedure,  but  this 
form  suffers  from  restrictive  stability  considerations.  In 
fact,  stability  restrictions  make  the  explloit  form  unsuitable 
for  the  olass  of  problems  requiring  solution  over  large  time 
increments  (e.g.  the  determination  of  temperatures  in  a 
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nuclear  reactor  after  saut  down/.  Thus,  one  must  resort  to 
tne  implicit  form  for  a  solution,  respite  relative  freedom 
of  choice  of  grid  size  and  time  steps,  if  the  number  of 
iterations  required  to  obtain  a  solution  for  each  time  step 
is  too  large,  the  implicit  form  may  be  equally  useless.  Then 
tne  only  nope  of  obtaining  a  solution  on  a  given  type  of  com¬ 
puter  i3  to  speed  tne  rate  of  convergence  of  the  iterative 
solution,  thereby  reducing  the  number  of  iterations  required. 
Aside  from  maxing  some  solutions  possible,  a  successful 
metnod  of  speeding  convergence  of  a  slower  metnod  of  solving 
the  system  of  simultaneous  equations  also  results  in  attrac¬ 
tive  economy  of  computer  operation.  <»ith  these  potential 
gains  as  impetus,  a  major  effort  within  the  general  field  of 
research  in  numerical  analysis  is  being  devoted  to  finding 
rapidily  converging  solution  procedures. 

rurpose  ana  ocope 

inis  study  is  an  attempt  to  add  to  the  practical  Xnow- 
ledge  of  the  use  of  techniques  for  accelerating  the  iterative 
solution  of  the  implicit  form  of  the  finite  difference  repre¬ 
sentation  of  the  transient  heat  transfer  problem,  opecifi- 
cally,  the  study  provides  a  comparison  of  two  tecnniques 
for  accelerating  the  metnod  of  successive  displacements; 
successive  overrelaxation  and  an  adapted  legstein's  accel¬ 
erating  technique.  The  primary  points  of  comparison  art 
rate  of  convergence  to  solution,  associated  error,  and  simpli- 
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city  and  applicability  of  the  method. 

The  comparison  of  the  two  methods  is  provided  by  the 
solution  of  a  sample  heat  transfer  problem  of  conduction 
in  a  cubic  space  whioh  has  a  constant,  distributed  source 
of  heat  energy  in  the  interior  region  and  boundaries  that 
are  maintained  r.x  a  constant  temperature.  The  ezaot  problem 
chosen  has  an  analytical  solution  which  provides  the  stand¬ 
ard  for  measure  of  error  magnitudes. 

Two  additional  points  in  the  general  area  of  acceler¬ 
ating  the  rate  of  convergence  of  the  method  of  successive 
displacements  are  investigated. 

1.  Using  an  experimental  approaoh,  an  attempt 
is  made  to  determine  the  effeot  on  the  rate  of  convergence 
of  scanning  the  finite  difference  grid  in  a  spiral  mode 
from  the  fixed  boundary  conditions.  The  rate  of  convergence 
is  then  compared  to  the  rate  obtained  using  a  conventional 
serial  scanning  sequence. 

2.  The  theoretical  basis  for  the  accelerating 
effect  of  successive  overrelaxation  and  the  adapted  Wegatein 
technique  is  reviewed  for  similarities  and  differences  of 
the  two  methods.  This  information  is  then  used  as  a  basis 
for  experimental  attempts  to  alter  the  basio  adapted  Wegstein 
technique  to  improve  the  speed  of  convergence  of  this  method. 
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XI.  Theory 

The  Partial  Differential  ociuation  of  Transient  Heat  Transfer 

From  tne  theory  of  heat  transfer,  tne  partial  differ¬ 
ential  equation  of  neat  conduction  with  an  internal  heat 
source  or  sin*  may  be  expressed  as  the  parabolic  partial 
differential  equation  (nef  1:45) 

sft=KV’KVa  +  > ki*'****  U) 

where  u(x,y,z,t)  is  tne  temperature  (°C) 

a.  la  the  conductivity  of  the  medium  ( cal/cm-sec-°C ) 
q  is  the  strength  or  tuo  source  or  sink  (oositive 
for  a  source — negative  for  a  sink)  in  units  of 
energy  generated  or  absorbed  per  unit  volume  time 
t  cal/onP-sec ) 

a  is  the  diffusivity  ( cm*/ sec )  waich  may  be  deter¬ 
mined  from  Cl=  K/<pCp  with  as  density  (gr/cm*  ) 
andCpas  the  heat  capacity  (cal/gr-*C) 
y  is  the  del  operator  of  vector  analysis  which  is 
expressed  as 

with  i,  t,  and  Tt  tne  usual  cartesian  unit  vectors 
>Vith  the  simplifying  assumption  that  the  conductivity  is  a 
constant,  equation  1  reduces  to 

iff  =  + 

fc. 
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&£  Finite  gl££&aflfie  ^presentation 

The  finite  difference  representation  of  equation  3  that 
is  used  in  this  study  is  obtained  by  replaoing  the  first 
partial  derivative  with  respect  to  time  by  the  approximation 

QOz  -  -  U( x, 4,.  jL)  (4) 

dt 

and  the  seoond  partial  with  respeot  to  a  given  space  variable 
by  an  approximation  of  the  form 


~  ~  2 17 (l,  f,A»,  i-fat)  (5) 

4* 


wH«re  4h  :  Ax  l •) 

to  obtain  the  lmpllelt  fora,  or  by 

i'u.  .  V(i.i".*,4..t)-ivc>L.,.-L,iu-u(iri.i.<.t)  ... 

3x*  -f,i 

to  obtain  the  explicit  fora,  riow  these  approx last ions  are 
foraed  using  truncated  Taylor  series  expansions  for  a  sample 
elliptic  and  parabollo  partial  differential  equation  and  the 
resulting  aatrix  forms  for  a  sample  serial  scan  are  shown  in 
detail  in  Appendix  A. 

•substituting  the  approximations  given  in  equations  4  and 
5  and  sialllar  expressions  for  the  partial  derivatives  with 
respeot  to  the  y  and  z  space  variables  into  equation  3,  one 
obtains  the  lmpllolt  fora  of  the  finite  difference  equation 
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OAJt  “  'nl 

+ 

+■  — 

/rota  this  equation,  one  can  aee  that  if  the  temperatures 
over  the  region  are  known  at  tine  t  and  the  boundary  values 
and  source  strength  q  are  fixed  for  any  time,  then  the  temp¬ 
eratures  at  the  interior  nodes  at  time  t  +  At  form  a  system 
of  N  simultaneous  equations  in  ft  unknowns;  one  unknown  for 
each  interior  node  point.  In  matrix  notation,  this  system 
may  be  expressed  as 

%X}=t  (9) 

Alternately,  one  may  use  equation  7  to  approximate  the 
seoond  partial  to  obtain  the  explicit  form 

&")  _  U(l+ 1, ■%,"&,£) >4  TJ (lr<s$, &,  Jt) 

Q.ut  ‘  15 

4-  ^  (*■>  1  +  •&.  t)  +  U  (L>4- 

4% 

l  Ufi.-a.-X+i.jtMWt.a 

W  - 

uo) 

Here,  eacn  interior  node  temperature  at  t-fAt  is  expressed 
explicitly  in  terms  of  known  temperatures  at  time  t  so  that 
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there  is  no  need  to  solve  a  system  of  simultaneous  equations. 

■solution  procedures  for  systems  of  linear  Equations 

It  has  been  shown  that  the  implicit  form  of  the  para* 
bolic  equation  *  leads  to  a  system  of  simultaneous  finite 
difference  equations.  The  following  forms  of  the  solution 
of  the  system  will  be  considered: 

1.  The  method  of  successive  displacements  (often 
referred  to  as  the  Gauss-sledel  metnod) 

2.  successive  overrelaxation 

3.  Adapted  •isgstein  convergence  technique 

successive  Displacements 


solution  Procedure 


For  simplicity  in  the  expression  of  relationships,  the 

(11) 


following  notation  is  adopted 

(o-O 

UU, 

where  n  =  the  number  of  the  iteration 


Consider  the  two  dimensional  parabolic  equation 

VHtr  4- 

CL  dt 


(12) 


in  the  ioplioit  difference  form  with  the  gird  size  equal 
in  the  x  and  y  direction  (l.e.  ^  -  “&»  )  or 


*  k  system  of  simultaneous  equations  also  results  from  the 
finite  equations  of  elliptic  boundary  value  problems  •-  see 
Appendix  A. 
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+U(  x, 1+1, *.+*£)  +  U^L,  *-»,*+*£)]  . 

+ 

where 

Using  the  method  of  successive  displacements  to  solve  tne 
system  of  equations  resulting  from  the  above  sample  equation 
involves  the  following  steps* 

1.  estimating  a  starting  value  of  U  for  the  in¬ 
terior  grid  points  li.e.  those  points  in  the  region  consid¬ 
ered,  excluding  the  known  boundary  points) 

2.  Improving  these  initial  estimates  according  to 
an  arbitrary  but  fixed  ordering  of  points. 

3.  Using  improved  values  as  soon  as  available 


*  There  is  a  similiar  procedure  called  simultaneous  dis¬ 
placements  which  is  generally  not  as  rapidily  converging  as 
successive  displacements,  though  it  may  converge  in  cases 
where  successive  displacements  diverges  (Ref  3:133).  The 
essential  difference  is  that  in  simultaneous  displacements 
no  improved  value  of  U  is  used  until  all  values  are  improved, 
oiace  convergence  is  generally  slower,  simultaneous  displace¬ 
ments  ia  not  u  popular  method  (oef  4:826). 


ad  ft 

if 
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4.  Continuing  the  iterative  cycle  until  the  abso¬ 
lute  value  of  U  of  the  present  iteration  minus  the  value  of 
U  of  the  previous  iteration  is  equal  to  or  less  than  an . 
established  criteria  put  into  the  computer  as  input  data. 

Thus,  if  the  iteration  notation  is  added  to  equation  13 
the  result  is 

U(t, ■$,£***)  -  ,+4/t 

+  £+&t)  +U(jL,f-i,£’+AJt)J 

+  ]  1161 

Convergence 


Probably  the  slmpliest  statement  of  the  necessary  con¬ 
ditions  for  convergence  of  the  method  of  successive  displace¬ 
ments  is  given  by  Forsythe  and  »asow.  They  state  tnat,  "If 
(the  matrix)  *  (of  equation  9)  has  diagonal  dominance  *  and  is 
not  reducible,  then  the  method  of  successive  displacements 
converges"  (Kef  4:236).  In  practioe,  decreasing  the  grid 
size  has  the  effect  of  decreasing  the  diagonal  dominance  by 
increasing  the  sum  of  the  off-diagonal  terms,  because  the 
diagonal  is  weakened,  the  convergence  is  slower. 

This  simple  statement  of  convergence  is,  however,  more 
restrictive  than  necessary,  a  less  restrictive  but  more 

*  See  "Definitions",  Appendix  C 
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complex  criteria  may  be  developed  as  follows  (f>ef  4:209). 
Partition  tne  matrix  A  into  three  matrices  such  that 

A -  (e  +  5  +  f  )  uci 


where  z.  is  the  lower  triangular  matrix  of  Rjj  where  i  >J 
am  there  are  zeros  elsewhere 

F  is  the  upper  triangular  matrix  of  a  where  i  <  J 

i  J 

and  there  are  zeros  elsewhere 

u  has  the  diagonal  elements  a^j  and  zeros  elsewhere 
It  is  tnen  possible  to  write  the  metnod  of  successive  dis¬ 
placements  in  tne  following  matrix  form  (hjf  4:2.56) 


117) 


fhea  tne  necessary  and  sufficient  condition  for  convergence 
is  that  all  the  eigenvalues  of  the  matrix  -(d+4  i'F  are  less 
tnan  cue  in  modulus,  further,  thc-se  eigenvalues  1 i )  are 
the  zeros  of  the  daterminantal  equation 

del  (v[l  +  /?5  +F  )-  0  lie) 

examples  of  the  forms  of  these  matrices  are  given  in  Appen¬ 
dix  B. 

oince  the  eigenvalues  of  the  large  matrices  encountered 
in  practical  problems  are  not  easily  determined,  this  cri¬ 
terion  is  not  very  useful  for  estimating  whether  a  given 
system  will  converge  for  the  metnod  of  successive  displace¬ 
ments.  The  size  of  the  maximum  eigenvalue  does,  however. 
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play  an  Important  role  in  accelerating  the  convergence  of 
the  metnod  of  successive  displacements.  This  role  is  dis¬ 
cussed  later. 

one  more  important  point  must  be  made.  The  method  of 
successive  displacements  depends  critically  on  the  order  in 
which  the  various  unicnowns  are  oomputed,  sinoe  the  size  of 
the  eigenvalues  depend  on  this  order  and  the  smaller  the 
maximum  eigenvalue,  the  faster  the  process  converges  (Ref 
4:218,257) . 

successive  Overrelaxation 
solution  Procedure 

successive  overrelaxation  is  the  first  method  consid¬ 
ered  of  «he  two  methods  of  accelerating  the  convergence  of 
successive  displacements  that  are  compared  in  this  study. 

In  this  method  acceleration  is  achieved  by  a  simple  modifi¬ 
cation  of  the  equation  for  solution  by  successive  displace¬ 
ments  (*ef  9:386).  i'he  modified  form  for  the  transient 
heat  equation  in  three  dimensions  is 

/«+i\  f  rt  T 

*■  +  &+At) 

+  WPS 

where  <a)  =  relaxation  faotor 

The  procedure  for  solution  is  the  same  as  for  the  method  of 
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successive  displacements.  Note  that  when(4>»l,  the  form 
reduces  to  the  method  of  successive  displacements. 

Applicability  of  tne  Method 

Young  in  his  doctoral  dissertation  showea  that  succes¬ 
sive  overrelaxation  improves  the  convergence  rate  of  the 
method  of  successive  displacements  for  a  class  of  matrices 
with  a  property  he  calls  Property  (A;  *  (hef  4: £43).  No 
attempt  will  be  made  here  to  describe  in  detail  Young's 
proof  of  the  theoretical  basis  for  the  matnod  of  successive 
ovarrelaxation.  hather,  interest  in  tne  theory  benina  the 
metaod  will  be  centered  on  the  determination  of  the  optimum 
relaxation  factor.  A  detailed  description  of  the  theoreti¬ 
cal  background  may  be  found  in  Forsythe  and  <*asow  beginning 
on  page  242  (Hef  4).  specifically,  the  method  is  applicable 
to  the  approximations  of  the  parabolic  partial  differential 
equation  of  heat  transfer  used  in  this  study  (nef  4:105). 

It  is  of  interest  to  note  that  overrelaxation  has  been 
so  successful  that  the  method  has  been  profitably  applied 
to  problems  with  matrices  that  do  not  have  rroperty  (A), 
inis  success  nas  led  to  work  toward  the  extension  of  Young's 
proof  to  a  more  general  class  of  problems  and  has  provided 
a  basis  for  the  speculation  that  the  method  may  be  useful 
for  any  problem  (hef  4:260,261). 

*  oee  "Definitions"  Appendix  C 
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determining  the  optimum  nelaxatlon  Factor 

*hen  the  optimum  relaxation  factor  is  known,  signifi¬ 
cant  savings  in  iterations  can  be  accrued  by  using  overre¬ 
laxation  in  the  solution  of  partial  differential  equations. 
There  is  unfortunately,  no  simple  direct  way  of  accurately 
determining  the  optimum  factor  prior  to  the  start  of  com¬ 
putations  except  for  simple  problems  involving  simple 
geometry.  Compounding  the  problem  is  the  rate  of  degrad¬ 
ation  of  savings  as  the 
cnosen  relaxation  faotor 
moves  away  from  the  opti¬ 
mum  value.  A  character¬ 
istic  plot  of  relaxation 
factor  versus  number  of 
iterations  to  solution  is 
shown  in  figure  1.  The 
result  is  that  determin¬ 
ing  the  optimum  faotor  is 
the  most  important  prob¬ 
lem  in  using  the  method  of 
successive  overrelaxation 
(Ref  4:257).  A  plot  such 
as  shown  in  flgurw  1  is  a  way  of  determining  the  optimum 
faotor,  though,  of  oourse,  a  very  impractical  way  unless  a 
large  number  of  elmlllar  problems  are  to  be  considered. 


Fif.  1 

Snplt  not  of  SHixitlon  Factor 
Versus  Iteration*  to  Solution 
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figure  1  demonstrates  additional  important-  properties 
of  the  relaxation  faotor.  First,  the  limits  on  tile  optimum 
factor  arei<(a)<2  where  (&  s  I  is  equivalent  to  the  .method 
of  successive  displacements,  factors*  less  than  one  are  term¬ 
ed  underreiaxation  and  are  known  tc  be  less  profitable  than 
U)  >  i  for  acceleration  i  no  f  4:3c8).  The  upper  limit  is 
based  on  the  maximum  eigenvalue  of  the  method  of  successive 
displacements,  This  relationship  will  be  described  lr.  more 
detail  later.  The  second  property  is  related  to  the  shape 
of  the  curve  in  the  vicinity  of  the  optimum  faotor.  The 
smaller  rate  of  change  of  slope  on  the  nigh  side  of  the 
optimum  demonstrates  tnat  it  is  better  to  overestimate  the 
relaxation  factor  than  to  underestimate  it.  A  theoretical 
basis  for  this  property  has  been  developed  (^ef  4:857). 

Various  means  of  estimating  the  optimum  fector  have 
been  proposed,  as  a  starting  point  in  reviewing  the  more 
prominent  metnod’s,  consider  the  means  for  determining  the 
factor  for  ^place’s  equation  over  a  regular  space.  For  the 
case  cf  equal  grid  size  (h)  and  the  finite  difference  repre¬ 
sentation  given  by  equation  74,  the  optimum  relaxation 
factor  may  be  found  from  (Ref  9:366) 

(jO  In  -  /  +  - -  ■ 

( i  +  rr~-\-)z  u**) 


j.6. 


where 
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(19b) 


(19c) 


The  variable  X  is  referred  to  a s  "the  spectral  radius  *  of 
tbe  linear  transformation  defined  by  the  Gauss  oiedel  (i.e. 
successive  displacements)  method"  (Kef  9:387).  Tbe  rela¬ 
tionship  of  A  to  tbe  progress  of  tbe  solution  by  successive 

displacements  may  be  demonstrated  as  follows.  Define 

^rm.0  =  u  <<*♦•>  _uf(n,  (1M, 

Then  X  may  be  described  as  tbe  limiting  ratio  of 
d<m)  or  (Ref  9:387) 


X  = 

na—woo 


cLM 


(19e) 


Young  goes  on  to  give  tbe  following  formula  for  computing 
X  for  Laplace's  equation  for  a  rectangle  with  sides 
0.-R-&  and  bsS’fc  where  H  and  <d  are  integers  (Ref  9:389) 

X  =  [l  (cos  i  +  Cos U9f) 

He  then  advises  that  for  regions  other  than  a  rectangle, 
may  be  estimated  for  a  rectangle  containing  tbe  region  under 
consideration. 

For  tbe  parabolic  equation  of  tbe  form 


*  8ee  "Definitions",  Appendix 
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r  H9gi 

3  kx 

ana  the  finite  aifferencs  approximation  of 

•u  A-[ua+i,h  tua-i.ii  +  ua+<,  ^*.1 

+  +[i-fl.]u(A,«  u“h' 

Young  provides  the  following  utfrr  bound  for  A  (xtef  9:403) 


where 


A. 

I  4  A. 


2 


( 191 ) 

(19J) 


Then  for  tne  two-diment ionel  par.itoiic  equation 

in  =  tUi  + 

3xl  3^ 

the  upper  bouna  for  A  becomes  l  r.e f  y;40b) 


1 19k) 


(191) 


of  course,  none  of  the  above  formulae  for  A  provide  a  pre¬ 
cise  expression  for  determining  ah  optimum  relaxation  factor 
for  the  parabolic  equation  oi  transient  neat  transfer.  They 
do,  however,  provide  an  insight  into  the  nature  of  the  de¬ 
pendence  of  &)  on  the  parameters  or  the  problem. 

wince  there  is  no  direct  method  of  computing  the  optimum 
relaxation  factor  prior  to  the  start  ox’  the  iterative  pro¬ 
cess,  a  numoer  of  possibilities  for  estimating  the  factor 
have  been  suggested,  principal  interest  in  this  study  is  in 
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tne  process  of  estimating  the  factor  by  running  the  computer 
program  withU>fcl  ( i.e.  by  successive  displacements)  (Hef  4: 
368,369).  now  this  leaas  to  the  optimum  factor  is  described 
as  follows.  One  may  define  the  optimum  factor  by 


(*>opt  r 


(20) 


i 

where  is  the  dominant  eigenvalue  of  tne  matrix  -(D+E)  ^ 

0 

of  the  method  of  successive  displacements 
Under  the  assumption  that  the  dominant  eigenvalue  l^is  real, 
a  condition  that  will  exist  for  example  when  the  a.  matrix 
is  symmetric  with  all  CL^>0  —  a  common  occurrence  with 
finite  difference  equations  (nef  4:252)  — the  value  of  r[‘^ 
may  be  determined  from  the  limiting  ratio  of  the  norms  of 
the  residuals*  of  two  successive  iterations,  or  (nef  4:369) 


no¬ 


where 


II  Y~|| 

(21) 

j(r-  ufr)2]Vi 

(22) 

is  r  the  number  of  unknowns  in  the  system 
This  process  is  subsequently  referreu  w  as  tue  norm  of  the 
residuals.  Alternately,  one  may  use  the  first  power  norm 

£  _  li  Y m-hll(t) 

'  II  Y«v  ||(o 

where 

*8ee  ’'definitions’1.  Appendix  o 


(23) 
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JL» 


(24) 


iL/vn,  Qm,  21  *?  ff) 

/VL-*  oO  (25) 

The  apparent  disadvantage  of  running  for  a  number  of  iter¬ 
ations  with  (4>sl  is  that  one  loses  the  effect  of  acceleration 
•«hile<a>is  one.  Unless  there  is  a  rapid  convergence  to  the 
vicinity  of  this  loss  of  acceleration  may  degrade  the 

metnod  to  a  point  where  very  little  acceleration  is  achieveu. 

Is  it  possible  then  to  determine  tne  optimum  factor 
waile  running  with  an  estimate?  To  oegin  this  discussion, 
the  effect  of  the  successive  overrelaxation  modification  of 
tne  successive  displacements  equation  on  the  pertinent  eigen¬ 
values  is  of  interest.  The  eigenvalues  of  the  overrelaxation 
equation  are  given  by  the  expression  («ef  4:247) 

daX  [(<V^3  *  E)^  +  (l- )  5  +  P]  =  0  (26) 

where  o,  t,  and  F  are  as  given  on  page  14 
This  compares  to  equation  lb,  the  equation  for  finding  the 
eigenvalues  of  the  method  of  successive  displacements.  If 
one  defines 

H(ou)  =  -  ((ii'b  +  E)  [f +  0-<b)~')d]  (g7) 

then  two  important  properties  of  H(<£)  can  be  cited.  "The 
(u)  * 

eigenvalues  r^of  iHOOj  remain  the  same  for  any  consistent 
order 9*"  («ef  4:251).  Order  nere  Implies  the  order  in  itiiich 
the  points  are  taken.  A  consistent  order  is  described  later. 
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Thus,  "when  the  matrix  A  has  Property  (A),  the  value  of 
is  independent  of  the  orderO'for  a  considerable  class  of 
orders  (called  consistent)"  («ef  4:261). 

The  problem  encountered  in  running  with  a  trial  relax- 
ation  factor  to  determine  the  optimum  is  the  "sequence  of 
values  of  IlYmll  in  any  norm  is  a  rather  irregular  sequenoe  for 
(a)  near  OUopt"  (hef  4:370),  and  there  is  difficulty  in  de¬ 
termining  which  trial  value  of  (ft)  is  the  best,  experiments 
by  Ortega  have  indicated  that  the  sequence  of  values  of  the 
ratios  of  the  first  power  norm  ^demonstrate  differences 
between  values  in  the  sequence  that  are  low  in  magnitude 
for  (ftXltiptand  relatively  high  in  magnitude  for  trial  near 
or  slightly  larger  than  bfeptl^ef  4:371).  The  results  thus 
far  have  been  Inconclusive,  however. 

The  remaining  point  to  be  covered  in  this  subparagraph 
is  the  definition  of  consistent  order.  According  to  Young, 
an  order  is  consistent  "if  and  only  if  each  elementary  square 
mesh  of  the  net  is  bounded  by  four  arrows  with  zero  circula¬ 
tion"  (iief  4:245,246).  figure  2  demonstrates  consistent  and 
nonconsistent  orders. 

A  pertinent  question  is  wnat  happens  when  one  uses  an 
nonconsistent  order  for  successive  overrelaxation?  Varga  has 
proven  that  for  LaPlaoe's  equation,  no  nonoonsistent  order 
has  a  rate  of  convergence  as  great  as  that  oommon  to  all  con¬ 
sistent  orders  (tief  4:259).  Powers  Investigated  this  point 
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for  the  birionlet  difference  equation  *  for  a  rectangle  and 
some  of  her  results  have  been  displayed  in  Table  1.  (Ref  7) 
These  results  are  generally  considered  inconclusive  since 
she  ohose  to  use  a  measure  of  convergence  criteria  of 


1,Jl|U(Aa,-Uu.4iJ  i  2 

which  proved  to  be  relatively  insensitive  (Ref  4:259). 
Convergence 

The  question  of  convergence  is  simply  covered  by  the 
following  quotation.  ttThe  power  of  Young's  method  lies  in 
its  acceleration  of  the  rate  of  convergence  of  an  already 
convergent  process,  not  in  any  ability  to  create  a  conver¬ 
gent  algorithm...."  (hef  4:254). 

The  Adapted  legate in  Convergence  Technique 


The  legate in  Technique 

The  true  legatein  technique,  which  is  the  basis  for 
this  second  method  of  accelerating  the  method  of  successive 
displacements,  was  originally  developed  for  accelerating 
the  iterative  solution  of  a  single  equatlon-not  a  system 
of  linear  equations.  If  the  problem  is  to  determine  the 
solution  of 

F(x):  O  (28) 

one  may  often  express  the  equation  in  the  form  (Ref  8:9) 

X  =  Six)  (29) 


See  "Definitions,"  Appendix  C 
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Table  i 


Sample  Results  of  Solution  of  the  Dirichlet 
Difference  Equation  by 
Consistent  and  Fonconsi stent  Orders 


Number  of  Iterations 

Order* 

Using  Optimum  <i3 

With  U)  *  1 

Consistent 

Serial b 

34 

146 

Even-Odd c 

30 

146 

Fonconsistent 

Ordered  by  rows, 
direction  reversed 
fov  each  row 

43 

148 

Spiral* 

41 

147 

Random  selection 

33 

147 

Spiral  border, 
random  interior 

36 

147 

This  data  is  extracted  from  reference 

*  The  problem  was  14  nodes  along  the  y  axis  and  6  nodes 
along  x 

b  This  scan  sequence  was:  U(l,l)  to  U(l,14),  then  U(2,l) 
to  U(2,14)»  and  so  on  to  U(6,l)  to  U(6,14) 

c Starting  with  U(l,l),  the  points  for  which  i+j  is  even 
were  computed  first  in  sequence*  Then  returning  to  U(2,l) 
the  remaining  points  we  re  computed  in  sequence. 

^Udjj)  was  taker,  for  j  «  1,...,14»  then  U(i,14)  i  e  2,3*4, 5 
then  U(6,j)  for  j  *  14,  13t»»»»  2,  1$  then  U(i,l)  for 
i  «  5,4,3, 2.  This  pattern  was  ther.  repeated  starting  with 
U(2,2)  to  U(2,13)  etc. 


26 


GA/Phys/63-6 


Then  an  iterative  solution  can  assume  the  following  algorithm 


(30) 


Alternately,  if  F0O=o  oannot  he  written  in  the  form  given 
by  equation  29,  the  algorithm 


-  x(m)+  rF(xn] 


(31) 


where  P ±0  is  some  appropriately  ehosen  constant 
can  be  used. 

Assuming  FM-O  has  a  solution,  the  normal  iterative 
process  is  to  assume  a  value  of  use  this  value  in 

equation  30  lor  31)  and  solve  for  Then  is  used  to 

find  X^and  so  on  until  the  absolute  value  of  the  residual 

(i.e.  )  is  less  than  some  specified  criterion. 

The  sequence  of  values  of  /'"will  show  any  one  of  the  follow¬ 
ing  characteristics  (*ef  8:9) 

1.  Oscillate  and  oonverge 

2.  Oscillate  and  diverge 

3.  Converge  monotonioally 

4.  Diverge  monotonlcally 

To  obtain  the  accelerating  technique,  the  basic  iter¬ 
ative  equation  is  modified  to  give  (Kef  8:9) 

0.g)x(m*0 

where  q  is  the  accelerating  faotor 

the  tilde  )  indicates  a  value  computed  by  the 
accelerating  technique 


(32) 


87 


GA/Phys/63-8 


tfegstein  states  that  with  the  appropriate  selection  of  "q", 
one  may  cause  convergent  cases  to  converge  rare  rapidly  and 
divergent  cases  to  beoome  convergent. 

The  procedure  used  is  best  described  by  a  diagram  of 
the  process  (Kef  8:11).  This  diagram  is  provided  in  figure  3. 
Starting  with  the  estimate  ,  compute  an  Improved 

value  using  x- 5 [£<'**].  On  this  first  iteration,  the 
modified  form  (equation  32)  is  not  used  and  the  values  for 
the  next  iteration  are  set  up  thru  the  center  path  on  the 
diagram.  On  subsequent  iterations,  the  process  loops  thru 
the  outside  loop  until  the  residuals  are  as  small  as  desired. 

The  derivation  of  "q"  may  be  described  in  the  following 
manner  (Kef  8:9-11).  "Geometrically,  the  solution  of 
amounts  to  the  problem  of  finding  the  point  of  intersection 
P  (see  figure  4)  of  the  curves  X  and  The 

Iteration  can  be  represented  graphically  as 

follows.  Pass  a  vertical  line  through  a  point  on 

s  X  bo  that  it  intersects  the  curve  at  some  point 

r  (m-H  ('tn-chl 

A  with  the  coordinates  IX  ?x  J  ." 

"The  ideal  location  for  x^'^of 
on  AB  would,  of  course,  be  the  intersection  point  C  with  the 
normal  to  AB  drawn  through  P.  Thus,  'q'  should  be  chosen 
suoh  that 
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Set: 


^(«*»)_  x^» 

.as  small  as 

\de*inw^/ 


SalVC 

CL  - 

xCna+*)  _  ^(<"0 

7  _  ^o*»-o 

fr  = 

a 

a-/ 

— 
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A  Graphical  Representation  of  the  Wegstein 
Accelerating  ?aotor 
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To  determine  *q*  approximately,  observe  that  PC  *  BC  and 

-a  =  PC  / /fc  (34) 

where  'a'  is  a  value  of  f’(x)  between  -a  and  P.  Thus, 


-a  - 


i 


2- 


or 


-8 


8  = 


a-/ 


(35) 


136) 


Since  a  more  convenient  expression  is  lacking,  'a'  can  be 
approximated  by  a  suitable  difference  quotient** 

a  «.  iM-ik**-'11]  _  x*"*0  - 

U  **  x(tn)  _  "  x  (^n)  _  l37* 

A  variation  of  the  above  is  to  compute  a  new  value  of 

"q"  after  each  iteration  using 

_  vtfn)  ,  , 

/u  ^  ^  (38) 

~  ^(mj  *(m-0  V  ; 

Adaptation  to  Simultaneous  equations 
Kaplan  and  Clark  used  the  basic  Wegsteln  technique  and 
through  a  series  of  experiments  empirically  adapted  it  to 
the  solution  of  the  system  of  equations  resulting  from  the 
lmplioit  finite  difference  representation  of  the  equation 
of  transient  heat  transfer.  The  forms  of  the  equations  used 
are  identical  to  equations  32,  36,  and  38 
!(<n40 


where 


%  - 


CL 


a- 1 


(39) 


(40) 
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Ul) 

The  basic  solution  procedure  used  is  the  method  of  successive 
displacements  (nef  6:80). 

The  steps  of  the  adapted  process  is  presented  in  simpli¬ 
fied  diagrammatic  form  in  figure  5  (Ref  5:10).  On  the  first 
iteration,  the  initial  values  are  set  up  but  no  acceleration 
is  attempted.  This  process  of  setting  up  but  not  accelerating 
the  solution  is  continued  until  the  solution  proceeds  to  a 
specific  number  of  iterations.  This  number  is  predetermined 
prior  to  running  tne  program  and  is  input  data  to  the  computer. 
When  the  appropriate  iteration  for  the  first  application  of 
the  technique  has  been  reached,  a  slope  "a*  is  computed  for 
each  node,  tested  to  insure  that  it  is  less  than  one  and 
greater  than  zero,  and  if  the  criteria  is  met,  an  acceler¬ 
ating  faotor  "q"  is  computed  (Ref  6:80).  This  value  of  "qn 
that  has  been  computed  for  eaoh  node  is  further  tested  to  in¬ 
sure  that  it  does  not  exceed  a  given  maximum  value  and  if  the 
value  exceeds  the  maximum,  "o.”  is  set  equal  to  the  maximum 
value,  ohould  "q"  not  exceed  the  maximum,  the  value  computed 
from  the  slope  is  used.  Then  the  value  of  "o"  resulting  from 
this  test  is  used  to  accelerate  the  solution  using  equation  39. 
The  accelerated  value  of  temperature  then  replaces  the  value 
computed  by  the  method  of  suooessive  displacements  in  the 
computer  memory  ana  the  process  proceeds  to  the  next  node, 
ohould  the  slope  test  fall,  the  accelerating  technique  is  not 
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Pig.  5 


Block  Diagram  of  the  Adapted  Vegateia  Teohnltpie 
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applied  and  the  value  of  temperature  computed  by  auoceesive 
displacements  is  used  in  subsequent  computations.  Follow¬ 
ing  this  initial  application  of  the  technique  to  selected  nodes, 
the  pertinent  tests  and  applications  of  the  accelerating 
factor  are  only  made  at  a  repetitious  interval  of  iterations 
specified  by  an  input  to  the  oomputer  (e.g.  once  every  five 
iterations).  For  those  iterations  wnen  the  technique  is  not 
applied,  the  value  or  temperature  oomputed  by  the  method  of 
successive  displacements  is  used  as  before  the  first  applica¬ 
tion. 

In  the  prooess  of  refining  this  adaptation  of  the  ftegstein 
technique,  Kaplan  and  Clark  made  the  following  observations 
on  the  size  of  pertinent  parameters. 

1.  The  correct  iteration  to  apply  the  technique 
is  when  the  value  of  the  slope  "a"  is  greater  than  zero  or 
less  than  one. 

2.  If  the  technique  is  applied  too  soon,  the  solu- 
tion  may  actually  be  slowed  down.  Kaplan  and  Clark  provide 
the  following  criterion  for  the  first  application  of  the 
technique.  The  initial  application  should  be  made  on  the 
iteration  which  is  equal  in  number  to  the  depth  of  the  deepest 
node  from  the  boundary  (Aef  5:6).  For  example,  a  one  dimen¬ 
sional  problem  of  30  nodes  has  the  deepest  node  15  from  the 
boundary,  so  the  initial  application  of  the  accelerating 
technique  le  made  on  iteration  number  15. 
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3.  The  solution  may  also  be  slowed  down  if  the  tech¬ 
nique  is  applied  too  often.  This  result  is  attributed  to  the 
fact  that  "the  slopes  of  all  the  functions  do  not  settle  down 
immediately"  (Kef  5:2).  The  criterion  suggested  for  frequency 
of  application  of  the  technique  is  given  as  one  half  the  dis¬ 
tance  to  the  deepest  node  from  the  boundary  l  Kef  5:6).  Thus, 
for  the  example  in  paragraph  2  above,  the  technique  would  be 
applied  every  seventh  iteration. 


a  Comparison  of  quo oe salve  Overrelaxation  and  the  Adapted 
legate in  Technique 


an  apparent  similiarity  between  the  method  of  successive 
displacements  ana  the  adapted  Aegstein  technique  can  be  demon¬ 
strated.  This  relationship  may  be  shown  as  follows, 
tfor  convenience  let 


T1  -r  (7*14  0 

u«d  -  ^ c*> 

as  computed  by  ths  method  ol  successive  displacements.  Then 
the  equation  used  for  overrelaxation  may  be  expressed  as 


U (xjtjll,  Jfc+a£)  =  fc>[l)so  1  +  ( 43) 
The  equation  for  the  adapted  legsteln  technique  may  be  re¬ 
written  in  the  form 

Assuming  a  correspondence  between  and 
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V  Cx.^,1  ,£  +  4t)  ,  it  can  be  seen  that 

=  I-  g  (45) 

oince,  furthermore,  it  is  known  that  (m)  is  restricted  to 
values  between  1  and  2,  it  would  follow  that  for  the  rela¬ 
tionship  between  ”qH  and  CO  to  hold,  the  range  of  values  of 
*q"  would  be 

"l<£<0  (461 

ihe  fact  that  q  is  negative  is  consistent  with  the  findings 
of  Kaplan  and  Clark  (Kef  5:5).  Kaplan  was  not  as  restric¬ 
tive  on  the  maximum  negative  value  of  "q",  however,  de 
presents  one  series  of  runs  where  "r"  was  restricted  to 
values  between  -ICC  and  0. 

further  evidence  of  a  potential  tie-in  between  the  two 
methods  of  acceleration  i?>  found  in  the  discission  of  esti¬ 
mating  the  relaxation  factor  by  running  with  (a)  -  I  to 
obtain  an  estimate  of  the  maximum  eigenvalue  .  *s  pre¬ 

viously  shown,  the  optimum  relaxation  factor  may  be  found 
from  the  equation 


(at)  - 


Jo 


/*  o-  ;'/4 


(20) 


||.o 

Yz 

.  (47) 

r  i 

+4 

r-N 

w  ' 

Now,  since  from  equation  20  abovs,  the  limits  on 


.<«) 

'll  can  bs  shown  to  b# 
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0 )‘/z  £  / 

li  ' 

or  since  0"^''}2ia  positive, 

O  6(  l-%)il 

I  <  iO 

I  J  ^  i  O 

By  comparison,  the  equation  for  tne  slope  test  of  the  adapted 


(48) 


vfegsteic  technique  is  given  by 


_  u  <**> 

^ ^  (*n)  rr  0 


(41) 


The  s iiit.il lari ty  between  tnis  •  *-.pc  equation  and  the  equation 
for  rf? is  obvious.  An  additional  simlliarlty  is  found  in 
the  fact  that  Kaplan  experimentally  determined  that  the  slope 
should  be  restricted  to  values  between  0  and  1,  which  is 
compatible  with  the  restriction  on  1^. 

Though  there  is  a  definite  simlliarlty  in  the  way  ^ and 
• a*  are  oomputed,  a  major  difference  lies  in  the  method  of 
application  of  the  factor  computed  from  "a"  and  In  the 

adapted  eegstein  technique  an  accelerating  factor  ”q"  is  oom¬ 
puted  separately  for  each  grid  point,  while  for  overrelaxa- 
tlon  the  same  relaxation  factor  is  used  for  all  grid  points. 

Further  evidence  that  the  two  methods  are  olosely  re¬ 
lated  can  be  derived  using  an  expression  for  accelerating  the 
convergence  of  the  method  of  successive  displacements  given 
by  V.  N.  Faddeeva  (Kef  3:241).  her  expression  is 


/  ~  fA./ 


(49) 
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where 


i  le  a  vector 


fX,  is  determinea  from  the  ratio  of  the  components 
of  tne  vectors  -  y^'and  X^r',)  -  y(,n"0 

oir.ce  tne  notation  is  in  vectors,  consider  an  n  component 
vector  for  X. 


(5G) 


Xi 
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** 

xf“ 

x3 
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xf. 
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1 

xi** 

yf> 
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r 

X>» 
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• 

Equating  the  first  component 


v  -  y^1 
x ,  *« 


!-H"  [ 


Om+t)  .On) 

X»  - 


■] 


rationalizing  and  rearranging 

„  .  +  *?""*  -*.** 

A  |  ^  —  SI  ■  -  - 

»  -  /X, 


(51) 


Now  1st 


(-(V  = 


M-* 


You 


then 


-^r 


fX'i  _ ^  _  J*Lt  ~/4. 1  +!  _ 

/4,  -  / 


to  = 


f 


/-  jJ,, 


1 52) 


(53) 


(54) 
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Substituting  these  expressions  into  equation  51,  one  obtains 

x,  ^  od  ■+  (l-<Q>) 


155) 


which  is  the  form  for  overrelaxation. 
Alternately,  let 


then 


/- 


% 


A/  -* 


/- 


A' 


so  that  one  may  write  equation  51  in  the  form 


x; 


(56) 


(57) 


158) 


whioh  is •  of  oourse,  the  form  of  the  adapted  >«egstein  tech¬ 
nique  without  the  tildes  (**). 

As  in  suoceeslve  overrelaxation,  JLL,ia  described  as  the 
largest  proper  number  (eigenvalue)  of  the  matrix  where 

%.  «  (!-6Vc  (5.) 


9.  z  x  +  F  let) 
and  B  and  C  are  the  triangular  constituent  matrices  of  the 
matrix  a  of  the  system 

x,=  Ax  +  f  _l61) 

Though  the  system  under  consideration  is  of  the  form  AlJ=£ 
it  osn  be  sisply  demonstrated  that  this  form  can  be  rewrit¬ 
ten  in  the  form  of  equation  60*  (see  Appendix  B) 
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A® turn  to  the  assumed  correspondence  between 

.  Ci  0*0  ,  , 

ll(*i  ^ , 2, >£  +  AJb)  and  U  .  As  pre¬ 

viously  described,  tne  tilde  l<\>)  indicates  a  value  computed 
by  the  ••egsteln  technique  as  differentiated  from  a  value 
computed  by  the  method  of  successive  displacements.  In  a 
sense,  the  value  of  cverrelaxation  should 

also  have  a  tilde  since  it  represents  a  value  computed 
using  the  accelerating  factor  Ca) .  This  analogy  is  not  as 
straight  forward  as  it  seems,  however.  The  We gstein  tech¬ 
nique  as  adapted  by  Kaplan  and  Clark  is  not  applied  at 
every  iteration,  whereas  the  overrelaxation  factor  is. 
oince  for  those  iterations  waere  the  >Vegstein  tecnnique  is 
not  applied,  the  O^'is  set  equal  to  the  Upvalue ,  if  the 
application  of  the  •Yegstein  technique  is  restricted  to 
every  other  iteration,  the  tilde  loses  its  significance, 
further,  if  application  of  the  technique  is  restricted  to 
every  fourtn  iteration,  tne  tilde  loses  its  significance 
in  the  equation  for  computing  the  slope  "a",  oinee  the 
frequency  of  application  is  usually  equal  to  the  number  of 
nodes  to  the  deepest  node  from  the  boundary,  and  a  problem 
only  four  nodes  deep  is  quite  small,  the  frequency  of  ap¬ 
plication  is  generally  greater  than  every  fourth  iteration. 

One  more  important  point  must  be  made.  The  Weggteln 
technique  for  a  single  equation  oan  foroe  oonvergenoe  on  a 
divergent  Iteration  prooess.  tfhen,  however,  the  method  was 
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adapted  to  solution  of  a  system  of  linear  equations,  no  de¬ 
termination  was  made  as  to  whether  convergence  could  be 
forced  on  a  diverging  set  of  linear  equations.  On  the  other 
hand,  the  metnoa  of  successive  overrelaxation  cannot  foroe 
convergence  on  a  diverging  set  of  linear  equations. 
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Hi.  Irocedurea 


Computer  rrogrananlng 

The  basic  computer  used  for  this  study  was  the  13k  1620 
with  associated  IBM  1620  additional  memory  unit  and  13k  1622 
card  reader/punch.  The  1620  proved  to  be  a  relatively  slow 
system  and  some  of  the  computer  runs  would  have  taxen  an  ex¬ 
cessive  amount  of  time  on  this  system.  To  circumvent  this 
situation,  the  IBM  7090  computer  of  the  aeronautical  systems 
Division  was  employed  for  the  longer  runs.  As  it  turned  out, 
the  IBM  7090  using  the  relatively  new  FOKTRAN  IV  programming 
system  proved  to  be  well  over  100  times  faster  than  tne  13m 
1620  for  the  programs  run.  Lack  of  funds  prevented  exclusive 
use  of  the  IBM  7090,  however. 

The  JFOHTRAU  language  system  and  its  variations  were  used 
for  programming;  FQ^TkaA  II  for  tne  three  dimensional  pro¬ 
grams  run  on  the  IBM  1620  and  FORTRAN  IV  for  the  IBM  7090. 
lertinent  programs  are  included  in  appendix  0.  It  should 
be  emphasized  that  these  programs  were  written  in  the  slm- 
pllest  possible  manner  consistent  with  the  requirement  to 
provide  an  accurate  specific  set  of  data.  No  attempt  was 
made  to  generalize  a  given  program  to  fit  a  number  of  sit¬ 
uations.  The  individual  programs  were  not  specifically 
reviewed  for  possible  revisions  to  reduoe  running  time  or 
required  memory  spaces.  The  single  exception  Is  the  program 
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for  the  analytical  solution  of  the  sample  problem.  Running 
time  was  particularly  critical  on  tnis  program,  so  it  was 
thoroughly  reviewed  to  insure  that  the  program  steps  pro¬ 
vided  the  fastest  solution  procedure  consistent  with  FORTRAN 
language . 

standard  Used  for  Comparison  o£  speed  of  Convergence 

From  the  discussion  above,  one  can  visualize  the  diffi¬ 
culty  in  using  time  as  a  measure  of  speed  of  convergence. 

The  problem  is  further  complicated  by  the  faot  that  different 
cooputere  may  do  a  given  operation  in  different  ways  and  in 
different  times.  Thus  one  must,  so  to  speak,  nondimension- 
allze  the  rate  of  convergence  and  maohine  time  to  solution 
to  have  a  meaningful  standard.  This  is  generally  aooompllshed 
by  using  the  number  of  iterations  to  solution  as  the  standard 
between  Iterative  procedures  rather  than  time.  An  iteration 
is  defined  as  the  computation  required  to  make  a  single  im¬ 
provement  in  the  values  of  all  unknowns.  This  study  uses  the 
number  of  Iterations  to  oompare  the  various  methods  in  rate 
of  convergence. 

PtffnrtBlafi  l&i  Sfflff*fH9tt  Factor  t£  Jin  ISL  successive 
gmalnailai 

dlnoe  there  is  no  general  method  for  computing  the  optimum 
relaxation  faotor  ( )  prior  to  the  start  of  the  Iterative 
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process,  there  is  a  problem  of  what  factor  to  use  for  over* 
relaxation  in  the  comparison  of  metnods.  Ihe  emphasis  in 
this  study  is  on  the  practical  application  of  a  given  method 
to  the  solution  of  the  transient  heat  transfer  problem.  Thus 
to  be  practical  the  method  of  obtaining  the  relaxation  factor 
((0)  should  be  automated  to  the  maximum  extent  possible  and 
should  require  a  minimum  amount  of  prior  knowledge  of  the 
successive  overrelaxation  process  by  the  using  engineer.  Of 
all  the  methods  suggested  for  estimating  CD  ,  the  estimation 
of  the  maximum  eigenvalue  by  running  withal  fits  the  above 
criteria  better  than  other  methods  discovered  in  the  liters* 
ture.  This  method  was  used  to  provide  the  number  of  itera* 
tions  to  solution  for  the  method  of  successive  over relaxation 

The  above  selection  does  not  preclude  an  interest  in 
what  successive  overrelaxation  can  do  under  optimum  condi¬ 
tions,  however.  To  determine  the  optimum  factor  needed  for 
optimum  overrelaxation,  a  plot  of  iterations  versus  various 
relaxation  factors  was  used  (see  page  lb). 

Arror  Measurement 

•linos  there  usually  is  no  exact  procedure  for  determin¬ 
ing  the  error  associated  with  a  given  finite  difference  ap¬ 
proximation,  a  common  procedure  is  to  select  a  problem  with 
an  analytical  solution  to  use  as  the  standard  for  comparison 
with  the  finite  difference  solution.  This  procedure  was 
adopted  in  this  study. 
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The  Sample  problem  Amnloyed 

Again  using  the  oriterion  of  practicality,  a  three- 
dimensional  problem  was  chosen.  An  analytical  solution 
for  a  three-dimensional  problem  for  heat  transfer  over  a 
rectangle  with  a  oonstant  internal  heat  source  was  found 
in  literature,  specifically,  if  heat  is  produced  for  t>0 
at  the  oonstant  rate  q  per  unit  time  unit  volume  in  a  rec¬ 
tangle  of  dimensions  x*b,  ys  c,  z*d  and  the  surface  is 
kept  at  zero  degrees,  then  the  temperature  over  the  region 

0*x<  b,  C *  y<  o,  0*  z<d  is  given  by  («ef  1:363) 

«o  oo  00 

“•  1 1 1 

xso  rm so  /n=o 

Sew  {(2&+i)^e]Sin[(2irn4ijT?j$Cn{(2m*i)^J(l-€  ) 

(Z& +i)  (Zrm+tKZ'n+OoC 

choosing  to  work  in  three  dimensions  did  limit  the  size 
of  the  problem  to  be  solved  due  to  the  limited  storage  oapa- 
oity  of  the  IBM  1620  system,  When  the  above  problem  was 
solved,  a  cube  one  oentimeter  on  a  side  was  used  for  the  re¬ 
gion  of  interest.  It  was  then  neoessary  to  use  symmetry  to 
permit  seleotion  of  grid  sizes  as  small  as  0.1  om. 
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Diffusivity  (a)  ana  conductivity  ( K )  were  chosen  as 

4%  _ 

1.0  cm  /sec  and  1.0  cal/cm-sec-  0  respectively.  These 
values  approximate  the  parameters  for  copper  whose  true 
values  are  1.14  cm*/ esc  Tor  diffusivity  and  0.93  cal/cm- 
sec-^  for  conductivity  (Ref  1:497).  The  source  strength 
(q)  was  limited  to  2000  cal/cm3-sec  to  permit  final  temper¬ 
atures  within  the  region  to  be  of  a  physically  reasonable 
size.  Unless  otherwise  stated,  the  above  values  of  a,  K, 
and  q  were  used  in  all  the  solutions  obtained. 

Comparing  Spiral  Versus  Serial  aoan 

If  the  size  of  the  residuals  are  observed  after  each 
Iteration,  one  immediately  notices  that  the  residuals  for 
points  close  to  the  fixed  boundary  conditions  meet  the  con¬ 
vergence  criterion  much  sooner  than  those  furtherest  from 
the  boundary.  Intuitively,  one  could  see  where  the  repeated 
use  of  at  least  one  "correct"  value  for  the  adjacent  nodes 
to  compute  the  value  of  the  point  being  considered  could 
oonoelvably  hasten  the  convergence  of  this  point  to  the 
desired  solution.  Thus,  it  seems  logical  to  assume  that  any 
scanning  procedure  which  would  use  all  the  boundary  values 
as  soon  as  possible  in  the  solution  sequence  would  converge 
faster  than  a  serial  scan.  The  optimum  scan  for  achieving 
this  condition  is  to  spiral  Inward  from  ths  boundary. 


46 


GA/?hys/63-6 

The  particular  spiral  scans  used  are  shown  in  figures 
6  and  7.  As  one  can  see  from  the  figures,  the  tnree-dimen- 
slonal  case  where  symmetry  is  used  can  be  classed  as  a  con¬ 
sistent  scan  (see  page  22),  while  the  full  spiral  used  in 
the  two-dimensional  case  is  nonoonsistent.  The  fiul  spiral 
case  was  restricted  to  two  dimensions  by  the  lack  of  storage 
capacity  of  the  IBM  1620  system  used  in  the  computations. 

The  particular  full  spiral  used  was  chosen  to  permit  maximum 
use  of  the  "DO  loops"  of  the  FORTRAN  language  system,  thereby 
achieving  a  degree  of  simplicity  in  the  program. 

While  testing  various  types  of  spirals  for  consistency 
prior  to  the  eholoe  of  the  one  to  use,  a  conclusion  was 
reaohed  that  all  full  spiral  scans  are  nonoonsistent.  This 
conclusion  is  based  on  the  following  fact,  as  one  proceeds 
normally  Inward  from  the  boundary  following  the  direction  of 
the  higher  numbered  points,  the  elementary  squares  of  the 
mesh  show  consistency  in  the  region  close  to  the  boundary. 
Wnen  the  oenter  of  the  grid  is  approached,  nonconsistencies 
develop  due  to  the  ohanging  direction  of  the  higher  numbers, 
on  the  other  hand,  consistent  orders  exhibit  the  property  of 
a  single  general  direction  for  the  progression  of  higher 
numbered  grid  points. 

2&1  Msiftqg  21  aaani fl£  SM  Iterative  Process 
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1.  By  sweep  for  Z  *  0.1  Plane  (Sean 
through  sequence  1  and  2  then  move 
back  through  Z  In  1-2  seqv  'nee*  etc) 


II  — ■ *■  ns-*  2  /  / 

f  f  f  I  !  f 

/O  74— *»  /*4-*>  /44-M9f-*-/f4 

f  r  m  f  f 

f  !  \  f  f  t 

« f  r  m  r  r 

7  -*»  47  48-*- 49-*-  70-*-  7/ 

J  f  f  f  f  *  ] 

- 

2.  By  individual  grid  roint  for  Z  *  0.1 
Plane 

Pig.  6 


Diagram  of  the  Three-Dimensional  Semi-Spiral  Scan  Used 
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The  Iterative  process  is  generally  concluded  by  measur¬ 
ing  the  size  of  the  residuals  and  terminating  the  program 
when  the  size  of  the  residuals  meet  some  predetermined  cri¬ 
terion.  Some  possible  criteria  that  may  be  used  are: 

1.  Hequiring  the  average  of  all  the  residuals 
to  be  less  tnan  soma  specified  number. 

2.  requiring  the  sum  of  the  residuals  to  be  less 
than  some  specified  number  (see  page  22). 

3.  ^squiring  the  Individual  residuals  to  be  less 
than  a  specified  number--if  one  falls  the  test,  all  fall* 

The  fact  that  all  residuals  do  not  converge  at  the  same 
rate  has  been  mentioned  on  page  44.  In  methods  1  and  2  above, 
any  major  lnbalance  in  convergence  rate  oould  result  in  a 
relatively  large  inaoouraoy  in  the  points  which  oonverge  at 
the  slowest  rate.  In  practical  applications,  method  3  is 
used  to  avoid  this  potential  inaccuraoy.  further,  method  3 
is  the  most  sensitive  measure  of  convergence . 

In  view  of  the  above,  the  convergence  criterion  used 
in  this  study  required  all  points  to  be  less  than  a  given 
number  or  a  new  iteration  was  begun. 

selection  o£  the  Initial  &stlaate  o£  Temperature 

A  good  initial  estimate  of  temperatures  throughout  the 
region  of  Interest  oould  reduoe  the  number  of  iterations 
required  to  solve  the  problem  by  the  method  of  suooesslve 
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displacements  and  the  two  accelerating  techniques  used  in 
this  study.  All  runs  were  Bade  with  an  initial  estimate 
of  zero  degrees  for  the  temperature  of  all  internal  grid 
points,  however.  Aside  from  reasons  of  standardization 
among  solution  procedures,  starting  at  zero  removes  any 
screening  of  the  oonwergenoe  of  the  running  computation  of 
relaxation  faotor  (using  (h>sl)  by  differences  in  accuracy 
among  the  estimates  of  temperatures  at  individual  points. 

gm&am  tsL  4mi ui  ss.  Isis±  Saamaasi  gala e  ths. 

*aggtt>a  Saaftalflaw 

A  quick  oheek  was  set  up  to  investigate  whether  the 
adapted  ffegstein  technique  oould  foroe  convergence.  The 
test  was  to  use  the  adapted  ftegsteln  technique  on  the 
following  set  of  equations  which  were  known  to  give  a  diver¬ 
gent  solution  when  the  method  of  suooessive  displacements 
was  used  in  an  atteapt  to  solve  the  system. 

X,  +2Xt-  2X3  =  I 

X|  +  X*t  Xs 

2Xi  +2Xt  +  X3  zS  I63) 

This  set  is,  incidentally,  a  case  where  convergence 
oan  be  obtained  using  the  method  of  simultaneous  displace¬ 
ments  (see  footnote  page  1£) . 
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IV.  hesults 

aerial  versus  spiral  Scan 

The  results  of  the  comparison  of  the  rate  of  conver¬ 
gence  of  a  spiral  and  serial  scan  are  displayed  in  tables 
2  and  3. 

There  was  no  observed  advantage  in  the  use  of  the 
consistent  semi-spiral  scan  using  symmetry  as  applied  to 
the  three  dimensional  problem,  since  the  number  of  itera¬ 
tions  remained  the  same  for  both  serial  and  spiral  scan 
even  when  the  convergence  criterion  was  tightened  until 
the  computer  out-off  rounding  procedure  caused  all  residuals 
to  become  zero. 

Ahen  a  two-dimensional  problem  was  solved  using  a 
nonconslstent  full  spiral,  some  small  savings  in  iterations 
did  occur  for  the  spiral  scan.  Furthermore,  as  the  grid 
size  (h)  was  decreased  to  cause  the  solution  to  converge 
more  slowly,  the  savings  in  iterations  snowed  a  slight 
increase. 

Determining  the  Optimum  Relaxation  Factor 

Plots  of  number  iterations  versus  relaxation  factor 
are  shown  in  figures  9  and  10.  Figure  8,  demonstrating  the 
progress  of  the  temperatures  at  the  center  of  the  cub#  with 
increasing  time,  is  also  pertinent  to  this  particular  dls- 
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Table  2 


Iterations  to  Solution  for  a  Serial  Scan  and 
a  Three-Dimensional  Semi-Spiral  Seen 


Scan 

Type 

Total  Number  of  Iterations 

K< 

riiinl  Si  r.m  I 

0.005 

0.00005 

Jk  (Hia 

Serial11 

91 

145 

195 

Spiral4 

91 

145 

195 

Residual  Size  ] 

1.0 

0.05 

0.005 

Serial* 

31 

51 

67 

mmm u 

31 

51 

67 

®  t  ■  2.0  seo; 
b  t  ■  0.2  860  J 
Solutions  by 


h,  -  hy  .  he  -  O.i  cm 

hx  ■  *  hz  ■  0,125  cm 

the  method  of  successive  displacemente 
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Table  3 


Iterations  to  Solution  for  a  Serial  Scan  and 
a  Two-Dimensional  Pull  Spiral  Scan 


Scan 

Total  Number  of  Iterations 

Type 

Residual  Siee 

moamm 

0.005 

0.00005 

WEBSESM 1 

29 

66 

102 

Spiral0, 

26 

65 

101 

mmrvLMmm 

69 

218 

E35UHI 

67 

215 

R!t!H 

80 

688 

WEBESEM 

78 

684 

a  t  *  0.2  aeoj  hx  »  hy  ■  0.1  on 

b  t  «  0.2  eeoj  hx  ■  hy  »  0,05  cm 

c  t  ■  0.2  eeof  hx  •  hy  ■  0,025  cm 

Solutions  by  the  aethod  of  euooeeelve  dieplaoemente 
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ousslon.  This  figure  shows  that  the  temperature  has  essentially 
reached  a  steady  state  value  after  only  0.2  seconds  for  the 
parameters  a,  K,  and  q  used  in  this  study. 

Figure  9  dearly  demonstrates  a  dependence  of  the  relax¬ 
ation  factor  on  the  time  variable  by  the  shift  in  the  optimum 
factor  for  the  times  of  0.01  and  0.5  seconds.  The  point  that 
the  optimum  faotor  is  the  same  for  0.5  and  2.0  seconds  is  at¬ 
tributed  to  the  fact  that  for  these  times  the  temperatures 
within  the  oube  are  at  their  steady  state  values.  This  is 
further  confirmed  by  figure  10,  where  there  is  only  a  small 
shift  in  the  optimum  between  0.2  (nearing  steady  state)  and 
2*0  seconds  (well  into  steady  state  conditions). 

A.  spaoe  dependence  of  the  relaxation  factor  is  demon¬ 
strated  by  a  combination  of  the  two  figures.  For  a  time  of 
2.0  seoonds  there  is  a  shift  in  the  optimum  factor  between 
the  oase  of  unequal  grid  sizes  for  the  three  coordinate 
directions. 

4gaiHiaU»g  g&aigii  al  &s.  Adapted  wens  tern 
overtaxation 

Figures  11  and  12  display  the  results  of  making  a  run¬ 
ning  estimate  of  the  relaxation  faotor  while  Iterating  with 
the  method  of  suooesslve  displacements.  The  relaxation  fao¬ 
tor  dp  was  computed  for  eaoh  iteration  using  the  ourrent 
estimate  of  the  maximum  eigenvalue  This  eigenvalue  was, 
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lTO 1.19  1.33  1.46  1.59  1.72  1.85 


Relaxation  Factor 
Fig.  9 

▼aviation  of  lumber  of  Iterations 
for  Varying  Magnitude  of  Relaxation  Faotor 
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In  turn,  computed  in  three  different  ways:  as  the  norm  of  the 
residuals  (equation  £1),  as  the  first  power  norm  of  the  resid¬ 
uals  (equation  £3),  and  us  the  average  of  the  Aegsteln  slopes 
computed  from  £  0.1  /number  of  interior  nodes.  Figure  11  is 
for  equal  grid  sizes,  while  figure  12  is  for  the  more  general 
case  of  unequal  grid  sizes. 

The  important  features  displayed  by  figures  11  and  1£  are: 

1.  The  use  of  tne  norm  of  the  residuals  to  find  the 
optimum  factor  converged  rapidily  and  smoothly  to  a  value  that 
was  less  than  the  known  optimum. 

2.  The  procedure  using  the  first  power  norm  con¬ 

verged  to  a  value  closer  to  the  optimum  factor  (once  to  a 
▼slue  greater  than  COM  and  once  to  a  value  less  than  <a> ogt  ) 
than  the  method  using  the  norm,  but  at  a  slower  rate.  In  the 
case  of  equal  grid  sizes,  for  example,  if  one  were  to  aooept 
a  value  of  <o>  as  the  relaxation  factor  when  <*>^n','<  .001 

then,  using  the  norm  would  provide  a  value  in  12  iterations 
while  the  first  power  norm  does  not  meet  this  criterion  until 
after  16  iterations, 

3.  The  average  value  of  the  slopes  does  provide  an 
and  therefore,  an  CO  which  converges  to  the  same  value  as 

the  method  using  the  norm.  The  convergence  is  slower,  however. 
To  meet  the  same  criteria  as  used  in  item  2  above,  21  iterations 
would  be  required. 
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At  tne  seme  time  the  above  data  was  computed ,  the 
individual  slopes  or  the  interior  grid  points  were  observed 
every  fiftn  iteration.  The  primary  observations  of  interest 
are : 

1.  All  the  individual  slopes  converged  to  the 
same  value  and  this  value  was  identical  to  the  maximum  eigen¬ 
value  obtained  by  use  of  the  norm  of  tne  residuals. 

2.  The  individual  values  of  the  slopes  converged 

in  aifferent  ways  and  at  different  rates.  Characteristically, 
the  points  nearest  the  boundary  which  normally  were  the 
fastest  to  arrive  at  a  solution  li.e.  meet  the  convergence 
criteria  for  the  residuals)  started  at  values  of  slope  con¬ 
siderably  below  tne  final  value  and  converged  rather  slowly, 
un  the  otner  naad,  the  points  furtherest  from  the  boundary 
which  are  normally  the  last  to  reach  a  solution,  started 
witn  values  of  slope  greater  than  the  final  value  and  con¬ 
verged  faster  than  tne  points  close  to  the  boundary.  This 
feet  is  displayed  in  figure  13  for  two  extreme  grid  points. 

observations  on  the  Characteristics  of  the  Convergence  of 
the  solution  for  the  ^hree  procedures  Used 

figure  14  is  a  plot  of  the  progress  of  the  solution  with 
increasing  iterations  for  the  node  in  the  center  of  the  cube. 
The  higher  rate  of  convergence  for  the  two  accelerating  pro¬ 
cedures  is  to  be  expected.  The  significant  point  displayed 
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by  this  figure  is  the  contrast  between  the  modes  of  convergence 
of  toe  two  accelerating  procedures,  as  in  the  case  of  successive 
displacements,  the  solution  by  successive  overrelaxation  con¬ 
verges  through  a  smooth  progression  of  values  of  temperatures, 
in  contrast  to  this  smootn  progression,  the*  convergence  of  the 
adapted  «egstein  technique  is  ragged  with  spurts  of  accelera¬ 
tion  followed  cy  a  smoothing  process.  This  raggefiness  is  due, 
of  course,  to  the  intermittent  application  of  the  Aegstein 
technique.  Note  that  as  the  solution  approaches  the  final 
value  eacn  application  of  the  Aegstein  accelerating  factor 
tends  to  cause  the  solution  to  overshoot  the  final  value. 

This  is  tnen  followed  by  a  period  of  use  of  the  method  of 
successive  displacements  wnicn  pulls  tne  value  back  toward 
the  final  temperature,  as  a  matter  of  fact,  fer  every  solution 
obtained  by  use  or  the  adapted  «egstein  teccnique  the  process  ut 
the  convergence  criteria  and  the  program  terminated  on  an  itera¬ 
tion  that  was  computed  by  the  metnod  of  successive  displacements 
and  was  about  two-thirds  or  more  of  the  way  along  toward  the 
next  application  of  the  teennique.  tor  instance,  in  the  solution 
displayed  in  figure  14  tne  process  stopped  on  iteration  number 
56  where  the  last  prior  application  of  the  technique  was  on 
iteration  number  54  with  the  next  application  due  on  Iteration  60. 

There  is  some  danger  of  misunderstanding  being  created  by 
snowing  only  one  grid  point  in  figure  14.  une  must  recognize 
that  the  over-acceleration  does  not  have  to  occur  on  the  individ- 

64. 


i 


Ct  A/Phye/63-n 


65 


G^/rhys /6c-8 


ual  grid  point  to  cause  an  overshoot,  The  computational  pro¬ 
cedure  uses  .adjacent  grid  points  to  calculate  the  temperature 
at  the  point  of  interest,  so  any  over-acceleration  of  a  nearby 
point  will  produce  a  significant  effect  on  adjacent  nodes  and 
a  diminishing  affect  as  one  moves  away  from  the  node  that  was 
over-accelerated.  Thus,  tnere  is  an  interrelated  cause  and 
effect  relatlonsnip  tnat  is  not  apparent  from  figure  14. 

One  adaed  point  must  be  cited.  »ith  the  parameters  se¬ 
lected  for  solution  of  the  sample  problem,  the  rate  of  con¬ 
vergence  to  solution  turned  out  to  be  higher  than  desired  to 
demonstrate  the  effects  of  acceleration,  To  circumvent  tnis 
situation,  two  basic  procedures  were  used  to  increase  the 
number  of  Iterations,  oince  the  available  memory  capacity  of 
the  IBM  1620  prevented  decrease  of  grid  size  below  0.1  cm, 
tighter  convergence  criteria  than  one  might  expect  to  choose 
in  a  practical  problem  was  used,  second,  the  span  between 
initial  estimate  and  final  temperature  was  maximized  by  using 
times  that  brought  the  temperature  near  or  well  into  steady 
state  conditions.  Increasing  the  source  strength  and  there¬ 
fore,  the  final  temperature  had  a  negligible  effect  on  in¬ 
creasing  iterations  due  to  the  Increased  contribution  of  the 
q  term  of  equation  10  to  the  size  of  the  temperature  step  for 
eaoh  iteration. 

effectiveness  of  the  Accelerating  Techniques 


Table  4  shows  the  results  of  solving  the  sample  problem 
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Table  4 


Results  of  Accelerating  the  Method  of  Successive 
Displacements  by  Using  Successive  Overrelaxation  or 
the  Adapted  Wegstei n  Technique 


Item 

Method^^s^ 

Iterations 

Convergence 

Criteria 

Savings 

Run  la 

Run  2b 

Run  1 

Run  2 

Run  1 

Run  2 

Successive 

Displacements 

80 

91 

145 

.005 

.005 

.00005 

— 

— 

OverrelaxatiOn 
with  known 

<4)*PYimu*tr  c 

29 

32 

48 

.005 

.005 

.00005 

64* 

65* 

67* 

- y 

0 verr elaxat i  on 
by  Estimating 
(0  While  Run¬ 
ning  with  4 Uml 

42 

47 

69 

.005 

.005 

.00005 

48* 

48* 

52* 

Normal  * 

Adapted 
Wegstein 
Technique 

58 

58 

89 

.005 

_ 

.005 

.00005 

28* 

36* 

59* 

aFor  these  runs,  t  «  0,2  sec;  At  ■  0,2  sec;  all  grid  sizes 
were  equal  at  0,1  om 


^For  these  runs,  t  »  2,0  sec;  At  *  2.0  sec;  all  grid  sizes 
were  equal  at  0,1  cm 

cThe  optimum  values  of  the  relaxation  factor  for  these  runs 
were  taken  from  figure  . 

^ When  <*)('*♦')_  the  running  computation  of  <*>  was 

terminated  and  the  current  value  of  fc^was  then  used  to 
oontinue  the  solution  by  regular  oven’s! axat ion,  (i)  was 
computed  using -the  norm  of  the  residuals  to  find  ffj0 , 

*  The  deepest  node  was  six  from  the  boundary,  so  the  tech¬ 
nique  was  first  applied  on  the  12ih  iteration  and  on 
every  sixth  iteration  thereafter,  q  was  restricted  to 
a  maximum  of  -100. 
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for  two  different  times  using  the  same  grid  size  for  both 
problems.  The  primary  items  of  interest  in  this  table  are: 

1.  successive  overrelaxation  using  tne  optimum  relax¬ 
ation  factor  showed  the  fastest  rate  of  convergence  witn 
savings  in  iterations  of  the  order  of  65/&  over  the  method 
of  successive  displacements. 

2.  significant  savings  did  result  from  the  more  practi¬ 
cal  procedure  of  making  a  running  estimate  of  relaxation 
factor  from  the  residuals  of  the  method  of  successive  dis¬ 
placements. 

3.  The  adapted  •'•'egstein  technique  was  only  about  57/0 
as  effective  in  accelerating  successive  displacements  as 
overrelaxation  with  the  running  computation  of  relaxation 
factor. 

4.  x.jrcept  for  the  cese  of  the  adapted  <»egstein  tech¬ 
nique,  savings  remained  essentially  constant  between  runs 
using  the  same  convergence  criterion,  in  the  case  of  the 
adapted  </egstein  a  significant  increase  in  savings  {8/a) 
occurred  for  the  run  with  the  slower  rate  of  convergence. 

This  is  due  to  the  ragged  nature  of  the  convergence  of  the 
adapted  tfegstein  technique. 

5.  All  methods  showed  a  slight  inorease  in  savings 
for  a  decrease  (tightening)  in  convergence  criterion* 

since  the  two  methods  of  computing  a  running  estimate 
of  the  optimum  relaxation  factor  (first  power  norm  and  norm 
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of  the  residuals  of  the  method  of  successive  displacements) 
had  approximately  equal  rates  of  convergence,  a  comparison 
of  the  two  methods  was  made.  Table  5  displays  tnese  results 
ana  shows  that  the  first  power  norm  further  increased  tne 
savings  in  iterations  over  the  procedure  using'  the  norm.  It 
is  obvious  that  though  the  first  power  norm  was  slightly 
slower  converging  to  a  final  value  of  relaxation  factor,  the 
fact  that  this  method  produced  a  better  estimate  of  the 
optimum  factor  than  the  norm  overcame  the  slower  start. 

■forcing  convergence  with  the  adapted  .Vega teln  Technique 

figure  15  displays  the  results  of  attempting  to  force 
convergence  on  the  set  of  tnree  simultaneous  equations  given 
in  equation  65.  only  one  of  the  three  unknowns  is  plotted 
since  the  iterative  solutions  exhibited  the  same  general  diver¬ 
gent  characteristics. 

Initially,  the  normal  procedure  for  using  the  adapted 
Wegsteln  technique  was  attempted  with  the  initial  applica¬ 
tion  on  the  fourth  iteration  and  every  other  iteration  there¬ 
after.  ±his  procedure  completely  failed  because  the  diver¬ 
gence  ms  such  that  the  slope  test  failed  in  every  attempt 
to  apply  the  technique.  The  prooess  was  terminated  after 
10  iterations. 

A  second  attempt  was  made  using  the  same  frequency  of 
application  as  in  the  first  attempt  with  the  slope  test  deleted. 
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Table  £ 


A  Comparison  of  Methods  of  Making  a 
Running  Computation  of  Relaxation  Factor 
for  Use  with  Successive  Overrelaxation 


Method  of  Making  the 
Running  Computation  of  &) 

Iterations 

*  Savings  (%) 

RN 

tm 

Run  1 

Run  2 

First  Power  Normc 

40 

43 

50 

53 

Norm  of  the  Residuals- 

42 

47 

48 

48 

a  For  these  rune,  t  ■  at  ■  0.2  seo;  all  grid  sizes  equal 
at  O.l  cm 

k  For  these  runs,  t  ■  at  ■  2.0  seo;  §11  grid  sizes  equal 
at  0.1  om 

c  See  equation  23 

*  See  equation  21 

*  All  savings  compared  to  the  method  of  successive  dis¬ 
placements 

All  oonvergenoe  criteria  were  0.005 
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Recruits  of  Attempting  to  Force  Convergence  with 
the  Adapted  Wegstein  Technique 


GA/rhys/ 63-6 


The  results  obtained  are  the  ones  displayed  in  figure  14. 

In  a  third  attempt,  the  technique  was  applied  at  the 
second  iteration  and  on  every  other  iteration  thereafter. 

The  slope  test  was  again  deleted,  ^hese  results  displayed 
the  same  oscillating  divergence  as  the  second  case. 

finally ,  an  attempt  was  tried  using  the  same  frequency 
jf  application  as  the  third  case,  but  this  time  the  slope 
test  was  used  again.  The  slope  test  still  prevented  the 
application  of  the  accelerating  technique  and  the  process 
was  terminated  after  10  iterations. 

Though  the  solution  continued  to  diverge,  there  was  a 
cnange  in  the  mode  of  divergence  when  the  adapted  Wegsteln 
technique  was  applied.  When  an  attempt  was  made  to  solve 
the  system  using  only  successive  displacements,  the  solutions 
diverged  monotonically.  With  the  use  of  the  adapted  «egstein 
technique,  these  solutions  became  oscillating  and  divergent. 

arror 


Table  6  compares  the  error  at  selected  grid  points  re¬ 
sulting  from  the  three  solution  techniques  used  in  this 
study.  The  important  points  are  demonstrated  by  this  table: 

l.  The  solutions  for  tne  successive  overrelaxation 
and  adapted  rtegstein  metnoda  are  more  accurate  than  the  method 
of  suooeaelve  displacements.  This  is  a  result  of  the  technique 
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used  for  terminating  the  iterations.  Using  the  residuals 
as  tae  measure  of  convergence  results  in  more  accurate 
answers  for  a  rapidiiy  converging  solution  than  for  a  slowly 
converging  one.  Figure  16  demonstrates  this  caange  of  con¬ 
vergence  error  due  to  a  enange  in  convergence  rate.  In  the 
results  shown  in  table  6,  this  effect  is  somewhat  moderated 
by  the  rather  tignt  convergence  criteria  of  0.005  used  to  ob¬ 
tain  the  solutions. 

2.  The  solutions  for  2.0  seconds  show  about  one- 
fiftn  of  the  error  encountered  for  0.2  seconds,  oince  tne 
time  step  for  2.0  seconds  is  ten  times  greater  than  0.2 
seconds  and  truncation  error  *  is  a  function  of  grid  and  time 
step  size,  one  would  expect  the  reverse  results.  As  pre¬ 
viously  cited,  at  2.0  seconds  the  solution  is  well  into 
steady  state  conditions,  wnile  0.2  seconds  is  just  upproaoh- 
ing  steady  state.  These  results  demonstrate  the  fact  that 
wnen  one  desires  to  use  a  transient  neat  transfer  program  to 
obtain  a  steady  state  solution,  the  best  results  are  obtained 
with  a  time  step  tnat  puts  one  well  into  steady  state  con¬ 
ditions.  Taking  the  larger  time  step  permits  more  complete 
convergence  of  the  iterative  solution. 

Testing  a  Revision  to  the  Adapted  Aeastein  Technique 

The  apparent  relationsnip  between  overrelaxation  and 
the  adapted  wegstein  technique  suggested  oertaln  potential 
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revisions  to  tne  adaptea  negstein  tecnnique.  a  limited  a- 
mount  of  time  was  available  to  look  into  possible  revisions, 
so  the  investigation  was  restricted  to  a  relatively  obvious 
revision  that  one  might  make  to  smooth  the  ragged  convergence 
of  tne  regular  technique, 

Tne  revision  attempted  was  to  replace  4)by  its  equi¬ 
valent  in  terms  of  ^I'lsee  equation  20)  in  equation  45  to 
obtain  a  new  expression  for  computing  nqw  of 

2 

/  +s/l-a.  (40a) 

Tnis  equation,  in  conjunction  with  the  slope  test  restrict¬ 
ing  values  of  "a"  to  0<a<i,  would  keep  the  acceptable 
values  of  "q"  within  the  limits  specified  by  equation  46,  It 
was  hoped  that  this  revision  would  also  take  advantage  of 
potential  offered  by  the  individuality  of  the  modes  of  con¬ 
vergence  of  the  slopes  for  individual  grid  points.  As  pre¬ 
viously  suggested,  the  nodes  with  the  slowest  rate  of  con¬ 
vergence  to  solution  have  slopes  that  converge  through  a  se¬ 
quence  of  values  that,  wnen  used  In  equation  40a,  would  give 
better  q’s  as  the  values  are  greater  than  optimum,  oinoe  the 
convergence  was  assumed  to  be  smoother  than  the  regular  adapted 
A'egstein  tecnnique,  the  acceleration  was  applied  on  every  iter¬ 
ation  beginning  with  Iteration  six  for  the  first  run  with  tnis 
revision,  a  second  run  was  then  made  with  the  technique  ap¬ 
plied  every  iteration  beginning  with  iteration  two--  the  earll- 
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lest  possible  Iteration  that  values  of  "a"  can  be  computed. 

The  results  of  this  revision  are  shown  in  table  7.  The 
complete  failure  of  this  revision  was  apparently  due  to  a 
combination  of  the  small  size  of  the  q's  used  and  the  number 
of  rejected  slopes  that  resulted  in  unacoelerated  nodes. 

Computer  Running  Times 

During  the  initial  attempts  to  obtain  an  analytical 
solution  using  the  IBM  1620  computer,  a  wholly  unexpected 
difficulty  arose.  The  time  required  to  obtain  an  analytical 
solution  using  the  IBM  1620  was  impossibly  long.  -For  in¬ 
stance,  it  took  8  minutes  and  20  seconds  per  grid  point  to 
do  the  triple  summation  using  equal  maximum  summation  indices 
of  only  5.  a  graphical  estimate  of  the  time  to  do  a  summa¬ 
tion  with  maximum  indices  at  20  showed  an  approximate  time  of 
two  nours  per  grid  point,  further,  another  plot  of  solutions 
versus  maximum  summation  Indices  indicated  thut  the  summa¬ 
tion  indices  should  be  at  least  20  to  obtain  the  three-place 
aocuracy  desired. 

This  situation  forced  the  use  of  the  faster  IBM  7090 
system.  The  same  problem  solved  on  the  7090  took  only  13. & 
seoonds  per  grid  point  for  maximum  indices  of  20,  and  36 
seconds  per  point  for  Indices  of  30. 

This  incident  prompted  interest  in  other  comparisons, 
for  insttnoe,  it  took  the  7090  only  49  seconds  to  solve  a 
two-dimensional  problem  of  400  nodes  that  took  218  iterations 
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Table  7 


Tabulation  of  the  Effects  of  a  Revision  to 
the  Adapted  Wegstein  Procedure 


Method 

Iterations 

Humber 

Savings®’  (#) 

Regular  adapted  Wegstein 

58 

36 

Revised:**  Apply  technique  on 
every  iteration  beginning 
with  iteration  6 

72 

21 

Revised » •  Apply  technique  on 
every  iteration  beginning 
with  Iteration  2 

70 

- J 

?" 

Savings  over  the  method  of  euoeesnive  dieplac  einent  a 
b  Using  equation  40a  in  lieu  of  equation  40  to  compute  q 


78 


GA/?hys/63-8 


whereas,  the  IBM  1620  used  roughly  one  hour  and  45  minutes 
to  solve  a  three-dimensional  problem  of  125  nodes  that  took 
only  91  iterations. 
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V.  Conclusions  and  Recommendations 


Conclusions 

It  is  apparent  that  the  adapted  ^egstein  technique 
and  successive  overrelaxation  have  the  same  theoretical 
basis,  Thefrixe  differences  are  in  the  method  of  appli¬ 
cation  or  the  accelerating  factor  and  the  used/permitted 
size  of  the  accelerating  factor. 

The  results  obtained  in  this  study  lead  one  to  con¬ 
sider  the  adapted  -vegstein  technique  as  a  gross  overre¬ 
laxation  procedure.  This  technique  apparently  obtains 
its  acceleration  power  from  a  series  of  gross  overrelex- 
ations  of  a  few  nodes  and  the  ability  of  the  method  of 
successive  displacements  to  distribute  this  large  accel¬ 
eration  to  other  adjacent  nodes  while  pulling  the  value 
of  tne  overrelaxed  node  bacit  in  line,  Ahe  problem  with 
this  procedure  is  the  lacx  of  control  of  the  spurts  of 
over  acceleration  in  the  late  stages  of  the  iterative 
process,  how  one  might  control  this  to  better  advantage 
is  an  unanswered  question.  The  results  obtained  when  an 
attempt  was  mads  to  smooth  the  acceleration  tentatively 
indicate  that  no  simple  revision  of  the  adapted  ^egstein 
technique  is  going  to  provide  a  significant  increase  In 
the  ability  of  the  procedure  to  reduce  iterations.  There 
is  no  intent  to  suggest  that  the  method  will  never  approaoh 
the  method  of  overrelaxation  in  total  savings  of  iterations* 


80 


GA/Thys/63-8 

The  results  obtained  cannot  support  such  a  conclusion.  The 
best  chanoe  for  improvement  seems  to  be  in  the  use  of  the 
individuality  of  the  convergence  of  the  slope  "a"  of  the 
nodes.  To  use  this  characteristic  effectively,  however,  one 
must  find  a  way  to  isolate  the  disturbing  effect  on  the 
slopes  of  having  the  adjacent  nodes  accelerated.  In  this 
way  the  slope  test  will  not  reduce  tne  number  of  nodes  that 
receive  acceleration. 

A  possible  procedure  for  isolating  the  slopes  would 
be  to  generate  a  matrix  of  q's  on  the  nth  iteration  without 
applying,  the  technique.  ‘Then  the  q's  generated  on  the  nth 
Iteration  would  be  applied  on  the  (n+l)st  iteration.  This 
would  be  followed  by  a  period  of  Iterations  by  successive 
displacements  to  smooth  the  slopes,  Ahen  the  same  procedure 
of  delayed  application  oould  be  used  again,  and  so  on.  Using 
the  normal  method  of  computing  q  (equation  40)  would  permit 
values  that  tend  to  overaccelerate,  so  that  this  aspect  of 
the  adapted  Aegstein  technique  would  be  retained.  A  disad¬ 
vantage  of  this  method  would  be  the  requirement  to  store 
another  matrix  equal  in  size  to  the  number  of  unknowns. 

The  results  using  the  optimum  factor  are  obviously  the 
maximum  savings  that  can  be  aohieved  by  the  process  of 
estimating  ID  for  successive  overrelaxation.  It  is  also 
likely  that  this  optimum  overrelaxation  represents  the  max¬ 
imum  attainable  savings  for  the  adapted  wegstein  technique. 
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Yet,  a  chance  series  of  events  could  cause  e  gross  over- 
acceleration  early  in  tne  process.  ±hen  after  a  few  iter¬ 
ations  of  successive  displacements,  tne  temperatures  could 
oe  orcugnt  bac&  to  a  value  near  tneir  final  value  oefore  the 
srrootner  successive  overrelaxation  process,  one  cannot,  of 
course,  build  a  procedure  on  the  basis  of  chance,  but  the 
results  ootained  permit  the  speculation  that  it  might  be 
possible  to  increase  tne  probability  of  occurrence  of  this 
cnance. 

fhe  adapted  <«eejstein  technique  suffers  from  two  other 
significant  disadvantages  wnen  compared  with  successive 
overrelaxation.  Ahe  technique  requires  a  more  complex  pro¬ 
gram  and  the  storage  of  one  additional  matrix  of  values 
equal  in  size  to  the  number  of  unjcnownB.  "ith  computers  the 
size  of  the  709G  with  a  capacity  of  40, COO  words  available, 
the  storage  space  required  may  not  be  a  critical  requirement, 
inis  does,  nowever,  restrict  the  use  of  the  method  to  rel¬ 
atively  large  computer  facilities.  In  view  ol  these  mech¬ 
anical  disadvantages  of  programming  and  storage,  the  adapted 
-'agate in  tecnnique  would  nave  to  be  able  to  produce  greater 
savings  tr.an  the  procedure  of  maxing  a  running  estimate  of 
the  relaxation  factor  b-fo  -r  the  technique  is  widely  accepted. 

*'or  successive  overrelaxation,  the  results  indicate 
that  the  practical  procedure  of  making  a  running  estimate 
of  relaxation  factor  does  result  in  significant  savings  of 
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iterations.  Phase  results  also  indicate  that  tue  best  pro- 
cedui'e  to  use  is  the  first  power  norm  when  estimating  the 
factor,  i'nis  is  fortunate  since  the  first  power  norm  re¬ 
quires  less  computation  tnan  the  norm  because  the  residuals 
are  not  squared  prior  to  summing,  to ore  conclusive  proof  is 
required,  nowever.  The  simplicity  of  the  sample  probiam 
could  nave  otimized  the  convergence  oi  the  first  power  norm. 
Despite  tne  results,  tue  nore  general  case  could  be  that  tue 
usa  of  tne  first  power  norm  degrades  the  procedure  compared 
to  the  norm  since  the  two  were  relatively  close  in  tne  re¬ 
sults  obtained,  supporting  this  assertion  are  the  facts 
that  in  the  two  cases  of  comparison  tabulated,  the  norm  did 
reach  a  final  value  sooner  and  was  consistently  lower  tnan 
optimum  wnile  the  first  power  norm  produced  a  factor  that 
for  one  run  was  greater  than  optimum  ana  for  the  other  less. 

Despite  the  simplicity  of  the  test  to  see  if  the  adapted 
•Vegstein  technique  can  force  convergence,  one  must  conclude 
tnat  it  cannot.  Tnis  follows  logically  from  the  fact  that 
tne  adapted  ••egstaln  technique  nas  the  same  theoretical  basis 
as  successive  overrelaxation  and  successive  overrelaxation 
has  been  shown  to  be  theoretically  incapable  of  forcing  con¬ 
vergence.  furthermore,  the  technique  uses  and  aepends  heav¬ 
ily  upon  the  method  of  suoceseive  displacements  which  is  a 
significant  difference  from  the  ^egstein  technique  applied  to 
a  single  equation.  Thus,  the  adapted  «egstein  technique  cannot 
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force  convergence  If  tne  method  of  successive  displacements 
diverges. 

oince  there  was  no  real  break  away  of  one  method  in 
the  results  obtained  for  the  comparison  of  serial  ana  spiral 
scans,  the  results  seem  inconclusive,  on  the  other  hana, 
the  results  can  support  s  conclusion  tiiat  no  cnange  in  scan 
procedure  will  produce  a  significant  acceleration  of  the 
solution  by  successive  displacements. 

respite  their  limitations,  rower’s  results  suggest  the 
second  possible  conclusion.  Her  results  show  only  a  small 
increment  of  change  in  the  number  of  iterations  by  successive 
displacements  for  the  various  scan  procedures  used,  mis  is 
consistent  witn  the  results  obtained  in  this  investigation. 

The  absence  of  a  source  of  acceleration  from  a  relaxation 
factor  or  eauivaieat  to  provide  the  large  temperature  incre¬ 
ments  per  iteration  in  the  initial  iterations  supports  the 
second  conclusion,  i'urtnermore ,  since  the  full  spiral  scan 
is  nonconsistent,  it  is  unlikely  that  using  successive  over¬ 
relaxation  or  the  adapted  ./egstein  technique  with  the  full 
spiral  f o  —  further  acceleration  will  prove  to  be  as  rood  as 
acceleration  with  a  serial  scan.  Thus,  the  facts  available 
from  this  study  seem  to  indi  cate  that  no  significant  advantage 
can  be  acnieved  oy  use  of  other  than  a  consistent  serial  scan. 
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Recommendations 

This  investigation  provides  a  theoretical  basis  1'or  the 
adapted.  -vegstein  technique  which  was  not  previously  defined, 
but  fails  to  use  this  knowledge  in  any  broad  experimentation 
designed  to  improve  the  technique.  There  is,  tnerefore,  a 
justification  for  further  study  and  comparison  to  improve  the 
technique  and  draw  a  firm  conclusion  as  to  the  better  metnod 
of  accelerating  the  method  of  successive  displacements. 

It  is  recommended  that  any  further  expansion  of  this 
work  be  done  with  sample  neat  transfer  problems  of  the  type 
solved  in  a  practical  situation  in  order  tnat  performance  may 
be  better  related  to  real  situations. 
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Appendix  A 


Constructing  Finite  EjJTgESSga  Equation 


Two  particular  types  of  partial  differential  equations 
are  to  be  considered  in  this  appendix,  the  elliptic  equation 
and  the  parabolio  equation. 

A  typical  elliptio  partial  differential  equation  is 
the  Laplaoian, 


V*U.  =  O  (64) 

U  S  f(x, n,l) 

The  boundary  conditions  specify  u  on  the  boundary  of  the 
region  considered. 

The  parabolic  partial  differential  equation  considered 
has  as  its  form 


with  the  following  boundary  conditions 


(65) 


on  the  boundary 

The  approximations  to  the  partial  derivatives  are 
formed  from  truncated  Taylor  series  expansions  in  the  fol¬ 
lowing  manner.  Expanding  about  point  (Xt.yt.i^to)  one  may 
obtain 
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u.6c0+ax,«$.,a,;0  * 


2.  3X 


(66) 


or 

u(x0- 


=  ax]*,*.*'** 


u. 


21  3*‘K*,a,A 


•  •• 


(67) 


Truncating  equations  66  nnd  67  at  the  second  partial  der¬ 
ivative  term  and  adding  the  twc  expressions,  one  obtains  the 
approximation  for  the  second  partial  derivative  term  of 


2tlL  ~  UiWM.'i. A.  t.) -2UU..4..U&J+U 

ax1  ax>*  (68) 

The  expression  for  the  first  partial  derivative  is  obtained 
by  subtracting  the  two  equations  to  obtain 

sv  U6<«»  ^•tg«.£o+A£)-TJ(X«/|e.3b,  £•)....  .  . 

at  At  (69) 


Note  that  as  the  intervals  Ax,  Ay,  Az  or  At  become  smaller, 
fewer  and  fewer  terras  are  required  to  obtain  a  good  approx¬ 
imation  since  the  interval  is  raised  to  successively  higher 
powers  with  each  new  tern  making  such  terms  as 


less  significant.  Thus,  if  one  reduces  the  time  and/or 
space  interval,  a  more  accurate  solution  should  result. 
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This  fact  does  constitute  a  simple  test  of  the  validity  of 
one's  solution*  For  the  test*  one  halves  the  grid  size  used 
for  the  solution  at  hand  and  determines  a  new  solution  with 
this  smaller  grid  size.  If  the  two  solutions  give  answers 
that  are  in  reasonable  agreement  with  each  other*  then  one 
has  some  assurance  that  the  original  solution  is  correct* 
Starting  first  with  the  elliptic  equation  (number  64), 
when  one  substitutes  the  expressions  given  in  equation  68 
for  the  partial  derivatives,  the  finite  difference  approxi¬ 
mation 


[M  (rf  GrtfJ  ^  rnx)* 


(70) 


4> 


(At? 


is  obtained.  One  may  generalize  this  equation  using  the 
following  notation 

AX  s  s  +  JL 4X  s  A*  JL *  3*3*  M* 

w  »•♦**•*  f (To 

4  1  sit,  *4,*  *•  +  ■*4*®—  S....OV 

Then  equation  70  may  be  written  in  the  general  form 


(72) 
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If  the  fora  simplifies  to 

UU.**)  =  f  [u(ui,3,-£)+ua-i,3M)+ua,3>Ml) 

+  t7(x,3'i,it)+U(x,3,il+i)+Ufr,3,^-i)J  (73) 

For  two  dimensions,  the  form  becomes 


=  ^[uCl+z^J+Wi-i^hIT^+O+U^i^)]  (74) 

Thus,  U(i, j,k)  is  expressed  in  terms  of  values  of  U 
at  nodes  one  point  away  in  the  direction  of  the  coordinate 
axes*  Since  the  given  boundary  values  only  specify  U's  on 
the  boundary,  each  TJ(i,J,k)  is  expressed  in  terms  of  other 
unknown  values  at  the  interior  node  points  or  a  combination 
of  unknown  interior  points  and  known  boundary  points*  The 
result  is  a  system  of  simultaneous  linear  equations  equal  in 
number  to  the  number  of  interior  nodal  points* 

As  an  example  of  how  the  system  of  equations  is  con¬ 
structed,  consider  the  two  dimensional  grid  with  mesh  points 
numbered  as  shown  and  boundary  values  lettered* 


44 


m 


r  »  »  a 

-1-T-T-* 

e  -4-1-4-  9 

n  i  n 

aTeT i 


With  boundary  values  prescribed.  Us  thru  Up  are  known  and 
may  be  considered  oonstant.  The  unknowns  oan  be  described 
by  the  following  nine  finite  dlfferenoe  equations  (for  the 
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nine  Interior  nodes)* 


4U« 

+  Ue  +U4. 

+  \Jp 

+  Ub 

=  O 

-4U»4  U8 

+  U 5 

+UC 

«  0 

Ua-4tfs 

+U4  4US+U/ 

=0 

u, 

-4V4 

+Ui 

+Ui  4  Uo 

SO 

Ua  +U4 

-40s  +U6 

■n Ja 

=  0 

Oe 

4  Us 

-404 

4U9 

=  0 

Ua  +Um.  +U* 

-4Ui  +  U§ 

=  0 

U* 

+U 5 

4Ui  “40®  4U9 

so 

Uj 

+U4 

4U8  -4V9 

=  0 

Rearranging  and  placing  the  constant  boundary  values  on  the 
right,  the  system  oen  be  written  in  the  matrix  form 


-4  I  O 
I  -4  / 

O  I  >4 
I  O  O 
O  (  0 

O  O  I 
OOO 
O  0  O 
OOO 


I  o 
0  I 
o  o 
•4  I 
I  -4 
O  i 
I  o 
O  I 

o  o 


o  o 
o  o 
I  o 
o  / 

I  o 

-4  O 
O  -4 
0  i 
I  o 


0  0 

m  m 

u, 

c, 

0  0 

Ua 

Cl 

0  0 

u. 

Cl 

0  0 

U4 

C4 

1  0 

Us 

s 

C, 

0  1 

U4 

Ci 

f  0 

U» 

Cl 

-4  1 

Ut 

Ct 

1  -4 

A 

LC4J 

(76) 


or  in  matrix  notation 

AU  =  C  (77) 

It  la  interesting  to  note  that  the  matrix  A  la  equal  to 

ita  transpose  (or  S  «  H  T  )  wfaioh  indioatee  that  A  is  a 
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As  in  the  case  of  the  elliptic  equation,  equation  68  is 
used  in  the  parabolic  equation  to  approximate  the  second  par¬ 
tial  derivative.  Then  equation  69  is  used  to  approximate  the 
partial  with  respect  to  time.  This  time,  however,  there  is 
a  choice  of  time  level  to  use  in  the  approximation  cf  the 
second  partial  derivative.  Which  level  of  time,  t  or  t  ♦  At, 
is  set  into  the  approximation  of  v2u  is  optional,  but  the 
two  levels  lead  to  two  different  forms  of  the  parabolio 
finite  difference  representation. 

To  obtain  the  explioit  form,  one  ohoses  the  level  t  to 
get 

17  ■tfl  I  U(l+h3t'&/4) 

’ 

+  ['  -  aa*A(* 

This  equation  is  known  as  the  explioit  form  sinos  U(x,y,s,t+  t) 
is  explicitly  expressed  in  terms  of  known  values  at  the  level 
t. 

Since  u  at  any  x,  y,  s  and  some  starting  time  t^  are 
given,  the  solution  oen  be  started  at  t  ■  t^,  where  the  U’e 
on  the  right  side  of  the  above  equation  are  knowns,  and  solve 
for  the  unknown  U(i,  j,k,t*At)  for  each  mesh  point  in  the  region. 
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Cnee  U  is  known  for  all  mesh  points  at  the  level  t^  ♦  At 
the  process  is  repeated  to  obtains  U’s  at  the  level  of 
t^  +  2 At,  and  so  on  to  the  time  desired. 

Though  this  explicit  form  does  not  lead  to  the  require¬ 
ment  to  solve  a  system  of  equations,  it  does  hove  the  dis¬ 
advantage  of  being  restricted  by  stability  considerations. 
Unless  the  following  criteria  are  met,  an  instability  exists 
that  causes  a  growth  of  errors  from  time  step  to  time  step 
which  ultimately  reaches  the  point  where  error  dominates  the 
solution  (Ref  4:92  ).  For  this  particular  finite  differ¬ 
ence  representation,  stability  will  exist  if 

ft  *  “g“  (80) 

for  the  three-dimensional  problem,  or 

^  <  T  (81) 

for  the  two-dimensional  problem j  where  for 

ft  =  (82) 

One  can  see  that  any  change  in  space  (or  time)  increment 
must  be  accompanied  by  v.  proportional  change  in  time  (or 
space)  increment. 

Alternately,  vau  oen  be  set  at  the  level  t  +  / »t  to  get 
the  implicit  form  (i.e.  "a  group  of  components  of  U  at  level 
t  +  mt  are  defined  simultaneously  in  such  an  interrelated 
manner  that  it  is  necessary  to  solve  a  linear  subsystem  for 
the  whole  subset  of  components  at  once  before  a  single  one 
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can  to  determined"  (Ref  :267) 

[' +  !aAtlii  +Tj+Zillja‘*‘'6MAt) 

TJ(  1+1, +  U  (/.,$- 1 

TIM  L  w  J 

For  the  case  of  equal  grid  size  in  the  three  coordinate 
directions,  one  may  write 


- 

Lx 

-W 


(83) 


D  +  */Jua,*,*i,i+4i) 

-  it  [Uteti,  3,  i,  t+Ai)+UU-i, £,i,i+4£) 

+  +U(jL,1+lt'£)Jb4At) 

_  (84) 

+  <^3.  it  <,  it  A*)  +ua,fr,i-/,i*4;fc)j =Ua,g,ii) 

where 


For  two  dimeneicne  and  equal  grid  aiaes,  one  may  write 


[l+4A,]UU,  j-.i+ai) 

•  H  [U^t  J,^,itA£)+UU-/, j.itai) 

(85 

+  Ua,}+‘,Jt+Aii  +  Uf/,3-i,^yAt)  =  UU,g,it) 

3inoe  U  at  any  x,  y,  a,  and  some  starting  tine  t^  la 

specified  by  given  initial  oonditlons,  the  aolution  can  be 
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commenced  at  t^  At  and  U(i,;j,k,t  )  would  be  known  and 
the  remaining  U's  form  a  system  of  simultaneous  equations 
at  the  level  tA  +  at.  This  need  to  solve  a  system  of  equa¬ 
tions  resulted  from  the  finite  difference  representation 
of  the  elliptic  equation,  end  a  comparison  of  the  matrix 
forms  of  the  elliptic  end  parabolic  equations  is  of  interest. 
Taking  the  same  sample  grid  system  that  was  used  for 
the  elliptic  equation,  the  system  of  equations  for  the  im¬ 
plicit  form  are 


O  +  A/DUi  -AUe.  “A  U4 

-  A  Ui  +Ci*+4A)17t  -  AUj  -  A  Ur 

-  AUa  +(/+4A)tri-  tlVi 

-  A  Ui  +  (I* -  JLJJt  -  *Uv 

-AUa  -AU4  +0*4A)U#-AU*- 
-/lUj  -AUf  4  -A.V9 

-AU4  +0+4A)Uv- A  Ug 
-4Uf  -  A Ur  4*  ('/+4A)Ue-AU4 
-AU4  -  AUI  +  04 4^t)U«J 


r-  (Vp+Vk)  4  U/(£ 
s-Uc4U*,jft 
s-flW+U r)+U** 
a  -U. 

AU«=U*;fc 

-  <86> 

=-  U4+U».i 
=  -  Ctt,+U4)+Uf,jt: 


where  UB|<t  ■  U(i, j,k,t)--lmown  from  initial  conditions 
or  previous  computations 

Uwn  *  U(i#  j,k, t+at) 

In  matrix  form  the  Bystem  beoomee 


95 


GA/Phys/63-8 

In  matrix  form  the  system  becomes 


0+4A) 

- 1 

0 

•A 

0 

O  O  0  o 

•  « 

u, 

•  ^ 
C, 

-  n 

(l-M/t)  -ft. 

o 

-a 

o  O  0  o 

% 

Pi 

0 

-/l 

0*44)  O 

o 

-A.  0  0  ° 

u> 

<* 

'fl 

0 

0 

(l+4A) 

-A 

o  -a  o  0 

a 

o 

-fl 

0 

-A 

0+4n)  -A  o  -ji  o 

• 

Cj’ 

0 

0 

•H 

0 

“A 

(l-MA)  o  o-^ 

% 

C* 

0 

o 

o 

•A 

0 

0  (iM 4)-/l  o 

T7> 

Cr 

0 

0 

o 

O 

o  -a  (H4A)  -  A 

% 

C| 

_  o 

o 

o 

0 

0 

-a  o  -fi 

to 

AJ 

which  in  matrix  notation  is 

I'u  S  C  '  (68) 

These  matrloes  have  the  same  general  form  as  those  for  the 
elliptic  equation  and,  thus,  the  implicit  form  may  be  thought 
of  as  equivalent  to  solving  an  elliptic  equation  at  each  time 
step. 
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Appendix  B 


Sample  Matrix  Forms 


The  system  under  consideration  is  of  the  form  ATJ  *  0 
where  sample  expanded  forms  are  given  in  equations  76  and 
86  in  Appendix  A.  The  example  parabolic  equation  used  in 
Appendix  A  will  be  used  as  an  example  in  this  Appendix  to 
develop  the  forms  of  the  matrices  that  comprise  the  expres¬ 
sion  for  determining  the  pertinent  eigenvalues  of  the  sol¬ 
ution  procedures  considered  in  this  study. 

To  begin  this  development,  divide  each  equation  of 
the  system  given  in  equation  86  by  the  common  diagonal  ele¬ 
ment  of  the  A  matrix  to  obtain  (Ref  3:128) 
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After  moving  the  off-diagonal  terms  to  the  right,  the  system 
can  be  written  in  matrix  notation  as 

U  =  Su  +  (go) 

where 
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The  matrix  3  may  be  subdivided  into  two  triangular  constit¬ 


uent  matrices  (Ref  1 :134) 
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Now,  note  that  the  B  matrix  contains  the  elements  associated 
with  the  (n  +  l)st  iteration  and  the  C  matrix  with  the  nth 
iteration  of  the  method  of  successive  displacements  for  the 
consistent  serial  scan  3hovm  in  Appendix  A. 
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One  may  then  write 

0(M*"  =  Bu  (**l>  +  CV™  +c‘  (95) 

which  may  be  rearranged  to 

ggc^w)  -  £0**°  +  c'  (96) 


(l  -  B)Uim*,>s  CU(m,)  +  C'  (97) 

,  (i-gfc  V**+(f-6J'c'  (98) 


which  is  equivalent  to  the  form 

X  =  fl,X  +  F 


(61) 


Prom  the  section  on  theory  one  can  find  two  descriptions 
of  the  partitioning  of  the  A  matrix  to  find  the  eigenvalues 
of  the  method  of  successive  displacements.  Equation  16 
partitions  the  A  matrix  into 


(16) 


to  arrive  at  the  detenninantal  equation  for  the  eigenvalues 


of 

d%t  (t\l  +t)  90  (18) 


In  this  case  the  forms  for  the  submatrices  are 
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Alternately, 

equation 

59  gives  the  form  of  partitioning 

of  the  A  matrix  as 

a=/r-s)"5  <»> 


For  this  form,  Mrs  Fadeeva  advises  that  the  eigenvalues 
may  he  found  from  the  determinantal  equation  (Ref  :135) 

det[c  -  (i-§)fj]  =  o  (102) 

The  forms  of  the  B  and  C  matrix  are  given  in  equations  93 
and  94*  The  I  matrix  is  an  identity  matrix  of  the  usual 
form  of 
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Note  that  the  relationships  between  the  D,  E,  and  F  and  I, 
C,  and  B  matrices  may  be  expressed  as 
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The  equivalence  of  Equations  18  and  102  can  be  demon¬ 
strated  in  the  following  manner.  First,  rearrange  18  to 

dml  +  =0  (107) 

Now  if  one  multiplies  this  expression  by  the  reciprocal  of 
the  diagonal  element  of  the  D  matrix,  solutions  of  the  deter- 
minantal  equation  are  uneffected.  This  action  results  in 

det{Wf+NB+^fB=0  <loe) 

which  by  the  use  of  equations  104  to  106  can  be  rewritten 
in  the  exact  form  of  equation  102. 
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Appendix  C 


Definitions 


Diagonal  Dominance : 

For  a  matrix  to  possess  the  characteristic  of  diagonal 
dominance,  the  following  relationship  must  exist  for  the 
elements  of  the  matrix: 

N 

iQill  >  £ 

property  (A) :  Ihef  4: £43) 

a  square  matrix  A  of  order  A 
(A)  if  there  exists  a  permutation 
is  diagonally  block  tridiagonal. 

■spectral  hadius:  (hef  10:94) 

The  spectral  radius  is  defined  as  the  maximum  of  the 
moduli  of  the  eigenvalues  of  a  given  matrix. 

Dlriohlet  Difference  acuatlon; 

The  finite  difference  equation  that  represents  the 
partial  differential  equation  V*U.  s  o  In  tne  region  id  th 
H8f  on  the  boundary. 


jD»  M 

1109) 


is  said  to  have  property 
matrix  such  that 
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heal dual: 

The  difference  between  two  auooesaive  improved  values 
at  a  given  node,  l.e. 
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SOLUTION  OF  THREE  SIMULTANEOUS  EQUATIONS  BY  THE  GAUSS 
SIEOEL  METHOD  WITH  OPTION  FOR  WEGSTEIN  ACCELERATION 

SPECIFY  DIMENSIONED  VARIABLES 

DIMENSION  A (3, 3)  ,B(3),XNEW( 3 ) , X ( 3 ) , DI FF ( 3 ) 

DIMENSION  MX (  3 ) , WXNEW (3) , WD I FF 13) 

READ  INPUT  DATA 

REAU.AI i, I) ,A< 1 , 2  > , A I  1, 2) ,A<2, l ) ,A( 2,2) ,A(2, 3) 
RfcAD,A(3,l),A(3,2),A(3,3) 

RCAD.Bd)  ,B(2>  ,B(  3) 

READ.XI 1 ) , X(2) ,X( 3)  ,ERR 
READ, I APLY, APPLY iWMAX.WSET 

SET  INITIAL  VALUES  AND  ESTIMATES 
DO  60  1  =  1  ,  3 
OIFFI I )=0.0 
60  WDI FF ( I ) =0.0 

COUNT  IS  AN  EXIT  CONTROL  THAT  COUNTS  THE  NUMBER 
OF  PASSES  THRU  THE  COMPUTATIONAL  SEQUENCE 
CQUNT=0.0 

SET  ITERATION  COUNTER 
I f NUM*0 

33  I  TNUM=I  TNUMd 

COUNT  =  COUNT ♦ I . 0 
PRINT  104 
104  FORMAT  I/) 

PRINT  6.ITNUM 

6  F0RMAT(NHITERATI0N,I4/J 

WHEN  A  GIVEN  RESIDUAL  FAILS  TO  PASS  THE  CONVERGENCE 
TEST,  ICHEK  IS  SET  *  I  AND  NO  MORE  RESIDUALS  ARE  CHECKED 
UNTIL  THE  NEXT  ITERATION 
ICHEK=0 

COMPUTATIONAL  SEQUENCE 
N  =  0 

PRINT  30 

30  FORMAT ( I3H  VALUES  OF  X  • 3X , I4HX I  N*  l  )  -  XIN)  /) 

46  N*N+1 

IFIN-l 11,1,40 

40  IF(N-2)2,2,41 

41  IFIN-3) 3,3,42 

1  XNEwm  =  l./AIl,l)*(B(l)-A(l,2)*X(2)-AIl,3)*X(3)  ) 

DIFFm  =  XNEW(l)-X(l) 

GO  TU  43 

2  XNEW  (  2)*l./A(2,2)*(0(2)-A(2,l)».X(l)-A(2,3)»X(3)  ) 
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DIFF(2)=XNEWI2)-X(2) 

GO  TO  43 

3  XNEW(3)=l./A(3,3)*(B(3)-A(3»l)*X(l)-A(3»2)*X(2)) 

DIFFC3)*XNEW( 3)-XC3l 
C 

C  PRINT  OUT  RESULTS  EVERY  ITERATION 

43  PRINT  62  »  XNEW ( N ) «  0  I FF ( N ) 

62  FORMAT  (2E14.8/I 
X { N ) =XNE  WIN) 

C 

C  BYPASS  THE  ACCELERATING  TECHNIQUE  BY  SETTING  SWITCH  ONE  OFF 

I  F  (  SENSE-  SWITCH  1)71,46 
C 

C  SET  UP  WEGSTEIN  VALUES  FOR  THE  FIRST  ITERATION 

71  IFl I TNUM- 1 ) 44, 44 ,45 

44  WO  1 FF I N ) *D I FF I N ) 

WX(N)*XNEW(N) 

X ( N ) =XNEW ( N ) 

GO  TO  46 

INPUT  IAPPLY  DETERMINES  THE  ITERATION  FOR  THE  FIRST 
APPLICATION  OF  THE  TECHNIQUE 

45  I F ( I TNUM- IAPLY)47, 53,53 

THE  INPUT  APPLY  DETERMINES  HOW  OFTEN  THE  TECHNIQUE  IS  APPLIED 

63  IFIC0UNT-APPLY)47,48,48 
47  WDIFFIN  ) *XNE Wl N ) -WXI N) 

WXIN)*XNEW(N) 

X ( N ) =XNE W I N ) 

GO  TU  46 

4B  1 F ( N-3 )  63,54,54 
54  COUNTED. 

THE  ADAPTED  WEGSTEIN  TECHNIQUE 
63  WTEST*DIFF(N)/WDIFFIN) 

WFACT*WTEST/(WTEST-1.  ) 

PRINT  i07,WTEST,WFACT 

107  FORMAT  I/BHSLOPE  *  , E 14 .8 , 2X,9HF ACTOR  «  ,E14.B) 

BYPASS  THE  SLOPE  TEST  BY  SETTING  SWITCH  2  OFF 
I F ( SENSE  SWITCH  2)  105,50 
105  IF(WTEST)47,47,49 

49  IFIWTEST-l. 150,47,47 

TEST  ACCELERATED  SOLUTIONS  AGAINST  ESTABLISHED 
CONVERGENCE  CRITERIA  WHEN  TECHNIQUE  IS  USED 

50  IF(A8S(WFACT)-WMAX)52,52,5l 

51  WFACT«-WSET 
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52  WXNLW(N)=WFACT*WX(N)*(  l .-WF ACT ) »XNEW( N ) 

PRINT  6 1 1  WXNEW ( N  ) 

61  FORMAT ( /llHWXNEW(N)  *  ,E14.8/) 

WO  I  FF (N)=WXNEW(N)-WX!N) 

WHEN  THE  WEGSTEIN  TECHNIQUE  IS  BYPASSED 
TEST  GAUSS  SIEDEL  SOLUTIONS  FOR  CONVERGENCE 
IF!  AHS(Wi)IFF(N)  )-ERR)  BC,B0,8l 

80  I C  HEK  =  l 

81  WX(N)=WXN£W!N) 

X  ( N )  =WXNEW ( N) 

GU  TO  46 

C 

C  TEST  TO  SEE  IF  A  NEW  ITERATION  IS  REQUIRED 

42  IF(ICHEK)  72,72,35 
72  DO  33  J *  1 , 3 

IF(ABS(DIFF(J) J-EftR)  33,33,35 
33  CONTINUE 
C 

C  PRINT  OUT  INPUT  DATA  FOR  REFERENCE  PURPOSES 

DO  100  J-  1 , 3 
DO  100  1=1,3 

100  PRINT  101, All, J), I, J 

101  FORMAT! /9HA ( I , J )  *,E14.8,6H  I  *  ,13, 6H  J  =  ,13) 
DU  102  1=1,3 

102  PRINT  103,  B(I) 

103  F0RKATI/9H  0(1)  *  ,E14.8) 

C 

13  STOP 
END 
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Compute  lKn+1)  and 
0(n+l)-0(n)  by  Sue*  Slap. 


Xnereaae  Zt. 
Counter 

0(n+l)  =  0(n) 
0(n+l)  =  U(b) 

0(b+1)  -  0(b) 
=  U(n)-U(n-1) 


Sat  aad  Store 
0(n+l) »  0(n+l) 
0(b)  *  «{b+1) 
U(b4-1)-3(b) 

0(n+l)-0(n) 


'till  O'a  N 

sheen  canp- 
\  uted? / 


/  la  \ 

ythii  tha  N 

\  tha  teeh» 
\8i?ia/ 


Sat  aad  Store 
O(n+1)-O(n*l) 
0(a)  =  0(b  1) 

0(b+1)-0(b)s 

0(b+1)-0(b) 


BLOCK  DXAOMM 
Of  TBS 


qr  a/a  -  1 


Zheseaae 

Zt*  eouBter 
0(n*l)»0(n) 
0(tt*l)«8(n) 
0(b+1)-3’(b)» 
B(b)-6’(b-1) 


9S"***  M 

»(«t 1)  qO(a) 
(l-q)0(ae.l) 
aad  $a*l)4f(a) 
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C  heat  transfer  in  cube-ncnsteady,  constant  source  by  gauss  sieoel 

l  USING  WEGSTE I N  CONVERGENCE  TECHNIOUE 

C 

DIMENSION  Ul 7,7,7) ,UOLO 17,7,7) , WOI FF l 7 , 7 , 7 ) , WU { 7, 7 , 7 ) 

C 

C  READ  INPUT  DATA 

RF AO  61,  LASTX.LASTY.LASTZ 
61  FORMAT  (313) 

READ53,UlFUS,COND,DELT  ,SORC 

53  FORMAT ( 4F  1 0 . 0 ) 

READ  54,LRR,XGRI0,YGR!C,ZGRID,TMAX 

54  FORMAT ( 5F  10 . 0 ) 

READ  56,IRITF,JRITF,KRITE,JITER,KITER 
56  FORMAT  t  D  I  4  ) 

READ  460,  IAPPl.Y,  APPLY, WMAX 
460  FORMAT! 14, 2F10.0) 

C 

C  SfcT  BOUNDARY  VALUES 

00  3  K*  1 , LAST L 
DO  3  J=1,LASTY 
DU  3  1  =  1 , L AS TX 
601  UOl.D  (  I  ,  J  ,  K  )  =0. 0 
U( l , J,K)=0.0 
3  WU ( 1 , J , K ) =U ( I , J , K ) 

C 

C  COMPUTE  CONSTANTS 

A=(2./(XGRID**2)  )M2./( YGRID»»2> ) ♦ ( 2. / ( ZGR ID**2 ) > 

8  =  0 1 FUS  » DEL  T 
C  =  H/ ( XCR 10**2) 

D=B/ ( YGR I D**2 ) 

£  =  B/ ( ZGR 10**2) 

G= l . / ( l  .  ♦ B*  A ) 

F=G*3»S0RC/C0ND 

C  SET  TIME  S1EPS  ANO  ITERATION  COUNTER 

T I  ME  =0 . 

14  T I Mt  =  I  1  MENDEL  T 

I TNUM-0 
COUN  T  =  0 . 0 

CHOOSE  EITHER  PUNCHED  OR  PRINTED  OUTPUT  BY  SETTING  SWITCH  ONE 
OFF  FOR  PUNCHED  OUTPUT  OR  ON  FOR  PRINTED  OUTPUT. 

1 F ( SENSE  SWITCH  1)81,80 

80  PUNCH  50, TIME 

GO  TO  l 

81  PRINT  50, TIME 

50  FORMAT  I / 1 8HNEW  TIME  INCREMENT, FiO. 3) 

l  ITNUM=ITNUMU 

COUNT “COUNT ♦  l  •  0 
I F ( SCNSE  SWITCH  1)83,82 

82  PUNCH  51,1 TNUM 
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GO  TO  84 

83  PRINT  51*1 TNUM 

51  FORMAT  I/15X,16H ITERATION  NUMBER, 14/) 

C 

C  COMPUTATIONAL  STEPS 

84  ILAST*LASTX-1 
JLAST»LASTY-1 
KLAST*L AS  TZ- 1 

KTEST  =  KLASTMLAST*JLAST 

ICHECK-0 

00  2  K*2 , KL AST 

00  2  J*2 , JL AST 

DO  2  I*2,ILAST 

KK«I*J»K 

P*G*C*(U< I ♦ i , J,K )*U( I -I, J,K )) 

0«G»D*(U( l, J*l,K)*U(I,J-l,K> ) 

R*G*E*(UII, J,K*l )+UCI,J,K-l) ) 

S»G«UOLO( I i J,K) 

UNEW*F«-P*Q*R  +  S 
01 FF»UNEW-U< l , J, K) 

ON  THE  FIRST  ITERATION  THE  INITIAL  VALUESS  OF  THE  WEGSTEIN 
VARIABLES  ARE  SET  UP  BUT  NO  WEGSTEIN  COMPUTATION  IS  DONE 
I F ( I TNUM— 1)125, 125,470 
125  WOIFF(I,J,K)«OIFF 
WUI l,J,K)«UNEW 
GO  TO  102 

INPUT  I  APPLY  TELLS  WHAT  ITERATION  WILL  BE  THE  FIRST  TO  APPLY  THE 
WEGSTEIN  TECHNIQUE  ON 
470  I F ( I TNUM- I APPLY)  452,450,450 

THE  INPUT  APPLY  DETERMINES  HOW  OFTEN  THE  TECHNIQUE  IS  APPLIED 
450  I F ( CQUNT-APPL Y )  452,109,109 

452  WDIFFI I , J , K ) *UNE W-WU ( I , J , K ) 

WU( I , J,K)*UNEW 
GO  TO  102 

THE  WEGSTEIN  TECHNIQUE 
109  I  FI WDIFF ( I , J, K ) )  709,289,709 
709  WTEST»DIFF/WOIFF( I, J,K) 

IF(WTEST)  289,289,100 

100  IF! WTEST-l# I  101,289,285 

101  WEGF AC »WTEST/ I WTEST-l. ) 

REJECT  VALUES  OF  THE  WEGSTEIN  ACCELERATING  FACTOR  THAT  ARE  ABOVE 
A  PRESET  VALUE  ANO  SET  THE  FACTOR  EQUAL  THE  PRE  SET  VALUE 
IF  I A8SF I WEGFAC )-WMAX )  453,453,454 
454  WEGFAC—WMAX 
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453  WUNEW=WEGFAC»WUII,  J ,K > ♦ t l .-WEGFAC ) »UNEW 
WOIFFI I , J,K)*WUNEW-WU< I  ,  J»K  > 

WU( I , J»K)*WUNEW 
U( I ,  J,  K ) =WUNEW 
C 

C  TEST  SOLUTIONS  AGAINST  ESTABLISHED  CRITERIA 

I F  t ICHECK)103,103» 116 

103  IF (ABSFI WOIFFI I. J,K) )-ERR 1116, 116,105 

105  I  CHECK  «  1 

116  IF(JITfeR-J)  110.110.290 
110  1 F ( K I TER-K )  112,112,290 

112  I F ( SENSE  SWITCH  1)114,113 

113  PUNCH  115,1,J,K, WUNEW, WC IFF(1,J,K),WTEST 

GO  TO  290 

114  PKINT  1 15 , I , J , K, WUNEW, WO IFF (I,J,K),WTEST 

115  FORMAT (5X,3I4,3E14.8) 

GO  TO  290 

289  U( 1 , J,K)*UNEW 

WOIFFI I , J,K)*UNEW-WU( I , J,K) 

WU(  I , J , K ) *UNE  W 

290  IF(KTEST-KK)  481,481,480 

481  C0UNT»0.0 

480  1 F I WT6ST-1.0)  2,107,107 

102  U( I , J,K)»UNEW 

I F ( ICHECK ) 106, 106, 107 

106  1F(ABSF(0IFF)-ERR)107,107,U8 

118  1CHECK-1 

107  IF( JITER-J) 119,119,2 

119  I F I K 1 TER-K )  120,120,2 

120  1 F ( SENSE  SWITCH  1)  122,121 

121  PUNCH  52,1, J,K,UNEW,DIFF 
52  FORMAT  ( 5X, 3I4.2E14.0 ) 

GO  TO  2 

122  PKINT  52,1, J,K,UNEW,01FF 
2  CONTINUE 

C 

C  SET  VALUES  ONE  GRID  POINT  BEYOND  PLANES  OF  SYMMETRY 
L«LASTX 

00  4  N«2,LAST2 
00  4  M*2,LASTY 

4  U(L»M,N)*U(L-2,M,N) 

M»LASTY 

00  5  L*2 »LASTX 
00  5  N»2,LASTZ 

5  U(L,M,N)*U(L, M-2 ,N I 
N*LASTI 

00  6  L*2,LASTX 
00  6  M*2,LAST£ 
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6  U(L,M,N) =U(L,M,N-2  ) 

C 

C  PRINT  RESULTS  OF  ITERATION 

IFI ITNUM-l) 1, 1,150 
150  1 F ( 1CHECK ) 10, 10,  1 

10  1 F ( SENSE  SWITCH  1)88,87 

87  PUNCH  60 

GO  TO  89 

88  PRINT  60 

60  FORMAT  ( // 10X , 26HF INAL  RESULTS  OF  ITERATION//) 

89  00  7  K=2,KLAST ,KRI TE 
00  7  J=2, JL AST , JRI Tfc 
DO  7  I*2,lLAST,IRlTE 

I F ( SENSE  SWITCH  1)91,90 

90  PUNCH  130,1, J,K,U(I«J,K) 

00  TO  7 

91  PRINT  130,I,J,K,U!I,J,K) 

130  FORMAT  ( 10X , 3  1 4, E 14 . 8 ) 

7  CONTINUE 
C 

C  CHECK  TO  SEE  IF  LAST  TIME  INCREMENT  HAS  BEEN  REACHED 
I F ( T  MAX- T I  ME ) 1 7, 17,  12 
C 

C  PREPARE  FOR  A  NEW  TIME  INCREMENT 

12  00  15  1 =2,L ASTX 

00  15  J»2,LASTY 

00  lb  K*2,LAST7 

15  UOLO ( I,J,K)*U(l,J,K) 

GO  TO  14 
C 

C  PRINT  OUT  INPUT  DATA 

17  I F ( SENSE  SWITCH  l  )94,93 

93  PUNCH  55 ,01 FUS,CONO,OEL T , SORC 
PUNCH  57, XGR10, YURI D*ZGR ID, TMAX, ERR 
PUNCH  58, IR[TE,JRITE,KRITE,KITER, J1TER 
PUNCH  59, IAPPLY, APPLY, MAX 

GO  TO  13 

94  PRINT  55 ,01 FUS ,COND,OEL  T , SORC 

55  FORMAT! 1 1HD IFFUS I V l TY,F 10.4 , 6X, 12HC0NDUCT I VI TY, F10. 4/ 

310H0ELT  A  T I  ME ,F 10.4 ,6X, 15HS0URCC  STRENGTH, F10. 4 ) 

PRINT  57,XGRID,YGR!0,ZGRID, TMAX, ERR 

57  FORMAT ( /6HXGR I0»,F10.4,6X, 6HYGRI D* »F10.4»6X,6HZGRID*,F10.4/ 
48HMAX  TIME.F10.4, 18HST0P  WHEN  ABSDIFF-,F 10.5, 15HIS  ZERO  OR  LESS) 

PRINT  58, IRITE, JRITE,KRITE,K!TER, JITER 

58  FORMAT ( /25MUUTPUT  CONTROLS  ARE  IR I TE, 1 3, 2X ,5HJR 1 TE , 1 3/ 

55HKRITE, 13, 7H  KITER,I3,7H  JITER, 13) 

PRINT  59, IAPPLY, APPLY, MAX 

59  FORMAT  (  /7H I  APPLY*  ,  14 , 4X» 6HAPPLY*  » F 10 .5', 4X »4HMAX*  , F 10. 5 1 

13  STOP 
END 
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HLAT  TRANSFER  IN  CU8E-NCNSTEADY,  CONSTANT  SOURCE  BY  GAUSS  SIEDEt 
WITH  PROVISION  FUR  COMPUTING  RELAXATION  FACTOR  AND  WCGSTEIN  DATA 

SPECIFY  DIMENSIONED  VARIABLES 

DIMENSION  U(7,7,7) ,UOLO(7,7,7) ,DIFF{7,7,7) 

READ  61,  LASTX.LASTY.LASIZ 
I  FORMAT  (313) 

READ63,DlFJS,COND,DELT,SORC 
3  FORMAT ( 4F 10 • 0 ) 

READ  54,ERR,XGRID, YGR I  0 , ZGR I D, TMAX 
54  FORMAT ( 5F 10 • 0 ) 

READ  56,IRITE,JRITE,KRITF, J I  TER , K I  TER , CHOSE , BOUND .CARD, COMP 
56  FORMAT ( 5I4.4F10.0) 

READ  124, I TSET.CONVRG 
124  FORMAT  (I3.FI0.0) 

SET  BOUNDARY  VALUES 
DO  3  I -  l , L AS  TX 
DO  3  J=l,LASTY 
DO  3  K*  l ,  l ASTZ 

CHOOSE  METHOD  OF  SETTING  INITIAL  VALUES  BY  LETTING  BOUND*0.0 
TO  READ  FROM  CARO  INPUT  OR  BOUND*l.O  TO  SET  ALL  POINTS  EQUAL  ZERO 
I F ( BOUND )  100,100,101 

100  READ  102, U< I, J,K) 

102  FORMAT  (EU.B) 

UOLD(  I  ,  J , K ) =U ( 1 ,J,<) 

D I FF ( l , J , K  J  *0. 

GO  10  3 

101  UOLD( I , J,K) =0.0 
D 1 FF l f ,J,K)*0. 

U( 1 , J,K)=0.0 

3  CONTINUE 
GRELAX=0. 

AVtA=0.0 

WKELAX*0.0 

POINTS*(LASTX-2)«(LASTY-2)*(LASTZ-2) 

l L AS  T*L AS  f  X- 1 

JLAST«LASTY-l 

KLAST-LASTZ-l 

WHEN  APPLICABLE  CAN  TAKE  INITIAL  VALUES  FROM  PUNCHED  OUTPUT  BY 
USING  CARO  *0.0  0R-1.0  ANO  SETTING  BOUNO  *  1.0 
IF (CARO)  2 00,200,201 
200  DO  202  K*2 , KL AST 
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L>0  202  J*2,JLAST 
00  202  I =2 • I  LAST 

202  READ  203,  III  I,  J,K) 

203  FORMAT  (22X.E14.8) 

L  *1 ASTX 

00  210  N*2,LASTZ 
00  210  M*2,LASTY 

210  U(L,M,N)*U(L-2,M,N) 

M»LASTY 

DO  211  L*2.LASTX 
00  211  N»2,LASTZ 

211  U ( L • M»  N ) *U( L  « M-2 , N ) 

N*L  AS  TZ 

00  212  L *2 , L ASTX 
00  212  M*2,LASTZ 

212  U(L,M,N)*U(L *  M,N-2 ) 

GSETA«0.0 

C 

C  COMPUTE  CONSTANTS 

201  A*(2./(XGRI0*»2) ) ♦ ( 2. /  ( YGR I  0**2 ) ) ♦ ( 2. / 1 ZGR ID**2 ) ) 

B=DIFUS»OELT 
C*8/  l XGR 1 0*  »2 ) 

D»B/( YGR10**2) 

E*B/(ZGRID*»2) 

G* l • /  1 1 • ♦B*A ) 

F=G»B»$QRC/CONO 

C 

C  SET  TIME  STEPS  AND  ITERATION  COUNTER 

T IME-O* 

14  TIME«T I ME^DELT 

I TNUM»0 
IC0UNT*0 

CHOOSE  EITHER  PUNCHED  OR  PRINTED  OUTPUT  BY  SETTING  CHOSE*0.0 
FOR  PUNCHED  OUTPUT  OR  CH0SE-1.0  FOR  PRINTED  OUTPUT 
IF (CHOSE) 80,80,81 

80  PUNCH  50, TIME 

GO  TO  1 

81  PRINT  50, TIME 

50  FORMAT! /18HNEW  TIME  I NCREMENT , F10. 3 ) 

1  I TNUM«I TNUMtl 

IC0UNT*IC0UNT*1 
TOT  ALA*0.0 
T0TGS*0.0 
IF (CHOSE  182, 82,83 

82  PUNCH  51, ITNUM 

GO  TO  84 

83  PRINT  51,1 TNUM 

51  FORMAT  (/15X, 16H ITERATION  NUMBER, 14/) 

C 
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C  COMPUTATIONAL  STEPS 

84  00  2  K*2 • KL AS T 
00  2  Ja2 , JL AST 
00  2  1=2,1LAST 

P*G*C«(U( !♦!, J,K }*UII-1  ,J,K)  ) 

0=G*U» ( U ( I ,J+1,K)+U( I »  J- I »  K ) ) 

R=G»E*(U(  I ,  J,K*1  )  ♦O  (  I  .  J.K-D) 

S=G*UOLO ( I , J,K) 

UNEW  =  F  +  P*0«-R*S 
I F ( SENSE  SWITCH  1)  190,191 

190  RUNL W  =  GR£L AX*UNE  W* (l.-G8ELAX)*U( I , J,K  ) 

UNEW-RIJNEW 

191  DIFOLO=OIFF(I,J,K) 

01 FF ( I , J , K ) =UNEK-U ( I , J  ,  K ) 

IFISENSt  SWITCH  1}  116,98 
98  IF(Cf)MP)  240,240,241 

240  T0Tl»S=r0T(iS*A8SF(DIFF!  I  ,  J,  K  )  ) 

GO  TO  242 

241  T0TGS-TQTGS*DIFF(I»J»K)»*2 

242  IF(ITNUM-l)  116,116,196 

196  1 F ( I T  SE  T- I  COUNT )  20,20,116 

20  WTEST=OIFFI I,J,K)/D!EOLC 

TOFALA=TOTALA«-WTEST 

WEGFAC=WTEST/(WTfcST-l.) 

IFtSCNSF  SWITCH  2)  192,116 

192  I F (CHOSE )  750,750,751 

750  PUNCH  951,1 ,J,K,WTEST,WEGFAC 
GO  TO  116 

751  PRINT  951,1 ,J,K,WTEST,WEGFAC 
951  FORMAT  IiI3,2614.8) 

116  1FIJ1TER-JI  110,110,2 
110  IF(KITER-K)8,8,2 
8  IF(CH0SE)85,85,86 

85  PUNCH  52 , I , J , K,UNE W , 01 FF ( I , J , K ) 

GO  TO  2 

86  PRINT  52,I,J,K,UNEW,DIFF(I,J,K) 

2  U( I , J,K)*UNEW 

I F ( I TSE T-1C0UNT )  975,975,976 

975  I C0UNT*0 
C 

C  SET  VALUES  ONE  GRID  POINT  BEYOND  PLANES  OF  SYMMETRY 

976  L*LASTX 

DO  4  N«2,LASTZ 
00  4  M*2,LASTY 

4  U(L,M,N)«UIL-2,M,N1 

m=lasty 

DO  5  L*2  ,  LASTX 
00  5  N*2 , LAST/ 

5  U(L,M,N)*U(L,M-2,N) 

N*LASTZ 
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DO  6  L*2  ,LASTX 
i)0  6  M*2 » LAST l 
6  U(L,M,N)»U(L,M,N-2) 

GSDlF*$QRTF(TOTGS) 
i  F  (  SENSE:  SWITCH  li  198,99 
99  IFUTNUM-l)  111,111,123 
123  GSETA*GSOIF/GOENOM 
ROLD*GRCLAX 

GRELAX«2.0/ ( l . *SORTF U .-GSE TA ) ) 

AVE A* TOT ALA/ POINTS 
WRELAX*2.0/( l.*SQRTF( l.-AVEA) ) 

111  GDENOM«GSOIF 

IFICHOSE)  243,243,244 

243  PUNCH  18l,GSETA,GSOIF,GRELAX,WRELAX,AVEA 
GO  TO  24* 

244  PRINT  18l,GScTA,GS01F,GRELAX,WRELAX, AVE A 

181  FORMAT  (8HGSETA  «  ,E 14 . 8, 2X ,9HGSDI FF  *  ,E14.8,11H  GRELAX  ■  ,E14.8 
2,/,llH  WRELAX  «  ,E14«8,4X, 7HAVE A  *  ,E14.8) 

245  I F ( SENSE  SWITCH  3)  743,198 
743  IF(ITNUM-l)  198,198,197 

197  I F ( ABSF ( GRELAX-ROLO l-CONVRG )  199,199,198 
199  PRINT  180 

180  FORMAT  ( 50H0VERREL AXAT I  ON  FACTOR  HAS  MET  CONVERGENCE  CRITERIA/ 
232HSW1  ON  TO  OVERRELAX,  PRESS  START) 

PAUSE 

C 

C  TEST  SOLUTIONS  AGAINST  ESTABLISHED  CRITERIA 

198  DO  10  K~2,KLAST 
DO  10  J*2 , JLAST 
00  10  1*2,1  LAST 

IFCABSFCOIFFl I,J,KI)-ERR)  10,10,1 
10  CONTINUE 
C 

C  PRINT  RESULTS  OF  ITERATION 

IF(CHOSE)87,87,88 

87  PUNCH  60 

GO  TO  89 

88  PRINT  60 

60  FORMAT  ( // 10X , 26HF I NAL  RESULTS  OF  ITERATION//) 

89  00  7  K«2,KLAST,KRITE 
00  7  K*2,KL AST ,KRl TE 
00  7  J*2  » JLAST , JRI TE 
00  7  1*2, IL AST, IR I TE 
IF (CHOSE ) 90, 90, 9 l 
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90 

91 
52 
7 

C 

C 

c 

c 

12 


15 


C 

c 

17 

93 


94 

55 


57 


58 


766 

13 


PUNCH  52 , I , J , K , U ( I i J , K ) • D 1 F F  t I  ,  J  »  K  ) 

GO  TO  7 

PRINT  52 i l ,  J ,  K ,  U  ( 1 1 J »  K ) , D I F  F { I ,  J  ,  K  ) 

FORMAT (  10X,3I4,2E14.8) 

CONTINUE 

CHECK  TO  SEE  IF  LAST  TIME  INCREMENT  HAS  BEEN  REACHED 
IF  I TMAX-TIME) 17, 1 7, 12 

PREPARE  FOR  A  NEW  TIME  INCREMENT 
00  15  I=2,lASTX 
DO  15  J=2,LASTV 
UO  15  K=2,IASTZ 
UOLUII. J,K)«U(I, J.K) 

GO  TO  14 

PRINT  OUT  INPUT  DATA 

IF(CH0SE)93,93,94 

PUNCH  55,01FUS,C0ND,0ELT »SORC 

PUNCH  57,XGRID,YGRI0,ZGKID,TMAX,ERR 

PUNCH  58, IRI TE, JR  I TE ,KR I TE , K 1  TER , J I  TER, CHOSE 

PUNCH  766 , CONVRG 

GO  TO  1J 

PRINT  55,01FUStC0ND,DELT , SORC 

FORMAT! 11HOIFFUSIVITY,F10.4,6X, 12HC0NDUC  TIVITY,F10.4/ 
310HUEL T  A  TIME,F10.4,6X, 15HS0URCE  STRENGTH, FI  0.4  I 
PRINT  57,XGRID,YGRI0,ZGR10,TMAX,ERR 

FORMAT ( /6HXGR ID* ,F 10.4 ,6X,6HYGRI D*,F10.4,6X, 6HZGR ID*,F10.4/ 
48HMAX  TIME.F10.4, 18HST0P  WHEN  ABSDI FF-, E 14.8 , 16H  IS  ZERO  OR 
PRINT  58, IRITE,JRITE,KRITE,KITER,JITER, CHOSE 
FORMAT ( /25HOUTPUT  CONTROLS  ARE  I R I TE , 1 3 , 2X ,5H JR  I TE , I  3/ 

55HKR I TE , I 3,2X,5HKITER,I 3 , 2X , 5H J I  TER , I  3, 2 X , 5HCH0SE , F6. 3 ) 
PRINT  766, CONVRG 

FORMAT  (38HWHEN  OMEGA ( N4  1 )  -OMEGA(N)  IS  LESS  THAN.E14.8/ 
253HSTOP  COMPUTATION  AND  START  OVERRELAXATION,  IF  DESIRED! 
STOP 
END 
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C  ANALYTICAL  SOLUTION-3  DIMENSIONS,  CONSTANT  SOURCE 

C 

C  SPECIFY  DIMENSIONED  VARIABLES 

DIMENSION  A  I  30 1 ,B ( 30  I ,C ( 30 ) 

C 

C  READ  INPUT  DATA 

READ  INPUT  TAPE  2,60, LMAX, MM AX ,NMAX 

60  FORMAT  1313) 

READ  INPUT  TAPE  2 , 6 1 , DI FUS , SORC , COND, XMAX , YMAX, 2MAX 

61  FORMAT ( 6F 10  «0 ) 

READ  INPUT  TAPE  2 , 62 , DE LX.DELY , DEL* ,DEL T 

62  FORMAT  I 4F 10*0  I 

READ  INPUT  TAPE  2 , 63 , TMAX, XGR1 D, YGR ID, ZGR I D 

63  FORMAT ( 4F 10  «0 ) 

C 

C  INITIALIZE  ANO  SET  DESIRED  TIME 

TIME«0. 

5  T IME  =  TIME«-DELT 

write  Output  tape  3,si,time 
51  FORMAT (/5X.I8HNEW  TIME  INCREMENT , F 10. 5 ) 

C 

C  PRINT  COLUMN  HEADER  FOR  RESULTS 

WRITE  OUTPUT  TAPE  3,53 

53  FORMAT ( //8X , 1HL, 3X, IHM, 3X , 1HN, 4X , 1 1HTEMPERATURE ,9X , 4HT  INC) 

C 

c 

C  COMPUTE  CONSTANTS 

P I »  3« 1415926 

0»64.«0IFUS»S0RC/IC0ND*Pl*PI*PI ) 

E«DIFUS»PI»PI 
DO  32  I*I»LMAX 

32  A{||t2*IHHl 
DO  33  J-l.MMAX 

33  BIJ)»2*(J-1)*I 
DO  34  X* 1 , NMAX 

34  ClK)»2*(K-l)*l 
C 

C  INITIALIZE  AND  SET  OESIREI)  X,Y,Z 
Z*0. 

35  Z»Z*DELZ 
Y«0. 

30  Y»Y*OELY 

X»0. 

8  X«X*OELX 

TEMP»0. 

C 

C  INITIALIZE  SUMMATION  VARIABLE  FOR  TEMPERATURE  INCREMENTS 

T  INC*0. 

C 

C  COMPUTE  TEMPERATURE  INCREMENTS 
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DO  10  N=  l , NMAX 
DU  10  M*  1 , MMAX 
1)0  10  L  =  l » LMAX 

F  =  A(U  *«2/XMAX»»2*B<M)*»2/YMAX**2*C(N)*»2/ZMAX»*2 
ALPHA=E*F 

D6N0M* A 1l)*H(M)*C(N)*AL PHA 

P=A<L)»PI*X/XMAX 

Q=B(M)*P1»Y/YMAX 

R*C(N)*PI*Z/ZMAX 

S*SINF(P)*$|NF(QJ*SINFIR) 

T*  AL PHA»  TIME 
C 

C  RESTRICT  EXPONENT  TO  PREVENT  SUBROUTINE  OVERFLOW 

I F ( T-30.  )  3,3,2 

2  V=1.0 
GO  TO  6 

3  V=l.-( l./EXPFIT)) 

6  TINC=D*S*V/DENOM 
C 

C  SUM  TEMPERATURE  INCREMENTS 

TEMP=T  EMP*T INC 
C 

C  PRINT  RE  SUL  1 S  WHEN  DO  LOOP  INDICES  ARE  EQUAL  FOR  MONITOR  ONLY 

1 F ( l-M) 10,7,10 

7  I F (K-N) 10,1,10 

1  1 =L- 1 

J  *M- 1 
K*N- 1 

WRITE  OUTPUT  TAPE  3 , 52 ,  I  , J , K, TEMP, T INC 
52  FORMAT  I  /5X,  3 1 4, 2 X, E 14. 8 ,2X, E 14. 8/ ) 

10  CONTINUE 

C 

C  PRINT  FINAL  RESULTS  OF  COMPUTATION 

WRITE  OUTPUT  TAPE  3,64 
64  FORMAT (  /5X, 1 3HFINAL  RESULTS/ I 

WRITE  OUTPUT  TAPE  3 , 50, X, Y , Z , TEMP 

50  FORMAT (5X,2HX*,F10.4, 4H  Y«,F10.4,4H  Z» , F 10 . 4, 3X , 5HTEMP« , E 14. 8/ ) 

CHECK  TO  SEE  IF  LAST  TIME  AND  SPACE  INCREMENTS 
HAVE  BEEN  REACHED 
IFUGRIO-X)  15,15,8 
15  1F(YGRID~*Y)20,20,30 

20  IF(ZGRI0-Z)25,25,35 

25  I F  t TMAX-TIME) 31, 31,5 

C 

C  PRINT  OUT  INPUT  DATA  FOR  REFERENCE  PURPOSES 
31  LL«LMAX-i 

MM«HMAX-l 
NN»NMAX-1 

WRITE  OUTPUT  TAPE  3,54,LL,MM,NN 
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54  FORMAT ( /5X, 25HMAX l  MUM  SUMMATION  INDEXES, 315) 

WRITE  OUTPUT  TAPE  3 , 55 , 0 l FUS, SORC 

55  FORMAT ( /5X, 1 IHDIFFUSI VI TY,F 10.4 ,3X, 15HS0URCE  STRENGTH, F 10.4 } 
Wft I f E  OUTPUT  TAPE  3,56,C0N0,XMAX,YMAX,ZMAX 

56  F0RMAT(/5X, 12HC0N0UCT  IVITY»F9.4,3X» 1 3HMAX I  MUM  X , Y, l , 3F 10. 5 ) 
CALL  EXIT 

END 
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C  HtAT  TRANSFER  IN  CUBE-NCNS TE AOY ,  CONSTANT  SOURCE  BY 

t  SUCCESSIVE  OVERREL AXAT I CN 

C 

01  MENS  ION  U!7,7,7),UOLD<7,7,7),DIFF{7,7,7) 

C 

C  RE AO  INPUT  DATA 

REAL)  61,  LAST X,L AS TY, LAST! 

61  FORMAT  1313) 

KEAD53,niFUS,COND,DELT,SORC 

53  FORMAT i 4F 10.0  ) 

READ  54.ERR.XGRI0, YGR l D ,2GR I D, TMAX 

54  FORMAT! 5F10.0) 

READ  56, 1 K I T6 , JR  I T E ,KR I TE , J I  TER, K 1  TER , CHOSE, BOUND 
56  FORMAT ( 514.2F10.0) 

READ  59,RFACT 
59  FORMAT  (E14.B) 

C 

C  SET  BOUNDARY  VALUES 

00  3  I*1,LASTX 
00  3  J  *  l ,  L  A  S  T  Y 
00  3  K*1,LAST2 
C 

C  THE  INPUT  BOUND  DETERMINES  METHOD  OF  SETTING  INITIAL 

C  ESTIMATES  OF  TEMPERATURE!  BOUND*  0.0  OR  LESS  AND  READ 

C  ESTIMATE  FROM  CARDS,  IF  BOUND  IS  GREATER  THAN  EERO  ALL 

C  POINTS  ARE  SET  TO  7ER0) 

IF l HOUND)  100,100,101 

100  READ  102 , U( I , J , K ) 

102  FORMAT  ( E  14. 8 ) 

UOLD! I , J,K)*U( I, J,K) 

GO  TO  3 

101  UOLD ( I , J , K ) *0.0 
Ul  I  ,  J,K)*0.0 

3  CONTINUE 

C 

C  COMPUTE  CONSTANTS 

A* ! 2. / ! XGRI0**2 ) )♦ ! 2. ! ! YGR I D»*2 ) ) ♦  ! 2. / ! ZGR I0#*2 ) ) 

B*D I FUS*DEL  T 
C*B/ ! XGR I D»  *2 ) 

D*tJ/  (  YGR  1 0*  »2  ) 
e«B/!2GKlD**2) 

G*l./I l.»B*A) 

F«G*B*SORC/CONO 

C 

C  SET  TIME  STEPS  AND  ITERATION  COUNTER 

T I  ME  *0. 

14  TIME*TIME*DELT 

I TNUM*0 
C 
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C  THE  INPUT  CHOSE  DETERMINES  WHETHER  THE  OUPUT  IS  PUNCHED 

C  OR  PRINTED  (SET  CHOSE  ■  0.0  OR  LESS  FOR  PUNCHEO  OUTPUT 

C  OR  GREATER  THAN  ZERO  FOR  PRINTED  OUTPUT) 

IF(CH0SE)B0,80,81 
BO  PUNCH  50, TIME 

GO  TO  1 

81  PRINT  50, T I  ME 

50  FORMAT ( /  IBHNEW  TIME  INCREMENT , F10.3 ) 

1  1 TNUM=I TNUM*1 
IF(C HO SF)B2, 82,83 

82  PUNCH  51,1 TNUM 
GO  TO  84 

83  PRINT  51,1 TNUM 

51  FORMAT  (/15X, L6H1TFRATICN  NUMBER, 14/) 

C 

C  COMPUTATIONAL  STEPS 

84  I LAST*L  AS  TX- l 
JLAST*LASTY-l 

klast*lastz-i 
DO  2  K*2,KLAST 
UO  2  Ja2 , JL AS T 
DO  2  I  s2 , ILAST 

P-G*C*(U( I*l,J,K)+U( 1-1 ,J,K) ) 

Q*G*0<MU(  I,  J*1,KI*UU,J-1,KI ) 

R*G*EMU( 1, J,K*1 )*U( I,  J,K-1 )) 

S=G*UOLO(l,J,K) 

UNEW=(  (F*P»0*R*S)*RFACT  )M(1.0-RFACT  )*U( I , J,K) ) 

OIFF(I*J»K) *UNEW-UC I , J, K ) 

C 

C  COMPLETE  OR  PARTIAL  OUTPUT  OF  TEMPERATURES  CAN  BE  OBTAINED 

C  FOR  EACH  ITERATION  BY  PROPER  CHOICE  OF  THE  INPUTS  JITER  AND  KITfcR 

116  IF(JITER-J) 110,110,2 

110  IF(KITER-K)8,8,2 

8  I F ( CHOSE ) 85 , 85 ,86 

85  PUNCH  52,1, J ,K,UNEW,DIFF  1 1 , J , K ) 

GO  TO  2 

86  PRINT  52,I,J,K,UNEW,DIFFU.J,K) 

52  FORMAT! 10X, 314, 2E14. 8) 

2  U< I,J,K)*UNEW 
C 

C  SET  VALUES  ONE  GRID  POINT  BEYOND  PLANES  OF  SYMMETRY 

L*LASTX 

DO  4  N-2.LASTZ 
00  4  M»2,LASTV 

4  U(L,M,N)*UIL-2,M,N) 

M-LASTY 

00  5  L*2 ,LASTX 
00  5  N*2 , LASTZ 

5  U(L,M,N)«U(L,M-2,N) 

N-LASTZ 
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on  6  L=2,LASTX 
00  6  M=2,LASTZ 

6  U(L,M,N)=U(L,M,N-2) 

C 

c  re S I  SOLUTIONS  against  established  criteria 

00  10  K-2  *  KL AS  T 
00  10  Ja2 » JLAST 
00  10  1=2,  H  AST 

IF ( ABSF (DIFFI 1 • J , K ) J-ERR )  10,10,1 

10  CONTINUE 

C 

C  PRINT  RESULTS  OF  ITERATION 

IF (CHOSt 187,87,88 

87  PUNCH  60 
GO  TO  89 

88  PRINT  60 

60  FORMAT  ( // 1  OX , 26HF INAL  RESULTS  OF  ITERATION//) 

89  DO  /  K-2.KI AST,KRI TE 
00  7  J=2, JLAST, JRlTfc 
1)0  /  I=2,ILAST,IRITfc 
IFICH0St)90,90,91 

90  PUNCH  52, I , J , K , U ( I , J , K ) ,DIFF( I , J , K  ) 

GO  TO  7 

91  PRINT  52«I»J,K,U(I,J,K),DIFF(I,J,K) 

7  CONTINUE 

C 

C  CHECK  TO  SEE  IF  LAST  TIME  INCREMENT  HAS  BEEN  REACHED 
IF(TMAX-TIME)17,17,12 

C 

C  PREPARE  FOR  A  NEK  T IMF  INCREMENT 

12  00  15  1*2, LASTX 

DO  15  J=2 , L AS T Y 
DO  15  K=2,l ASTZ 

15  UOLDI I , J,K)»U(I , J,K) 

GU  TO  14 

C 

C  PRINT  OUT  INPUT  OATA 

17  IFICh0SF)93,93,94 

93  PUNCH  55,DIFUS,C0ND,UEIT,S0RC 

PUNCH  57,XGRI0, YGRID.ZGRIO, TMAX.ERR 
PUNCH  53,IRITE,JR1TE,K(UTE,KITER,JITER,CH0SE 
PUNCH  62,  RF ACT 
GO  TO  13 

94  PRINT  55, 01 FUS ,CONO,OELT , SORC 
55  FORMAT ( I IHD I FFUS IVl TY,F 10.4,6X, I2HC0N0UCI I Vl TY.F10.4/ 

310HDEL  T  A  T I  ME , F 10. A , 6X , 15HS0URCE  STRENGTH, FI  0 . A ) 

PRINT  5 7 , XGR IO,YGRID»ZGRID»TMAX,ERR 

57  FORMAT ( /6HXGR ID> ,F10.4,6X,6HYGR10*,F10.4,6X,6HZGRID*,F10.4/ 

48HMAX  T1ME,F10.4,18HSTOP  WHEN  ABSOI FF-,E 14.8 , 16H  IS  ZERO  OR  LESS) 
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PAINT  58,1  RITE, JR  I TE ,KR I TE , KI TER . J l TER, CHOSE 
58  FORMAT ( /25HOUTPUT  CONTROLS  ARE  I R I TE , 1 3 , 2X ,5HJR I TE , 1 3/ 
55HKRI TE , I  3, 2X , 5HK ITER, I  3 , 2X , 5HJ I  TER , 1 3, 2 X , 5HCHOSE , F6. 3 ) 
PRINT  62,  RFACT 

62  FORMAT (/22HRELAXATION  FACTOR  USED, 3X , E 14 . 8 > 

13  STOP 
ENO 
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C  He  A  T  TRANSFER, CONSTANT  INTERNAL  SOURCE,  TWO  DIMENSIONS 
C  SOLUTION  UY  GAUSS  SIcOEL  WITH  SERIAL  SCAN 
C 

DIMENSION  U(4  W41  )  ,UOLO  C  4 1 , 4 1  )  ,  Q  l  FR  (  4  1  <  4  1  > 

Re  AO  Input  TARE  2 , 700 ,LASTX .LAST Y , I R I TE , JRI TE , JI TEK 
700  FORMAT  (514) 

READ  INRUT  TARE  2 , 70  I , 01 FUS , COND , SORC , XGR 1 0 , YGR I D 
7C 1  FORMAT  ( 5F 12.01 

READ  I N  R  lj  1  TAPE  2,702,1  MAX, CELT, ERR 
702  FORMAT  (if  12.0) 

READ  INPUT  TAPE  2 , 70  3  .BOUND  ,  P  I  TER 
70)  FORMAT  (2F1C.0) 

C 

C  PRINT  UUT  INPUT  DATA  FOR  REFERENCE  PURPOSES 

WRITE  OUTPUT  TAPE  3 , 53 , DEL T , SORC 

S3  FORMAT  (11W0E4.TA  T  IMC  =  ,F  1  0.4.6X,  UHSOGftCE  S  T8ENGTH*  ,  F  10 .  */  ) 
WRITE  OUTPUT  TAPE  3, 55,DIFUS,CONQ 
51*  FORMAT  (  /12HDI  FFUSI  VI  TYa,Fl0.4,6X,  13HCCN0UCTI  VI  TYstF  10.4/  ) 

WRITS  OUTPUT  TAPE  3,56 

56  FORMAT ( U  7UTUIS  RUN  IS  20  GAUSS-S 1E0EI/ > 

WRITE  OUTPUT  TAPE  3, 57, XGKI D.YGRI D, TMAX 

57  FORMAT  ( 6HXGR ID- ,F10.4,3X , 6HYGR IO=,F10.4, 3X , 9HMAX  T I ME  =  , FI  0 . 4 / ) 
WRITE  OUTPUT  TAPE  3 , 58,  IR  I  TE , JR  I TE , J I  TER 

58  FORMAT  (  /  6H IU I  T  E  =  ,  I  4.8H  J!UTE-^I4,5H  JJIJER=,14/) 

C 

C  IF  DOUNO  -  C.O  OR  -1.0, ROAD  INITIAL  VALUES  AS  INPUT  DATA,  OR  IF 

C  BOUND  =  1.0, SET  ALL  POINTS  EQUAL  ZERO. 

DO  3  J=  l  , l AS  TY 
DO  3  1  =  1,1. ASTX 
IF  (SOUND)  100,100,101 

100  REAP  INPUT  TAPE  2,704,U(I,J) 

704  FORMAT  'F12.C) 

pir<(  i ,  i)  =o.o 

UUL  0 ( 1 , J)=U< 1 , J) 

(.1!  TU  3 

101  UOLD ( I , J ) =0.C 
n i f p ( i,j)=c.u 

U( 1. J  1=0.0 

3  CONTINUE 

n 

G 

C  COMPUTE  CONSTANTS 

ILAST=LASTX-1 
JLAST=LASTY-1 

201  A= ( 2./ I XGR I D» *2 ) )F(2. /(YGR 10**2)) 

B=P1FUS*DELT 

C=B/(XGR10«*2) 

D=H/ ( YGR 10**2) 

G=I./U.*B*A) 

F=G*B*SORC/CONO 

C 

C  SET  TIME  STEPS  AND  ITERATION  COUNTER 

T I MF=0. 

14  T I MC  =  T I  ME  *DEL  T 
ITNUM*0 

WRITE  OUTPUT  TAPE  3, 50, TIME 

184. 


•A/PHIS/43-8 

50  FORMAT ( /18HNEW  TIME  INCREMENT  *  F 10 . 3 ) 

1  I  TAUJKsi  TNUK*1 
C 

C  WHEN'  PRIN11NG  OUT  DATA  EACH  ITERATION  SET  PITER  =0.0,0«  -UO 
C  OTHERWISE  SET  PITER  =  1.0. 

IF  (PITER?  304,306,84 
306  WRITE  OUTPUT  TAPE  3,51,ITNUM 

51  FORMAT  {/15X, 16H1TFRAT10N  NUMBER,I4/J 
C 

C  COMPUTATIONAL  STEPS 

84  DO  2  J  =  2 » JL  AS  T 
UO  2  1=2, ILAST 
P=G«C*(U( I+I, J)+U(I-I,J) ) 

0=G*0* (Ul  I,J+1)+U(I*J-1)) 

S=0*U0LUU  ,  J) 

UNEW=F*P*C+S 

OIFR ( I , J ) =UNE W— U ( I , J ) 

116  I F ( J I  TER- J )  110,110,2 

110  WRITE  OUTPUT  TAPE  3,52,1 , J , UNE W , 0 1 FR ( I , J ) 

2  U  l  l ,  J  )  =UfiLW 
C 

C  CHECK  CONVERGENCE  OF  SOLUTION  AGAINST  ESTABLISHED  CRITERIA 

24  00  10  J=2,JLAST 

00  10  1=2, ILAST 

IF(ABSF(U1FR( I ,J)  )«ERR)  10,10,1 
10  CONTINUE 

C 

C  PRINT  OUT  FINAL  RESULTS  OF  ITERATION 

WRITE  OUTPUT  TAPE  3,60 
WRITE  OUTPUT  TAPE  3,51,ITNUK 
WRITE  OUTPUT  TAPE  3, 54, ERR 

54  FORMAT  l  /34HSTCP  ITERATION  IF  ABS  (I(h!*l)-T  (N)- ,  E14.0 , 10H*0  OR  LESS) 

60  FORMAT  (//3IH  FINAL  RESULTS  OF  ITERATION//) 

WRITE  OUTPUT  TAPE  3,602 

602  FORMAT  (46H  I  J  TEMPERATURE  DI FF , T ( N+l )- T ( N ) ) 

89  CO  7  J=1  ,LAST Y, JRITC 
DO  7  1*1, LAST*, IRITE 

7  WRITE  OUTPUT  TAPE  3 , 52, I , J,UIJ , J ) ,DI FR ( l , J ) 

52  FORMAT  (214 , 4X , E 1 5. 8 , 4X , E 15. 8 ) 

C 

C  IF  DESIRED, RE  AO  IN  NEW  CONVERGENCE  CRITERIA  AND  START  WHERE  LEFT  OFF 

REAO  INPUT  TAPE  2, 705, ERR 
705  FORMAT  (FIO.O) 

IF ( 100. -ERR )  302,300,300 
300  WRITE  OUTPUT  TAPE  3, 305, ERR 

305  FORMAT  ( //27HNE.W  CONVERGENCE  CRITERIA  *  ,EI4.8//) 

GO  TO  l 
C 

C  CHECK  TO  SEE  IF  LAST  TIME  INCREMENT  HAS  BEEN  REACHED 
302  IF ( TMAX-T IM£)  17,17, L2 
12  DO  15  I*2,LASTX 
OU  15  J*2,LASTY 
15  UOLD ( I , J ) *U ( I , J ) 

GO  TO  14 
C 

17  CALL  EX  1 1 
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C  HC  A  T  TRANSfcR.CONSTANT  INTERNAL  SCURCE,  TWO  DIMENSIONS 

C  SOLUTION  BY  GAUSS  SIEHEl  WITH  SPIRAL  SCAN 

DIMENSION  0(41,41) .HOLD (41, 41) ,DI FR(41 ,41  ) 

READ  INPUT  TAPE  2 , 700 ,L AS TX ,1  AS T Y , I R I TE , JK I T F , J I  TER 
70C  FORMAT  (514) 

READ  INPUT  TAPE  2 , 701 , 0 1 FUS ,COND, SORC , XGR 10 , YGR 1 0 

701  FORMAT  (5F12.0) 

READ  INPUT  TAPE  2, 702  ,  TMAJC,  DEL  T  ,  ERR 

702  FORMAT  (3F12.0) 

READ  INPUT  TAPE  2 , 703 .BOUND , P I  TER 

703  FORMAT  (2F10.0) 

C 

C  PRINT  OUT  INPUT  DATA  FOR  REFERENCE  PURPOSES 

KRITc  OUl PUT  TAPE  3, 53.DELT »SURC 
53  TORMAT  (11HDELTA  T I  ME = , F 1 0 . 4 , 6X , 1 6HSOORCE  STRENGTH*, F 10.4/ ) 
WRITE  OUTPUT  TAPE  3 , 55, D I FUS ,COND 

55  FORMAT (/12HCIFFUS I VI TY*,F 10.4, 6X, 13HCCNDUCT I VI TV*,F10.4/) 

WRITE  OUTPUT  TAPE  3,56 

56  FORMAT ( /44HTHIS  RUN  IS  2D  GAUSS-SIEDEL  K I TH  FULL  SPIRAL/) 

WRITE  OUTPUT  TAPE  3, 57, XGAID, YGR Ifi.TMAX 

57  FORMAT  ( 6HXGR 1 0* ,F 10. 4, 3X .6HYGR1 D* , F 1 C . 4, 3X , 9HMAX  T I  ME* , FI  0 .4/ ) 
WRITE  OUTPUT  TAPE  3, 56 , I R I TE , JR  I TE , J I  TER 

58  FORMAT  ( / 6H I R I T E= , 1 4 , 8H  JRIIC  =  ,I4,8H  JITEK«,  14/ ) 

C 

C  IF  BOUND  =  C.O  OR  -l.O.READ  INITIAL  VALUES  AS  INPUT  DATA,  OR  IF 

C  BOUND  =  1.0, SET  AU.  POINTS  EUUAL  ZERO. 

DO  3  J*  1 , LASTY 
DU  3  1  =  1, LAS  T  X 
IF  (BOUND)  100,100,101 

100  READ  INPUT  TAPE  2,704,U(I,J) 

704  FORMAT  (F12.C) 

UIERll, J)=C.G 
UULD ( I , J ) *0 ( 1 , J ) 

GO  TO  3 

101  ooint i , j) =o.o 
0  I  FP.  I  I  ,  J  T  =0.0 
U( I  ,  J  )  =  0 . 0 

3  CONTINUE 
C 

C  COMPUTE  CONSTANTS 

1LAST=IASTX-1 
JL  AS  F  =1 AS  I Y- 1 

201  A* ( 2. / ( XGR 1 U»  *2 ) ) ♦ ( 2 . / ( YGR I D**2 ) ) 

B=D1FUS»DEL T 
C*B/(XGR1D**2) 

D*B/ I YGR 10**2) 

G  =  I  . / ( 1 . ♦B*A ) 

F=G«B*SORC/CONO 

C 

C  SET  TIME  STEPS  AND  ITERATION  COUNTER 

T1ME=G. 

14  TIME=  TIME  +OEL  T 
I TNUM=0 
C 

WRITE  OUTPUT  TAPE  3, 50, TIME 
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50  FORMAT ( / 1 8HNE W  TIME  I NCREMENT , F1P, 3 ) 

L  I  TNUAU1  TNUM>1 

C 

C  WHEN  PRINT  I NG  OUT  CAT  A  EACH  14-ERA-T10N  SET  PlI-ER  =Rl.U,Oft  -UO 

C  OTHERWISE  SET  PITER  =  1.0. 

IF  (PITER)  3C6, 306,84 
306  WRITE  OUTPUT  TAPE  3,5l,lTNUM 

51  FORMAT  1/15X, 16HIJERAT10N  NUM6ER, 14/) 

C 

C  COMPUTATIONAL  STEPS 

84  11*1 

J  Js2 

1 1 1  =  I  LAST 
J JJ= 1LAST 
CCUNT=0.0 

30  COUNT *CQUNT +1 .0 

IF (COUNT- 1.0)  31  , 31, 32 

31  J=JJ 

1 1*  1 1  ♦  1 
GO  TO  34 

32  I  F  ( COUNT-2 . G )  33,33,36 

33  III=IJI-1 

45  J=JJJ 

34  DO  35 

P=G*C* ( U{ I+l, J)+U(  I-1,J)  ) 

0=G»D*(U(1,J-*L)+U(  l,J-l) ) 

S=G*UOLO( I , J) 

UNEW=F*P+U+S 

OIFR( I , J)=UNEW-U( I , J) 

116  IFUITER-J)  110,110,35 

110  WRITE  OUTPUT  TAPE  3, 52 , 1 , J ,UNEW, 0 IFR  I  l  ,  J  ) 

35  U ( 1 , J ) =UN£W 

IF(  I  I  —  III )  46,29,46 

46  GO  TO  30 

36  1 F ( COUNT- 3. 0 )  37,37,38 

37  I  =  1 1 

JJ  =  JJ*l 

JEND= J J J 
GO  TO  39 

38  I*III41 
JEND*JJJM 

39  00  44  J*J J , JENO 
p=G*c*(ud+i,j)*ui  i-i,  on 
0*G*0*(U( I, J+1)*U(  I , J- 1 )  ) 

S*G*UOLO(I, J) 

UNEW*F4Ptli4S 

DIFR ( I , J ) *UNE W-U( 1 , J ) 

40  I F ( J I  TER- J )  41,41,44 

41  WRITE  OUTPUT  TAPE  3,52,1  ,J,UNE«,DIERU,J> 

44  U ( I , J ) *UN£W 

l F( COUNT- 3,0  30,30,47 

47  COUNT*0.0 
GO  TO  30 

C 

C  CHECK  CONVERGENCE  OF  SOLUTION  AGAINST  ESTABLISHED  CRITERIA 
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29  DO  10  J  =  2,JLAST 
OU  10  1=2#  JLASI 

IF(AQSF(01FR( I,J) )-ERR)  10,10,1 
10  CONTINUE 
C 

C  PRINT  OUT  FINAL  RESULTS  OF  ITERATION 
WRITE  OUTPUT  TAPE  3,60 
WRITE  OUTPUT  TAPE  3,51,1TNUM 
WRITE  OUTPUT  TAPE  3, 54, ERR 

54  FORMAT! /  34HST0P  ITERATION  IF  AP.S  (  T  I N*  1 )  -  T I  N )- ,  E  14.8 , 1  0H*C  OR  LESS) 

60  TORMAT  ( / / 3 1H  FINAL  RESULTS  OF  ITERATION//) 

WRITE  OUTPUT  TAPE  3,602 

602  FORMAT  (46H  1  J  TEMPERATURE  0 1 FF , T ( N+ 1 ) -T ( M ) 

89  DO  7  J=1»LASTY,JKITE 
DO  7  1=1 ,1 ASTX, IRITF 

7  WRITE  OUTPUT  TAPE  3 , 52 , 1 , J ,U ( I , J ) ,0 1 FK ( 1 , J ) 

52  FORMAT  ( 2 l 4 , 4 X , E 1 5 . 8 , 4X ,£ 1 5 . 8 ) 

C 

C  IF  DESIRED, READ  IN  NEW  CONVERGENCE  CRITERIA  AND  START  WHERE  LEFT  OFF 

READ  INPUT  TAPE  2, 705, ERR 
705  FORMAT  (FlU.G) 

I F ( 100. -ERR)  302,300,300 
30C  WRITE  OUTPUT  TAPE  3, 305, ERR 

30  j  FORMAT ( //27HNEW  CONVERGENCE  CRITERIA  =  ,E14.8//) 

GO  TO  1 
C 

C  CHECK  TO  SEE  IF  LAST  TIME  INCREMENT  HAS  BEEN  REACHED 

302  IF ( TMAA-T IKE )  17,17,12 
12  DO  15  I=2»LASTX 
00  15  J  =  2  »  L AS  TY 
15  UGLD ( I ,1 ) =U ( I , J ) 

GO  IU  14 
17  CALL  EXIT 

E  MO (lfGfOfOfOiCt I#C»0flf0t0f0#G»0) 
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OVEKREIAXATION  DATA  RUN  TIME  *  0.2  SECtH  »  O.I  CM, ERR  »  0.005 


NEW  TIMC  INCREMENT  .200 

ITERATION  NUMBER  I 

node  -temper*  bure  residual 

2  6  6  . 1 183S672E+02  .118356726+02 

366  . 18639951E+02  . 186399516+02 
466  .224053656+02  . 22405365E+02 
5  6  6  . 24404534E  +02  . 24404534E  +  02 

666  .25423319E+02  .254233I9E+02 

ITERATION  NUMBER  2 

266  .19837496E+02  .800182406+01 
366  .31652511E+02  . 130 12560E+02 
4  6  6  « 38395C91E+02  . 15989726E+0 2 

566  .42103991E+02  . 17699457E+02 
666  .43346459E+02  . 17923140E+02 

ITERATION  NUMBER  3 

266  .25360193E+02  . 55226970E+01 
366  .40854991E+02  . 92024800E+01 
466  .50007493E+02  . 11612402E+02 
566  .54920496E+02  . 12816505E+02 
666  . 56663 7616+02  . 13317302E+02 

ITERATION  NUMBER  4 

266  .29337214E+02  . 39770210E+01 
366  .47617934E+02  .67629430E+01 
466  . 58493375E+02  . 84858820E+01 
566  . 64348760E+02  .94282640E+01 
666  .66207560E+02  . 95437990E+01 

ITERATION  NUMBER  5 

266  . 32258438E+02  .29212240E+01 
366  .52541071E+02  .49231370E+01 
466  .64693040E+02  .61996650E+01 
566  .71143600E+02  .67948400E+0I 
666  . 732774 10E+02  . 70698500E+01 

ITERATION  NUMBER  6 

266  .34390562E+02  .21321240E+01 
366  .56157207E+02  . 36161 360E+01 
466  .69202630E+02  .45095900E+01 
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566  . 76175C30E+02  . 50314300E*01 
666  .78357930E+02  . 50805200E+01 

ITERATION  NUMBER  7 

266  • 35958S67E+02  . 15684050E+01 
366  • 58782838E+02  . 26256310E*01 

4  6  6  . 725 17520E ♦02  . 33148900E+01 

566  . 797895 10E+02  . 36148800E*0I 
666  .82129770E+02  . 377 18400E+01 

ITERATION  NUMBER  8 

2  6  6  . 37093 l 14E*02  .11 341470E+01 

3  6  6  .60712752E+02  . 19299 140E+01 

466  .74915920E+02  . 239840006*01 
566  .82479230E+02  . 268932006*01 
6  6  6  .84836000E*02  . 270623006*01 

ITERATION  NUMBER  9 

266  .37930510E*02  . 83739600E*00 
366  .62109970E*02  . 139721806*01 
466  .76687140E*02  . 17712200E*01 
566  .843986506*02  . 191942006*01 
666  .863477006*02  .201 17000E*01 

ITERATION  NUMBER  10 

2  6  6  .385328386*02  .602328006*00 

366  «63139090E*02  . 102912006*01 
466  • 77958190E*02  . 12710500E*01 

5  6  6  ,85836100E*02  . 143745006*01 

6  6  6  .88286590E*02  .  143889006*01 

ITERATION  NUMBER  11 

266  .389802826*02  .447444006*00 
366  .638808806*02  . 741 79000E+00 
466  . 789081406*02  .949950006*00 
566  .868551706*02  .101907006*01 
6  6  6  .893611706*02  . 107458006*01 

ITERATION  NUMBER  12 

266  .39299367E*02  .319085006*00 
366  .644326606*02  . 551 780006*00 
466  .795802006*02  .672060006*00 
566  .876236406*02  .768470006*00 
666  .901 252106*02  . 764040006*00 
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ITERATION  NUMBER  13 

2  6  6  . 395*02 7 IE *02  . 24090400E+00 

366  .64824630E+02  . 39197000E+00 
466  .80090450E+02  . 51025000E+00 
566  • 88164880E+02  . 541 24000E+00 
666  «90700150E+02  . 57494000E+00 

ITERATION  NUMBER  14 

266  .39707595E+02  . 16732400E*00 
366  »65121470E+02  . 29684000E+00 
466  .80445510E«02  . 35506000E+00 
566  • 88575520E+02  .41064000E+00 
6  6  6  .91105250E+02  . 405 lOOOOE+OO 

ITERATION  NUMBER  15 

2  6  6  .39838593E*02  . 13099800E  +  00 

366  *65328Q00E+02  . 20653000E+00 
466  • 80719400E+02  .27389000E+00 
566  .88862960E+02  . 28744000E+00 
666  •91413090E+02  . 30784000E+00 

ITERATION  NUMBER  16 

266  • 39925239E>02  . 86646000E-01 
366  • 65487 780E *02  . I5978000t*00 
466  ,80907200E*02  . 18780000E+00 

5  6  6  ,89082480E*02  . 2 1952000E+00 

6  6  6  .91627760E+02  . 2 1467000E»00 

ITERATION  NUMBER  17 

266  • 39997249E«,02  . 72010000E-01 
366  •65596810E >02  . 10903000E+00 
466  .8l054l30fc*02  . 14693000E«00 
566  »89235080E+02  . 15260000E*00 
666  .91792690E*02  . 16493000E+00 

ITERATION  NUMBER  16 

266  •40041580E+02  . 44331000E-01 
366  .65682530E*02  . 85720000E-01 
466  .8ll53440E*02  . 993 lOOOOE-OI 
566  •89352440E+02  . I 1736000E*00 
666  .91906360E«02  . I 1367000E*00 
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ITERATION  NUMBER  19 

266  .4008L412E+02  . 398 32000E-01 
366  .65740310E*02  . 57780000E-0I 
466  .81232J00E+02  . 78860000E-01 
566  .89433450E+02  . 81010000E-01 
6  6  6  • 91994 770E ♦02  . 88410000E-01 

ITERATION  NUMBER  20 

266  .40103881E*02  . 22469000E-0I 
366  .657861L0E+02  . 45800000E-01 
466  .8I284740E+02  . 52440000E-01 

5  6  6  . 89496200E *02  . 62750000E-0I 

666  .92054930E«02  . 60160000E-01 

ITERATION  NUMBER  21 

2  6  6  .40125995E+02  . 22 1 14000E-01 

366  .65816670E+02  . 30760000E-01 
466  .8I327I40E*02  .42400000E-01 
566  .8953920QE«-02  . 43000000E-01 
666  .92102350E+02  . 47420000E-01 

ITERATION  NUMBER  22 

2  6  6  .60137283E *02  . I I288000E-01 

366  .65841240E+02  . 24370000E-01 
466  .81354770E+02  . 2 7630000E-01 
566  .89572 760E*02  . 33560000E-01 
666  .92134160E+02  . 318I0000E-0I 

ITERATION  NUMBER  23 

266  .40149585E+02  . I2302000E-01 
366  .65857690E+02  . 16450000E-0! 
466  .81377620E*02  . 22850000E-0! 
566  . 89595570E+02  . 228 10000E-01 

6  6  6  .92159600E*02  . 25440000E-01 

ITERATION  NUMBER  24 

266  .40155211E+02  . 56260000E-02 
366  .65870600E+02  . I2910000E-0I 
466  .81392120E+02  . 14500000E-01 
5  6  6  .89613520E4-02  . 17950000E-01 
666  .921 76410E+02  . 16810000E-0I 
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ITERATION  NUMBER  25 

266  .40162C44E+02  . 68330000E-02 

3  6  6  « 65879430E+Q2  . 88300000E-02 

4  6  6  . 8  140449QE+02  .  12370000E-01 

5  6  6  . 89625630E+02  . 12 l lOOOOE-Ol 

666  .92190C80E*02  . 13670000E-01 

ITERATION  NUMBER  26 

266  . 40164849E+02  . 28050000E-02 
366  .65886260E+02  . 6B300000E-02 
466  .81412050E+02  . 75600000E-02 
566  • 89635220E*02  . 95900000E-02 

6  6  6  .92198960E*02  . 88800000E-02 

ITERATION  NUMBER  27 

266  .40168617E+02  . 37680000E-02 
366  .65891COOE+02  .47400000E-02 
466  .81418780E*02  . 67300000E-02 
5  6  6  •  8964 1660E  *-02  .  64400000E-02 

666  «92206290E*02  . 73300000E-02 

ITERATION  NUMBER  28 

266  .40170C28E+02  . 141 10000E-02 
366  .65894620E*02  . 36200000E-02 
466  •81422690E+02  . 39100000E-02 

5  6  6  . 89646 760E *02  .  5 1000000E-02 

666  .92210980E*02  .46900000E-02 

ITERATION  NUMBER  29 

266  «40172096E+02  .20680000E-02 
366  .65897 150E+02  .25300000E-02 
466  .8 1426350E+02  . 36600000E-02 
566  .896502 10E+02  . 34500000E-02 

6  66  .92214920E*02  . 39400000E-02 


FINAL  RESULTS  OF  ITERATION 

node  temperature  residual 

222  ,11337487E*02  . 22420000E-02 
322  .16195678E+02  . 19200000E-03 
422  . 18495888E+02  . 20890000E-02 
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5 

6 
2 

3 

4 

5 

6 
2 

3 

4 


2 

2 

3 

3 

3 

3 

3 

4 
4 
4 


2  . 19558  /23E ♦02 
2  .  19871819E  *02 
2  . 16195676E+02 
2  . 2407684 7E  *02 
2  .28071  184E *02 
2  .29983233E+02 
2  . 30552  b  16E  *02 
2  . 184958876*02 
2  .280  71 1H6E+02 
2  • 33 1 37749E ♦02 


.  6 100Q00QE-04 
•  16750000E-02 
.  1 8700000E-03 
.79000000E-03 
. 779000006-03 
. I1780000E-02 
. 94500000E-03 
. 20860000E-02 
. 77900000E-03 
.23800000E-02 


5  4 

6  4 

2  5 

3  5 

4  5 

5  5 

6  5 

2  6 

3  6 

4  6 

5  6 

6  6 

2  2 

3  2 

4  2 

5  2 

6  2 

2  3 

3  3 

4  3 

5  3 

6  3 

2  4 

3  4 

4  4 

5  4 

6  4 

2  5 

3  3 

4  5 

5  5 

6  5 

2  6 

3  6 

4  6 

5  6 

6  6 

2  2 


2  . 35623870E>02 
2  . 363735006*02 
2  . 19558  725E  +02 
2  . 2998323 IE  +02 
2  . 35623672E  >02 
2  *  38433894E  *02 
2  . 3928601 7E *02 
2  . 1987182 IE  +02 
2  .30552516E+02 
2  . 36373503E+02 
2  . 392860 1 7E ♦02 

2  .401720946+02 

3  . 16195675E+02 
3  .24076848E+02 
3  .2807 i 187E  +02 
3  • 2998323 IE >02 
3  . 305525206*02 
3  . 24076648E+02 
3  . 37132 153E  *02 
3  .44140288E+02 
3  .47593100E+02 
3  .48633861E+02 
3  .28071 1876*02 
3  .44140283E+02 
3  . 530913836*02 
3  .57599982E+02 
3  .58972680E+02 
3  . 29983232E *02 
3  .475931056*02 
3  .575999816*02 
3  .627059806*02 
3  .642705306*02 
3  • 305525176*02 
3  .48633858E+02 
3  . 589726766+02 
3  .642705306+02 

3  .65897  1406*02 

4  . 18495887E+02 


.  942000006-03 
.21 3600006-02 
.61000000E-04 
.11 720000E-02 
.94600000E-03 
. 174400006-02 
. 12150000E-02 
. 16760000E-02 
. 94400000E— 03 
.21390000E-02 
. 121300006-02 
. 206400006-02 
. 188000006-03 
. 79200000E-03 
. 78200000E-03 
. 1 1730000E-02 
. 95200000E-03 
.79300000E-03 
. 12610000E-02 
. 16250000E-02 
. 17370000E-02 
. 184 50000 E -02 
. 78000000E-03 
. 16160000E-02 
.188400006-02 
.22970000E-02 
. 21910000E-02 
.117300006-02 
. 1 7460000E-02 
.22940000E-02 
. 24100000E-02 
. 25800000E-02 
. 94400000E-03 
.18410000E-02 
•21850000E-02 
. 25800000E-02 
.25101)0006-02 
. 20900000E-02 
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3  2 

4  2 

5  2 

6  2 

2  3 

3  3 

4  3 

5  3 

6  3 

2  4 

3  4 

4  4 

5  4 

6  4 

2  5 

3  5 

4  5 

5  5 

6  5 

2  6 

3  6 

4  6 

5  6 

6  6 

2  2 

3  2 

4  2 

5  2 

6  2 

2  3 

3  3 

4  3 

5  3 

6  3 

2  4 

3  4 

4  4 

5  4 

6  4 

2  5 

3  5 

4  5 

5  5 

6  5 

2  6 

3  6 

4  6 

5  6 


4  .280711856+02 
4  . 33137748E+02 
4  . 35623870E +02 
4  .363735006+02 
4  .28071 185E+02 
4  .441 40289E  +02 
4  . 53091 380E+02 
4  .575999856+02 
4  . 58972683E+02 
4  . 33137744E+02 
4  • 53091 387E+02 
4  .64588  170E  +02 
4  . 70494550E+02 
4  . 72310C30E+02 
4  .35623874E+02 
4  . 57599976E+02 
4  . 70494550E  +02 
4  . 7 7198940E  +02 
4  . 79270700E+02 
4  • 3637349BE  +02 
4  . 58972678E  +02 
4  • 72310C20E  +02 
4  • 79270690E  +02 

4  . 81426340E+02 

5  • 19558725E  +02 
5  .29983231E+02 
5  .356238706+02 
5  • 38433898E+02 
5  . 39286C16E +02 
5  .29983233E+02 
5  .47593105E+02 
5  .575<»9985E+02 
5  .62705970E+02 
5  . 64270530E+02 
5  . 35623874E+02 
5  .57599977E+02 
5  • 70494550E+02 
5  .77198940E+02 
5  .79270700E+02 
5  • 38433898E  +02 
5  . 62 705980E  +02 
5  .771989406+02 
5  .848186706+02 
5  .871863806+02 
5  .392860186+02 
5  .642705306+02 
5  .792706806+02 
5  .871863706+02 


.77900000E-03 

•  23770000E-02 
• 94500000E-03 
.21 320000E-02 
. 779000006-03 
• 16230000E-02 
.187200006-02 
.229900006-02 

•  2 1930000E-02 
.237100006-02 
.188900006-02 
. 34500000E-02 
.252000006-02 
.348000006-02 
.951000006-03 
.22810000E-02 
. 25200000E-02 
. 32400000E-02 
.298000006-02 
.213000006-02 
. 21890000E-02 
. 34700000t-02 
.29700000E-02 
. 36400000E-02 
.630000006-04 

•  1 1730000E-02 
.944000006-03 
.175200006-02 
.121300006-02 
.117700006-02 
. 174300006-02 
.229900006-02 
.239000006-02 
. 25800000E-02 
.951000006-03 
.228400006-02 
.252000006-02 
. 32400000E-02 
.298000006-02 
.174900006-02 
.241000006-02 
.325000006-02 
.331000006-0? 
.363000006-02 
.121500006-02 
.258000006-02 
.295000006-02 
.362000006-02 
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665  .89650210E+02  .345000006-02 

2  2  6  • 1 98 7 1 82 IE  +02  . 16770000E-02 

326  .305525196+02  . 95000000E-03 
4  2  6  • 36373501E+02  .21 340000E-02 

526  .392860176+02  . 12160000E-02 
626  .40172C986+02  . 20700000E-02 
236  .305525176+02  .94400000E-03 
336  .48633865E+02  . 185 10000E-02 

4  3  6  *  589726  79E  *02  . 2 1 840000E-02 

5  3  6  .64270530E *02  . 258000006-02 

6  3  6  .65897  1 30E  +02  .249000006-02 

2  4  6  . 36373502E  +02  .213700006-02 

346  .589726836+02  .219300006-02 
446  . 72310C30E+02  . 34800000E-02 
546  • 79270700E+02  .298000006-02 
646  .81426340E+02  .364000006-02 
256  .392860186+02  .121500006-02 
356  .642705306+02  .258000006-02 
456  .792706806+02  . 29500000E-02 

5  5  6  .871 86370E  +02  . 36200000E-02 

6  5  6  . 896502 10E  +02  . 34400000E-02 

266  .40172C96E+02  .206800006-02 
366  .65897 150E+02  .253000006-02 
466  .814263506+02  . 36600000E-02 
566  .896502106+02  . 34500000E-02 

666  .92214920E+02  .394000006-02 

DIFFUSIV1TY  1.0000  CONDUCTIVITY  1.0000 

DELTA  TIME  .2000  SOURCE  STRENGTH  2000.0000 

XGRID*  .1000  YGR ID*  .1000  ZGRID*  .1000 

MAX  TIME  . 2000STOP  WHEN  A8SDIFF-  . 50000000E-02  IS  ZERO  OR  LESS 

OUTPUT  CONTROLS  ARE  IRITE  1  JRITE  1 
KR1TE  1  KITER  6  JITER  6  CHOSE  0.000 
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DATA  RUN  GAUSS  SIEDEL  T  =  0.2  SEC  ERR=0.005  H  *  0.1 

NEW  TIME  INCREMENT  .200 

ITERATION  NUMBER  1 

flode  temperature  residual 

266  .49352900E+01  .49352900E+01 
3  6  6  .61509471E+01  .61509471E+0] 

466  .64493826E+01  .64493826E+01 
566  .65222818E+01  .65222818E+01 
666  .65399740E+01  .65399740E+0 1 

ITERATION  NUMBER  2 

26  6  .88801862E+01  . 39448962E+0 1 

3  6  6  .  1 1 733575E+02  .55826280E+01 

466  .12595459E+02  .61460770E+01 
566  .12844798E+02  .63225170E+01 
666  . 12909294E+02  .63693200F+01 

ITERATION  NUMBER  3 

266  . 12 185228E+02  .33050420E+01 

3  6  6  .  16783346E+02  . 504977 1 0E  +  0 1 

4  6  6  • 18 385662E+02  .57902030E+01 
566  .1891 0285E+02  .60654870E+01 
666  .190521 05E+02  .614281 10E+01 

ITERATION  NUMBER  4 

266  . 15032849E+02  .28476210E+01 
366  .21351303E+02  .45679570E+01 
466  .23794226E+02  . 54085640E+0 1 

5  6  6  . 2467 1974E  +  02  .57616890E+01 

666  • 249168 19E+02  .58647140E+01 

ITERATION  NUMBER  5 

266  . 1 7527906E+02  .24950570E+0 1 
366  .2548fc291E+02  .4 1 349880E+0 1 
4  6  6  .28814 165E+02  . 50 1 99390E+0 1 

566  . 30098025E+02  .54260510E+01 
666  .30464459E+02  . 554 76400E+0 1 

ITERATION  NUMBER  6 
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266  .19736819E+02  .220891 30E+01 

3  6  6  .2923  1800E+02  « 3 7455090E+Q 1 

4  6  6  •  3  345  09  7 1 E  +  02  .46368060E+01 

5  6  6  .  35  171  236E+02  . 507 32 1 1 OE+O 1 

666  .  35670706E  +  02  • 52062470E+0 1 

ITERATION  NUMBER  7 

2  6  6  .21 705484E+02  •  1  9686650E  +  0 1 

366  . 32626222E+02  • 3 3944220E+0 1 
466  .37718101E+02  . 426 71 300E+0 1 
5  6  6  .39886695E+02  .  471 54590E  +  0 1 

666  .40524592E+02  • 485 38860E+0 1 

ITERATION  NUMBER  8 

266  .23467933E+02  • 1 7624490E+0 1 
366  .35703599E+02  . 307 73770E+0 1 
466  .41633842E+02  . 391 57410E+01 

5  6  6  .44248680E+02  .43621850E+01 

6  6  6  .4502 5903E  +  02  .4501 31 10E+01 

ITERATION  NUMBER  9 

266  .25050775E+02  . 1 5828420E+0 ] 
366  .38494323E+0 2  • 2 7907240E+0 1 
4  6  6  .452191 52E+02  . 35853 1 OOE+O 1 

566  .48268970E+02  . 402 00900E+0 1 
666  .49182403E+02  .41565000E+01 

ITERATION  NUMBER  10 

266  .26475626E+02  . 14248510E+01 

3  6  6  .4102 5641 E+02  . 253 1 3 1 80E+0 1 

466  .484962 04 E +02  . 32770520E+01 
566  .51 962594E+02  . 369 36240E+0 1 
666  • 53007382E+02  . 38249790E+0 1 

ITERATION  NUMBER  11 

266  .27760516E+02  . 12848900E+01 
366  .43322061E+02  .22964200E+01 

4  6  6  . 5 1 48 7 394E+02  .299 1 1900E+01 

566  . 55348087E+02  . 33854930E+01 
666  .565 1 7688E+02  . 351 03060E+01 
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ITERATION  NUMBER  12 

2  6  6  .28920777E+02  .  1 16 026 1 OE +  0 1 

3  6  6  .45405671E+02  • 206 36 1 OOE  +  O  1 

4  6  6  .54214692E+0 2  . 272 72980E  +  0  1 

5  6  6  *584452  13E+02  . 309 7 1 260E+0 1 

6  6  6  .59737212E  +  02  . 32 1 45240E  +  0  1 

ITERATION  NUMBER  13 

2  6  6  .29969620E+02  . 10488430E  +  0  1 

3  6  6  *472964 C6E  +  02  . 18907350E  +  0  1 

466  • 56699244E+02  . 24845520E+0 1 

5  6  6  .61274204E  +  02  . 282899 10E+0  1 

6  6  6  .62670778E+02  . 29385660E  +  0 1 

ITERATION  NUMBER  14 

266  .30918556E+02  .94893600E+00 
366  .4901 2283E+02  . 1 71 58770E+01 

4  6  6  .58961  106E+0?  . 226  1  8620E+0  1 

566  *638551 23E+02  . 268091 90E+0 1 
666  .65353378E+02  .26826000E+01 

ITERATION  NUMBER  15 

2  6  6  *  3 1 777684E+02  . 8  59  1 2800E  +  00 

366  .50569593E+02  . 1 55 731 OOE+O 1 
466  *610191 15E+02  .205B0090E+01 
566  .66207443E+02  .23523200E+01 
666  *677996 1 1E+02  . 24462330E+0 1 

ITERATION  NUMBER  16 

266  .32555933E+0?  . 77824900E+00 

3  6  6  *51 98  3078E+02  • 141 34850E  +  0 1 

4  6  6  .62890829E+02  .  1 87 1 7 1 40E  +  0 1 

5  6  6  .68349745E+02  . 2 1423020E+0 1 

666  .70028350E+02  . 22287390E+0 1 

ITERATION  NUMBER  17 

266  .33261224E+02  . 705291 OOE+OO 
366  •53266087E+02  • 12830090E+01 
4  6  6  *645925 16E+02  .  1 70 1 6870E+0 1 

566  .70299569E+02  . 1 9498240E+0 1 
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6  <1  6  .7P0575U5L+02  .20291550E+01 

I  If  RAT  ION  NUMBER  18 

2  6  6  • 3  39006  24  E  +  0  2  .63940000E+00 

366  .54430714E+02  . 1 16462 70E+0 1 

4  6  6  . 66 1 3  9 1 97E  +  02  .  1 54  668  1 OE  +  O  l  * 

566  .72073331E+02  . 1 77 37620E+0 1 

6  6  6  »  7390  39 1 2  E+02  . 1 8464070E  +  0 1 

ITERATION  NUMBER  19 

2  6  6  • 3448 0452 E+02  . 57982800E+00 

366  .55487917E+02  . 10572030E+0 1 
466  .67544674E+0?  . 14054770E+01 
566  #736B6290E+02  . 16 1 29590E+0 1 
666  .75583279E+02  . 16793670E+0 1 

ITERATION  NUMBER  20 

266  • 35006383E+02  . 52593 1 OOE+OO 

3  6  6  .5644  76  35  E+02  .95971800E+00 

466  .6882 1617E+02  . 12769430E+01 

5  6  6  .75  152556E+02  . 14662660E+0 1 

666  • 77 1 1 0 165E+02  . 1 5268860E+0 1 

ITERATION  NUMBER  21 

266  .35483510E+02  .47712700E+00 
366  .57318872E+02  . 8 7 1 23700E+00 
466  #6998 1606E+02  . 1 1599890E+01 
5  6  6  • 764851  35E+02  . 1 3325790E+01 

666  .78498025E+02  . 1 38 78600E+0 1 

ITERATION  NUMBER  22 

266  .35916428E+02  .43291800E+00 
366  #581 09802€+02  . 79093000E+00 

4  6  6  • 7 1035 240E+02  . 10536340E+01 

5  6  6  #7  769  5965  E+02  . 1 2 1 08300E+C 1 

6  6  6  .79759219E+02  . 126  1 1940E+01 

ITERATION  NUMBER  23 

2  6  6  •  36  309280E+02  . 39285200E+00 

3  6  6  #5882 7842E+02  . 7 1 8 04000E+00 
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4  6  6  • 71992  185E+02  . 95694500E+C0 

5  6  6  » 78 795996E+02  *  1 1000310E+01 
666  .80905099E+02  .  1 1458800E  +  0 1 

ITERATION  NUMBER  24 

266  • 36665808E+02  . 35652SOOE+OC 
366  • 594797 1 0E+02  .65 1 86B00E+00 
466  •  7286  1 2  39E+G2  . 86905400E+00 
566  . 79795230E+02  . 9992 3400E+00 
666  .81 946054E+02  . 10409550E+01 

ITERATION  NUMBER  25 

266  • 36989397E+02  . 32358900E+00 
3  6  6  .6007 1513E+02  .591803G0E+00 

466  .7*650439E+02  . 78920000E+00 
566  .80702814E+02  . 90758400E+00 
666  .8289 1585E+02  . 945 5 3 1 OOE+OO 

ITERATION  NUMBER  26 

266  .37283106E+02  . 293 70900E+00 
366  .60608790E+02  . 53727700E+00 
466  .74367081E+02  .71664200E+00 
5  6  6  .81 52  70B0E+02  • 82426600E  +  00 

666  .83750353E+02  . 858 76800E+00 

ITERATION  NUMBER  27 

2  6  6  . 375497 10E+02  . 26660400E+00 

366  .6 1096567E+02  .48777700E+00 
466  .7501 78 17E+02  .65073600E+00 
566  .82275630E+02  . 74855000E+00 
666  .84530266E+02  .77991300E+00 

ITERATION  NUMBER  28 

266  .37791716E+02  .24200600E+00 
366  .61539409E+02  . 44284200E+00 
466  • 75608685E+02  . 59086800E+00 
566  .82955381E+02  .67975100E+00 
666  .852385 12E+02  . 70824600E+00 

ITERATION  NUMBER  29 


141. 
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2  6  6  •  3  80 1 1 403  E  +  02  .  2  1968 700E  +  00 

3  6  6  .61941450E+02  .402041 OOE+OO 
466  .76145 185E+02  . 5 36 50000E+00 
566  .83572632E+02  .61725100E+00 
666  .8588 1652E+02  . 643 1 4000E+00 

ITERATION  NUMBER  30 

2  6  6  .382108  36E+02  . 1 9943300E  +  00 

366  .62305459E+02  . 36500900E+00 

4  6  6  . 76632  308E  +  02  .48712300E+00 

566  .841331 07E+02  . 56047500F+00 
666  . 86465647E+02  . 58399500E+00 

ITERATION  NUMBER  31 

2  6  6  . 38 39 1 882E+02  . 181 04600E  +  00 

3  6  6  .6263  7847E+02  . 33 1 38800E  +  00 
466  . 770745B9E+02  .44228100E+00 

5  6  6  .846420 14E+02  . 50890700E+00 

666  .86995924E+02  . 5 302 7700E+00 

ITERATION  NUMBER  32 

266  . 38556245E+02  . 16436300E+00 
3  6  6  .629387 02E+02  . 30085500E  +  00 

466  .77476 1 53E+02  .401 56400E+00 
566  .85104090E+02  .46207600E+00 
666  .87477400E+02  . 48 1 47600E+00 

ITERATION  NUMBER  33 

266  • 38705457E+02  . 14921200E+00 
366  .6321 1850E+02  . 273 14800E+00 
466  . 77840747E+02  . 36459400E+00 
566  .85523636E+02  .41954600E+00 
666  .8791 4566E+02  .437 16600E+00 

ITERATION  NUMBER  34 

266  • 38840922E+02  . 1 3546500E+00 
3  6  6  .634598 35E+02  . 24798500E  +  00 

466  .78 1 71 769E+02  . 331 02200E+00 
566  .85904559E+02  . 38092 300E+00 
666  .8831 1493E+02  • 39692700E+00 


14*. 


0A/Phy«/63-8 


ITERATION  NUMBER  35 

2  6  6  *3896 3905  E+02  .12298300^+00 
366  .63684976E+02  . 225 14 1 00E+00 
466  .78472308E+02  . 3005 3900E+00 
566  .86250416E+02  . 34585700E+00 
6  6  6  .8867 1880E+02  . 36038700E+00 

ITERATION  NUMBER  36 

2  6  6  • 3907 5556E+02  . 1 1 1651 OOE+OO 
366  .638893B2E+02  . 20440600E+00 
466  .78745170E+02  .27286200E+00 
566  .86564423E+02  . 31400700E+00 
6  6  6  .88999084E+02  . 327-20400E  +  00 

ITERATION  NUMBER  37 

266  • 39 1 76923Et02  . 101 36700E+00 
366  .640 74 960 E+02  . 1 8557800E+00 
466  • 78992909E+02  .24773900E+00 
566  .86849521E+02  .28509800E+0U 
666  .89296165E+02  .29708100E+00 

ITERATION  NUMBER  38 

266  .39268951E+02  .92028000E-01 
366  .6424 3445 E+02  . 16848500E+00 
466  • 792 1 7827E+02  , 22491800E+00 
566  .87108364E+02  . 25884300E+00 
666  .89565886E+02  .26972100E+00 

ITERATION  NUMBER  39 

26  6  • 39352499E+02  .83548000E-01 
366  . 643964 10E+02  • 1 5296500E+00 
466  • 79422030E+02  . 20420300E+00 
566  .87343369E+02  . 2 3500500E+00 
666  .8981 0770E+02  . 24488400E+00 

ITERATION  NUMBER  40 

266  .39428353E+02  . 75854000E-0 1 
366  .64535285E+02  • \ 3887500E+00 
466  .79607429E+02  . 18539900E+00 
5  6  6  .875567 30E+02  .21336100E+00 


141. 


GA/ffays/63-6 


6  6  6  .90033097E+02  .  222  327.00E  +  0G 

ITERATION  NUMBER  41 

266  .39497219E+02  • 68866000E-0 1  , 

3  6  6  .64661367E  +  02  .  12608200E  +  00 

466  .79775748E+02  . 16831900E+0G 
5  6  6  .87750442E  +  02  •  19371200E+00 

666  .90234955E+O2  . 201 85800E+00 

ITERATION  NUMBER  42 

266  .39559741E+02  .62522000E-01 
366  .64775836E+02  . 1 1446900E+00 

4  6  6  .79928569E+02  . 15282100E  +  00 

5  6  6  •87926315E  +  02  .  1 7587300E+00 

666  .90418219E+02  . 18326400E+00 

ITERATION  NUMBER  43 

266  • 396 16504E+02  . 56763000E-0 1 

3  6  6  .64879765E+02  .  10392900E  +  0G 

4  6  6  .80067315E+02  •  1  38 74600E+0U 

5  6  6  .88085990E+02  .  1 5967500E+00 

6  6  6  .90564606E+02  .  166 38900E  +  00 

ITERATION  NUMBER  44 

2  6  6  •39668038E+02  . 5 1 5 34000E-0 1 

366  .64974121E+02  . 943 56000E-0 1 
466  .B0193279E+02  . 12596400E+00 
566  .88230956E+02  . 14496600E+00 
666  .90735667E+02  . 15105900E+00 

ITERATION  NUMBER  45 

266  .39714828E+02  .46790000E-0 1 
366  •65059780E+02  ,856590006-01 

4  6  6  .80307644E+02  .  1 1436500E+00  i 

566  .88362575E+02  . 1 31 61900E+00 
6  6  6  .90872821E+02  .  1  37 15400E+00 

ITERATION  NUMBER  46 

266  • 39757305E+02  .42477000E-0 1 
366  .65137556E+02  .77776000E-01 


144, 


GA/Phya/63-8 


4 

5 
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2 

3 

4 

5 

6 


2 

3 

4 

5 

6 


2 

3 

4 

5 

6 


2 

3 

4 

5 

6 


2 

3 

4 

5 
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6  6  .8041 1472E+02 
6  6  .88482070E+02 
6  6  .90997336E+02 

ITERATION  NUMBER  47 

6  6  .39795870E+02 
6  6  .65208  166E+02 
6  6  .80505739E+02 
6  6  .885905  58E  +  02 
6  6  .911 1 0385E+02 

ITERATION  NUMBER  48 

6  6  • 39830885E+02 
6  6  •65272270’E+OP 
6  6  .80591 321 E+02 
6  6  .8868905 4E+02 
6  6  .91213024E+02 

ITERATION  NUMBER  49 

6  6  .39862671E+02 
6  6  .65330470E+02 
6  6  .80669023E+02 
6  6  .88778474E+02 
6  6  .91 306204E+02 

ITERATION  NUMBER  50 

6  6  .39891531E+02 
6  6  .653833 13E+02 
6  6  .80739565E+02 
6  6  .8885 9661 E+02 
6  6  • 9 1 390804E+02 

ITERATION  NUMBER  51 

6  6  .39917735E+02 
6  6  .6543 1285E+02 
6  6  .808036 12E+02 
6  6  .88933371E+02 
6  6  .91467613E+02 

ITERATION  NUMBER  52 


. 10382800E+00 
•  1 1949500E+00 
•12451 500E+00 


• 38565000E-0 1 
.7061 OOOOE-O 1 
• 94267000E-0 1 
. 10848800E+00 
•  1 1304900E+00 


•  350 1 5000E-01 
.64104000E-01 

•  85582  OOOE-0 1 
•98496000E-01 
• 10263900E+00 


•  31 786000E-0  1 
•58200000E-01 
•77702000E-01 
•89420000E-01 
.93180000E-01 


• 28860000E-0 1 
• 52843000E-0 1 
•  70542000E-0 1 
•81 187000E-01 
• 84600000E-0 1 


•26204000E-01 
•47972000E-01 
•64047000E-01 
•  7  37 10000E-0  1 
•76809000E-01 


I4C. 


<JA/Phy«/65-e 


2  6  6  .39941523E+02  . 2378 8000E-0 1 

366  .65474840E+02  .43555000E-01 
466  .80861762E+02  . 58 1 50000E-0 1 
566  • 89000295E+02  .66924000E-01 
666  .9153734SE+02  .69735000E-01 

ITERATION  NUMBER  53 

266  *39963 120E+02  . 2 1 5 97000E-C 1 
366  *655 14382E+02  . 39542000E-0 1 
466  •80914552E+02  . 52790000E-0 1 

5  6  6  *8906  1050E  +  02  . 60755000E-0 1 

666  .91600660E+02  .63312000E-01 

ITERATION  NUMBER  54 

2  6  6  *399827  30E+02  • 1 96 1 OOOOE-0 1 

366  .65550283E+02  . 35901000E-0 1 
466  .80962483E+02  . 4793 1 OOOE-O I 
566  *891 16210E+02  . 55 1 60000E-0 1 

6  6  6  *91658  140E  +  02  . 57480000E-0 1 

ITERATION  NUMBER  55 

2  6  6  .40000530E  +  02  .  1 78 COOOOE-0 1 

366  .65582878E+02  . 32595000E-0 1 
466  .81005999E+02  . 435 1 6009E-0 1 
566  .89166291E+02  .50081000^-01 
6  6  6  *917 1 0325E+02  • 52185000E-01 

ITERATION  NUMBER  56 

2  6  6  .40016691E+02  . 16161000E-01 

366  .65612470E+02  . 29592000E-0 1 

4  6  6  ■ 8  1 045 506E+02  . 39507000E-0 1 

5  6  6  *8921 1.7  58E+02  . 45467000E-0 1 
666  .91757701E+02  .47376000E-01 

ITERATION  NUMBER  57 

266  *4003 1 366E+02  . 14675000E-C1 
366  .65639336E+02  * 26866000E-0 1 
466  *8108 1 373E+02  . 35867000E-0 1 
566  .89Z53035E+02  .41277000E-0 1 
666  .91800718E+02  . 430 1 7000E-0 1 


1 44* 
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•  6 


2 

3 

4 

5 
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2 

3 

4 

5 

6 
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3 
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2 

3 

4 

5 

6 


3 

4 

5 


ITERATION  NUMBER  58 

6  6  .40044689E+02  . 1 332 3000E-0 1 
6  6  .65663728E+02  .24392000E-01 
6  6  *811 13936E+02  . 32563000E-0 1 
6  6  *892905 13E+02  . 37478000E-01 
6  6  •  9 1 8  39766E  +  02  . 39048000E-0 1 

ITERATION  NUMBER  59 

6  6  .40056783E+02  .  1 2094000E-0 1 
6  6  .65685872E+02  .22144000E-C1 
6  6  .81143500E+02  .29564000E-C1 
6  6  .89324538E+02  .34025000E-01 
6  6  *9187$226E+02  . 35460000E-0 1 

ITERATION  NIJMBFR  60 

6  6  .40067765E+02  . 10982000E-01 
6  6  «65705979E*02  .20107000E-01 
6  6  .81170341E+02  .26841000E-01 
6  6  .89355427E+02  . 30889000E-01 
6  6  .91907413E+02  .32187000E-01 

ITERATION  NUMBER  61 

6  6  .40077734E+02  . 99690000E-02 
6  6  .65724233E+02  . 18254000E-01 
6  6  .81194713E+02  . 243 72000E-0 1 
6  6  .89383472E+02  .28045000E-01 
6  6  .91936636E+02  .29223000E-01 

ITERATION  NUMBFR  62 

6  6  .400867B6E+02  . 90520000E-02 
6  6  .65740805E+02  *  165 72000E-01 
6  6  *812166 33E+02  .22120000E-01 
6  6  .89408934E+02  .25462000E-01 
6  6  *91 963 168E+02  .26532000E-01 

ITERATION  NUMBER  63 

6  6  .40095003E+02  .82170000E-02 
6  6  .65753850E+02  .15045000E-01 
6  6  *8l23692*lE+02  .20088000E-01 
6  6  .89432052E+02  .231 18000E-01 


1 47, 


GA/Phys/63-8 


666  .91987258E+02  . 24090000E-0 1 

ITERATION  NUMBER  64 

266  .40102464E+02  * 746 10000E-02 
366  .6576951QE+02  . 1 3660009E-0 1 

466  •  8  1 2  5  5  1 56E  +  02  . 18235000E-01  * 

566  .89453039E+02  . 209 8 7000E-0 1 
666  .92009131E+02  . 2 1 8 73000E-0 1 

ITERATION  NUMBER  65 

266  .40109237E+02  .67730000E-02 
3  6  6  *6578191  1E+02  . 124G10G0E-01 

466  *812717 17E+02  . 16561000E-01 
566  .89472095E+02  . 19056000E-01 
666  .92028985E+02  . 19854000E-01 

ITERATION  NUMBER  66 

2  6  6  *4011 53f  7E+02  . 6 1 5 00000E-02 

366  *65793 169E+02  . 1 1258000E-01 
466  .81286745E+02  . 1 5028000E-0 1 
566  .89489396E+02  . 1 7301 OOOE-0 1 
6  6  6  *9204 7015E+02  . 18030000E-01 

ITERATION  NUMBER  67 

266  .40120970E+02  . 558 30000E-02 
366  .65803392E+02  . 1O22300OE-01 
466  .81 300392E+02  . 1 364 7COOE-0 1 
566  *895051 00E+02  . 1 5704000E-0 1 
666  .92063380E+02  . 16365000E-01 

ITERATION  NUMBER  68 

266  .40126038E+02  . 50680000E-02 
366  .65812672E+02  . 92800000F-02 
466  .81312784E+02  . 12392000E-01 
566  *895 1 9360E+02  . 1 4260000E-0 1 
666  .92078239E+02  . 14859000E-01 

ITERATION  NUMBER  69  f 

266  .40130640E+02  .46020000E-02 
366  .65821096E+02  .84240000E-02 


>48 


GA/fchya/63-6 


4 

5 

6 


*  2 

3 

4 

•  5 

6 
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3 

4 

5 

6 


2 

3 

4 
3 
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2 

3 

4 

5 

6 


5 

6 


6  6  *81 324033E+02  .  1 1 249000E-C 1 
6  6  .89532308E+02  .  1 2948000E-0 1 
6  6  .92091727E+02  . 1 3488000E-0 1 

TERAT ION  NUMBER  70 

6  6  .40 1 3  48 18E+02  .41780000E-02 
8  6  .65828749E+02  . 76530000E-02 
6  6  *81 334246 E+02  . 1 0? 1 300PE-0 1 
6  6  .89544059E+02  . 1 1 75 1 OOOE-0 1 
6  6  *92 1 03976E+02  .  1 2249000E-0 1 

TERAT ION  NUMBER  71 

6  6  *401386J 1 E+02  . 37930000E-02 
6  6  .65835693E+02  . 69440000E-02 
6  6  *81 343517E+02  . 927 1 OOOOE-02 
6  6  .89554730E+02  . 10671000E-01 
6  6  .9211 5094E+02  . 1 1 1 18000F-01 

TERAT ION  NUMBER  72 

6  6  .40142055E+02  . 34440000E-02 
6  6  *6584 1997E+02  .63040000 E-02 
6  6  .81351936E+02  .84190000E-02 
6  6  .B9564420E+02  . 96900000E-02 
6  6  *921251 92E+02  . 1009800PE-0 1 

TERAT ION  NUMBER  73 

6  6  .40145183E+02  . 31280000E-02 
6  6  .65847723E+02  . 57260000E-02 
6  6  *81 359579E+02  . 76430000E-02 
6  6  *395732 16E+02  . 87960000E-02 
6  6  *92 134357E+02  .91650000E-02 

TERAT ION  NUMBER  74 

6  6  .40148021E+02  . 28380000E-02 
6  6  .65852923E+02  . 52000000E-02 
6  6  *81 3665 1 7E+02  . 69380000E-02 
6  6  .89581 199E+02  .79830000E  -02 
6  6  *92 142679E+02  .83220000E-02 

TERAT ION  NUMBER  75 


1 


0A/Phy»/63-8 


266  .40150599E+02  . 25780000E-02 
366  .65857640E+02  . 4 7 1 70000E-02 
466  .81372816E+02  . 62990000E-02 
566  .89588450E+02  . 725 1 0000E-02 
6  6  6  .92  150230E+02  . 755 1 0000E-02 

ITERATION  NUMBER  76  , 

266  .40152938E+02  .23390000E-02 

366  .65861925E+02  .42850000E-02  , 

466  .81 378536E+02  . 57200000E-02 
566  .89595029E+02  .6579GO00E-02 
666  .92157088E+02  .68580000E-02 

ITERATION  NUMBER  77 

266  .40155063E+02  .2 1250000E-02 
366  .658658 12E+02  . 3887COOOE-02 
466  .81383728E+02  . 5 192000PE-02 

5  6  6  .89601005E+02  . 59  760000F-02 

666  .92163316E+02  .62260000E-02 

ITERATION  NUMBER  78 

266  .40156991E+02  . 192B0P0CE-02 
366  .65869343E+02  . 353 1 OOOPE-02 
466  .81388442E+02  .47140nOCE-02 
566  .89606429E+02  . 5424C000E-C 2 
666  .92 1 68968E+02  . 56520000E-02 

ITERATION  NUMBER  79 

266  .40158743E+02  . 1 7520000E-02 
366  .65872551E+02  . 3208C000E-C2 
4  6  6  .813927 22E+02  .42800000E-02 

566  .8961135 9E+02  . 49290000E-02 

6  6  6  .92  174101E+02  . 5 1 3 30000E-02 

ITERATION  NUMBER  80 

266  .40160331E+02  . 1 5880000E-02  * 

366  .65875463E+02  .29120000E-02 
4  6  6  .813966  1 0E+02  . 38880000E-02 

566  .8961 5831E+02  .44730000E-02  * 

6  6  6  .92 178763E+02  .46620000E-02 


i  Co. 


t 


GAA  hya/63-8 


FINAL  RESULTS  OF  ITERATION 


node  temperature  residual 


2 

3 

,  4 

5 

6 

,  2 

3 

4 

5 

6 
2 

3 

4 

5 

6 
2 

3 

4 

5 

6 
2 

3 

4 

5 

6 
2 

3 

4 

5 

6 
2 

3 

4 

5 

6 
2 

*  3 

4 

5 

6 
2 
3 


2  2  • 1 1 3355  33E+02  *  1 8500000E-03 

2  2  • 16 193666E+02  .* 33900000E-03 

2  2  » 1 8492 290E+02  .45300000E-03 

2  2  .19555513E+02  . 52000000E-0 3 

2  2  • 19867755E+02  .54300000E-03 

3  2  •  1 6 1 9  3666  E  +  02  . 3 3900000E-0 3 

3  2  .24072778E+02  .61800000E-03 

3  2  *28065678E+02  .8280000PE-C3 

3  2  .29976634E+02  .95200000E-03 

3  2  .30545704E+02  . 99300000F-03 

4  2  •  1 649.2290E+02  .45300000E-03 

4  2  .28065678E+02  .82800000E-03 

4  2  • 33  129558E+02  . 1 10 70000E-02 

4  2  .35615056E+02  . 12730000E-02 

4  2  .36363747E+02  . 1 3260000E-02 

5  2  *195556  1 3E+02  . 52000000E-03 

5  2  .29976634E+02  . 95200000E-03 

5  2  *3561 5056E+02  .  1 2730000E-02 

5  2  .38423278E+02  . 1465C000E-02 

5  2  *3927513 1 E+02  . 1 5240000E-02 

6  2  • 19867755t+02  . 5430C000E-03 

6  2  .30545704E+02  .993C00G0E-03 

6  2  • 36363747E+02  . 1 3260000E-02 

6  2  *392751 31E+02  . 1 5240000E-02 

6  2  .40160331E+02  . 1 5880000E-02 

2  3  • 1 6 193666E+02  . 33900000E-03 

2  3  .24072778E+02  .61800000E-03 

2  3  .28065678E+02  .820OOOOOE-O3 

2  3  .29976634E+02  .95700000E-03 

2  3  .30545704E+02  .99300000E-03 

3  3  .24072779E+02  .61900000E-03 

3  3  •  37 124431E+02  • 1 1 340000E-02 

3  3  *44 129680E+02  . 1 5 1 50000E-02 

3  3  .4758071 1E+02  .  1  "»4  30000E-02 

3  3  .48620848E+02  . 181 70000E-02 

4  3  •  28065678E+02  .82800000E-03 

4  3  .44129680E+02  . 1 5 ] 50000E-02 

4  3  • 5 3076940E+02  . 2O24OOO0E-O2 

4  3  .57583003E+02  . 23290000E-02 

4  3  .58954961E+02  . 242 50000E-02 

5  3  .29976634E+02  . 95200000E-03 

5  3  *475807 1 1E+02  .17430000E-02 


f  SI. 

! 


GA/Phya/63-8 


4  5 

5  5 

6  5 

2  6 

3  6 

4  * 

5  6 

6  6 

2  2 

3  2 

4  2 

5  2 

6  2 

2  3 

3  3 

4  3 

5  3 

6  3 

2  4 

3  4 

4  4 

5  4 

6  4 

2  5 

3  5 

4  * 

5  5 

6  5 

2  6 

3  6 

4  6 

3  6 

6  6 

2  2 

3  2 

4  2 

5  2 

6  2 

2  3 

3  3 

4  3 

5  3 

6  3 

2  4 

3  4 

4  4 


3  .57583003E+02 
3  .62686205E+02 
3  .64249776E+02 
3  .30545704E+02 
3  *4862 0848E  +  02 
3  .58954961E+02 
3  • 6424  97  78E+02 

3  .6587 5463E+02 

4  « 1 8492289E+02 
4  .28065678E+02 
4  .33129558E+02 
4  .3561  505 6E+02 
4  .36  36  3747E  +  02 
4  .28Q65678E+02 
4  .44 12  9680E+02 
4  . 5  3076940E+02 
4  .57583003E+02 
4  . 5 8  95496 1 E+02 
4  .33  129558E+02 
4  . 53076940E+02 
4  .64567954E+02 
4  .70471552E+02 
4  .72285538E+02 
4  .3561 5056E+02 
4  .57583003E+02 
4  • 70471 552E+02 
4  .771719 12E+02 
4  .79242543E+02 
4  .36363747E+02 
4  .58954961E+02 
4  .72285538E+02 
4  .79242543E+02 

4  .81 3966 1 0E+02 

5  .1955551 3E  +  02 
5  .29976634E+02 
5  .  356 1 5056E+02 
5  .38423278E+02 
5  *392751 31E+02 
5  .29976634E+02 
5  .4758071 1E+02 
5  .57583003E+02 
5  .62686205E+02 
5  .64249778E+02 
5  .35615056E+02 
5  .57583003E+02 
5  .70471552E+02 


23290000E-02 
26760000E-02 
27950000E-02 
99300000E-03 
181 70000E-02 
24250000E-02 
27950000E-02 
29120000E-02 
45200000E-03 
8?800nOOE-03 
1 1070000E-02 
.  12730000E-02 
.  13260000E-02 
82800000E-03 
.151  50000E-02 
20240000E-02 
23290000E-02 

•  242  50000E-02 
.1  1070000E-02 

•  20240000E-02 

•  2  70  30000E-02 
•3111 OO0OE-O2 
.  32390000E-02 

•  1 2730000E-02 
.23290000E-02 
.311 10O0OE-02 
•358 1 0000E-O2 
•373 1 0000E-02 
.13260000E-02 
• 24250000E-02 
•32390000E-02 
•3731 0O00E-02 
. 38880000E-02 
.521 00000E-03 

•  95300000E-0  3 
. 12730000E-02 
. 14650000E-0? 
.15240000E-02 
.95300000E-03 
. 17430000E-02 
. 2  3290000E-02 
• 26780000E-02 
.27950000E-02 
. 12730000E-02 
•23290000E-02 
.311 10000E-02 


* 


r 


i 


ifaw 


OA/Phya/63-8 


DIFf-USIVlTY 
OELTA  TIME 


5  4  5 

6  4  5 

2  5  5 

3  5  5 

4  5  5 

5  5  5 

6  5  5 

2  6  5 

3  6  5 

4  6  5 

5  6  5 

6  6  5 

2  2  6 

3  2  6 

4  2  6 

5  2-6 

6  2  6 

2  3  6 

3  3  6 

4  3  6 

5  3  6 

6  3  6 

2  4  6 

3  4  6 

4  4  6 

5  4  6 

6  4  6 

2  5  6 

3  5  6 

4  5  6 

5  5  6 

6  5  6 

2  6  6 

3  6  6 

4  6  6 

5  6  6 

6  6  6 

1*0000 

•  2000 


771719 12E+02 
79242543E+02 
38423278E+02 
62686205E+02 
77171912E+02 
84787288E+02 
87153425E+02 
392751 31 E+02 
64249778E+02 
79242543E+02 
8  715  3425E+02 
896  15831  E  +  02 
19867755E+02 
30545704E+02 
36363747E+02 
392751 31E+02 
40160331E+02 
30545704E+02 
4B620848E+02 
5895496 1 E+02 
64249778E+02 
65875463E+02 
i  36363747E+02 
,5895496 1 E+02 
,  72285538E+02 

>  79242543E+02 
,  8 1 3966  10E+02 

>  392751 3 1 E+02 
.64249778E+02 
« 79242543E+02 
.87153425E+02 
.89615831E+02 
.40160331E+02 
.65875463E+02 
1 8 1 3966 10E+02 
>8961 5831E+02 
>92 1 78763E+02 

CONDUCT  IV 
SOURCE  STRl 


.35810000F-02 
.3731 J000E-C2 
.  14650000E-0? 
•26780000E-02 
•358 1 OOOOE-02 
•41220000E-0? 
•42930000E-02 
• 15240000E-02 

•  27950000E-02 
•37310000E-02 
.42910000E-02 
•44730000E-02 
•54300000E-03 
•9930C000E-03 
p 1 326C000E-02 
• 15240000E-02 
•15880000E-02 
•9930CO00E-03 
• 18170000E-02 
• 24250000E-02 
•27950000E-02 
•29120000E-02 
•13260000E-02 
•24250000E-02 
•32390000F-02 
•3731 0000E-02 
•38880000E-02 

•  1 5240000E-02 
•27950000E-02 
•3731 0000E-02 
•42930000E-02 
•44730000F-02 
•15880000E-02 
•29120000E-02 
•38880000E-02 
.44730000E-02 
•46620000E-02 
TY  1.0000 
NOTH  2000.0000 


XGR 1 0*  .1000  YGR I D«  .1000  ZGRID*  .1000 

MAX  TIME  .2000STOP  WHEN  ABSDIFF-  .50000000E-02  IS  ZERO  OR  LESS 


OUTPUT  CONTROLS  ARE  IRITE 
KRITE  1  K  ITER  6  JITER 


1  JRITE  1 
6  CHOSE  0.000 


IS*. 


GAAhys/6?5-e 

DATA  R  IJ»4  REGULAR  WEGSTEIN  —  T  -  0.2  SEC,  GRIDS  EQUAL  AT  0.1  CM 


NEW  T I  ME  INCREMENT  .200 

ITERATION  NUMBER  1 

NODE  TEMPERATURE  RESIDUAL  SL0PF. 

2  6  6  .49352900E+01  .4935290OE+O1 

3  6  6  .615094/1E+01  .61509471E+01 

4  6  6  .64493826E+01  .b4493826E+01 

5  b  6  .6522281 8E+01  .65222818E+01 

666  .65399740E+01  .65399740E+01 

ITERATION  NUMBER  2 

266  .88801 862E+01  .3944896 2 E+01 

3  6  6  . 1 1733575E+02  .55826280E+01 

466  . 1 2595459E+02  .61460770E+01 
566  . 1 2844798E+02  .632251 70E+01 
666  .12909294E+02  .63693200E+01 

ITERATION  NUMBER  3 

266  . 1 2185228E+02  .33050420E+01 
366  . 16783346E+02  .50497710E+01 
46b  . 1 8385662E+02  .57902030E+01 

5  6  6  . 18910285E+02  .60654870E+01 

666  .1 90521 05E+02  .614281 10E+01 

ITERATION  NUMBER  4 

266  . 1 5032849E+02  .28476210E+01 
366  .21351303E+02  .45679570E+01 
466  .23794226E+02  .54O85640E+O1 
566  . 24671 974E+02  .57616890E+01 
666  .24916819E+02  . 58647 140E+01 

ITERATION  NUMBER  5 

266  .17527906E+02  .24950570E+01 
366  .25486291 E+02  .41 349880E+01 

4  6  6  ,288l4l65E*f02  .50199390E+01 

566  .30098025E+02  . 542605 10E+01 
666  . 30464459E+02  .55476400E+01 

ITERATION  NUMBER  6 

2  6  6  .1 97368 19E+02  . 220891 30E+01 

366  . 29231 800E+02  .37455090E+01 
*♦  6  6  .33450971E+02  .  46368060E+01 

9  6  6  .35171 236E+02  .507321  OE-MH 

6  6  6  .35670706E+02  .52062470E+01 

104. 


GAAhys/65-8 


ITERATION  NUMBER  7 

266  .21705484E+02  . 19686650E+01 
366  .32626222E+02  . 33944220E+01 
466  .37718101E+02  .42671300E+01 
566  .39886695E+02  .47154590E+01 
666  .40524592E+02  .48538860E+01 

ITERATION  NUMBER  8 

266  .23467933E+02  . 17624490E+O1 
366  .35703599E+02  .3O773770E+O1 
466  .41633842E+02  .391 5741 OE+Ol 
566  .44248880E+02  .43621 850E+01 
6  6  6  .45025903E+02  .450131  OE+Ol 

ITERATION  NUMBER  9 

266  .25050775E+02  . 1 5828420E+01 
366  .384943 2 3E+02  . 27907 240E+01 
466  .452191 52E+02  .35853100E+01 
566  .48268970E+02  .40200900E+01 
666  .49182403E+02  .41565000E+01 

ITERATION  NUMBER  10 

266  .26475626E+02  .1424851 OE+Ol 
366  .41 025641 E+02  .25313180E+01 
466  .48496204E+02  .32770520E+01 
566  .51962594E+02  .36936240E+01 
666  .53007382E+02  .38249790E+01 

ITERATION  NUMBER  11 

266  . 277605 16E+02  . 1 2848900E+01 
366  .43322061 E+02  .22964200E+01 
466  .51487394E+02  .2991  900E+01 
566  .55348087E+02  .33854930E+01 
666  .56517688E+02  .35103060E+01 

ITERATION  NUMBER  12 


ITERATION  NUMBER  13 

266  .44852152E+02  .101  0451E+02 
366  .86393980E+02  .2341 7471 E+02 
466  .78428443E+02  .15310434E+02 
966  .73398767E+02  .1 054081 7E+02 
666  .70871016E+02  .89507020E+01 


105. 


GAAtys/63' 


-e 


ITERATION  NUMBER  14 

266  .51 36857 1E+02  .65164190E+01 
366  .85345 189E+02-.1 048791 0E+01 
466  .83703602E+02  .52751 590E+01 
566  . 80649055 E+02  .72502880E+01 
666  .79693043E+02  .88220270E+01 

ITERATION  NUMBER  15 

266  .51339767E+02-.28804000E-01 
3  6  6  .77532625E+02-.781 25640E+01 

466  .84251 769E+02  .54816700E+00 
5  6  6  .85209332E+02  . 456027 70E+01 

666  . 855495 36E+02  .58564930E+01 

ITERATION  NUMBER  16 

266  .47631 880E+02-.37078870E+01 
3  6  6  .74328679E+02-. 32039460E+01 

466  . 847505 29E+02  .49876000E+00 
566  .883091 44E+02  .309981 20E+01 
666  .89347939E+02  .37984030E+01 

ITERATION  NUMBER  17 

266  .45810482E+02-.18213980E+01 
366  .72268586E+02-.20600930E+01 
466  .84915766E+02  .16523700E+00 
566  .90204296E+02  . 18951520E+01 
666  .91824796E+02  .24768570E+01 

ITERATION  NUMBER  18 

3  6  6-. 1341  200E+03-. 20638058E+03  .99188337E+00 

ITERATION  NUMBER  19 

266  .85797480E+01-.35130178E+02 
366  . 463691 07E+02  .180481  0E+03 

4  6  6  .75833670E+02  .2504171 3E+02 

566  .88148656E+02  .24035820E+01 
666  .89534468E+02-. 10191980E+01 

ITERATION  NUMBER  20 

266  . 34557 193E+02  .25977445E+02 
3  6  6  .51 224439E-KJ2  .48553320E+01 

466  . 74969 384E+02-.86428600E+00 
566  .86551 262E+02-.15973940E+01 
666  .89934059E+02  . 399591 OOE+OO 

156. 


C A/lhy s/ 6?-8 


ITERATION  NUMBER  2 1 

2  6  6  .336741 01 E+02-.88309200E+00 

3  6  6  .59306274E+02  .80818350E+01 

466  .776460/7E+02  .26766930E+01 
566  .87125141 E+02  .57387900E+00 
666  .89426491 E+02-.50756800E+00 

ITERATION  NUMBER  22 

266  .371 24091 E+02  .34499900E+01 
366  .61684080E+02  .23778060E+01 
466  .78234057E+02  .58798000E+00 

5  6  6  .86822557E+02-.30258400E+00 

6  6  6  .89560956E+02  . 1 3446500E+00 

ITERATION  NUMBER  23 

266  . 378501 49E+02  .72605800E+00 
366  .62959550E+02  .12754700E+01 
466  .78700857E+02  . 46680000E+00 
566  .87051002E+02  .22844500E+00 
666  .895241 59E+02-.36797000E-01 

ITERATION  NUMBER  24 

366  .67402930E+02  .44433800E+01  .77697O84E+O0 

ITERATION  NUMBER  25 

266  .3971 2428E+02  .10285740E+01 
366  .64827 1 60E+02-. 25757700E+01 

4  6  6  ,79679853E-^02-.87063000E-01 

5  6  6  .87572706E+02  . 24242500E+00 

666  .90034450E+02  . 25863 100E+00 

ITERATION  NUMBER  26 

266  .39474969E+02-.23745900E+00 
366  .64796037E+02-.31  23000E-01 
466  .79854291 E+02  .17443800E+00 
5  6  6  .87786705E402  .21399900E+00 

666  .90260794E+02  .22634400E+00 

ITERATION  NUMBER  27 

266  . 39567 337E+02  .92368000E-01 
366  .64803561 E+02  .75240000E-02 
466  .79942841 E+02  .88550000E-01 
566  .87940229E+02  .15352400E+00 


157 


Cil/thj-e/6  3-C 

66b  .90443049 E+02  .18225S00E+00 

ITERATION  NUMBER  28 

266  .39596002E+02  .28665000E-01 
366  .6486591 4E+02  .62353000E-01 

4  6  6  .80058366E+02  .1  552500E+00 

966  .88090688E+02  . 1 5045900E+00 
666  .90593803E+02  . 1 5075400E+00 

ITERATION  NUMBER  29 

266  .39644080E+02  .48078000E-01 
366  .64943841 E+02  .77927000E-01 
466  .801 70791 E+02  .1  242500E+00 
566  .88221 41 7E+02  . 1 3072900E+00 
666  .90733244E+02  .1 39441 OOE+OO 

ITERATION  NUMBER  30 

6  6  6  .10174910E+03  .1 '015860E+02  .98750008E+00 

ITERATION  NUMBER  31 

266  .40534931E+02  .49760100E+00 
366  .65530987E+02  .33360400E+00 
466  . 80543 187E+02  . 20694300E+00 

5  6  6  .90331534E+02  J9631770E+01 

666  .91917293E+02-.98318100E+01 

ITERATION  NUMBER  32 

266  . 4044908 3E+02-.85848000E-01 
366  .65608256E+02  .77269000E-01 
466  .80981496E+02  .43830900E+00 
566  .89076081E+02-.12554530E+01 
666  .92268235E+02  .35094200E+00 

ITERATION  NUMBER  33 

266  .40187156E+02-.26192700E+00 
3  6  6  .65619279E+02  .P023000E-01 

466  i 80895 394E+02-. 86 102000E-01 
966  .89281 976E+02  . 20589500E+00 
666  .91747795E+02-.52044000E+00 

ITERATION  NUMBER  34 

266  .401 1 3039E+02-.741 17000E-01 
366  .6S607S79E+02-.1  700000E-01 
466  ,8101 23H7E+02  .1  699300E+00 
566  .89179240E+02-.10273600E+00 
666  .91798949E+02  .51  54000E-01 

_18B. 


GA/Fhys/63-8 


ITERATION  NUMBER  35 

2  6  6  . 400692 29E+02-. 438 10000E-01 

3  6  6  .65632031 E+02  .24452000E-01 

4  6  6  .81031  34E+02  . 18747000E-01 

566  .89235540E+02  .56300000E-01 
6  6  6  .91 775925E+02-, 23024000E-01 

ITERATION  NUMBER  36 

266  .40047945E+02-.21284000E-01  .32695731E+00 
366  .65659270E+02  .27239000E-01  .52699165E+00 

5  6  6  .892921 20E+02  .56580000E-01  .501261 10E+00 

ITERATION  NUMBER  .  37 

266  . 40055 292E+02  .73470000E-02 
366  .65676297E+02  . 17027000E-01 
466  .811 34991 E+02  ,4l078000E-0i 
566  .8931 441 8E+02  .22298000E-01 
666  .91870999E+02  .39107000E-01 

ITERATION  NUMBER  38 

266  .40065759E+02  . 10467000E-01 
366  .657001 44E+02  .23847000E-01 
466  .81 160295E+02  .25304000E-01 
566  .89347487E+02  .33069000E-01 
666  .91898458E+02  .27459000E-01 

ITERATION  NUMBER  39 

266  .40075656E+02  .98970000E-02 
366  .65718252E+02  . 18108000E-01 
466  .811 8691 9E+02  .26624000E-01 
566  .89374843E+02  .27356000E-01 
666  .91928420E+02  .29962000E-01 

ITERATION  NUMBER  40 

266  .40084553E+02  .88970000E-02 
366  .65735925E+02  .1767300OE-OI 
466  .81209764E+02  .22845000E-01 
566  .89401 222E+02  .26379000E-01 
666  .91955065E402  .26645000E-01 

ITERATION  NUMBER  41 

2  6  6  ,40092966E-K)2  .84130000E-02 

366  .65751456E+02  .15531000E-01 
466  .81230797E+02  .21033000E-01 
5  6  6  .8942495 1E-MJ2  .23729000E-01 

666  .91 97991 4E+02  . 24849000E-01 

159. 


GA/Fhys/63-8 


ITERATION  NUMBER  42 

266  .40664700E+02  .571 73400E+00  .98549863E+00 
ITERATION  NUMBER  43 
266  .40155975E+02-.50872500E+00 

366  .65941 433E+02  .81778000E-01  i 

466  .81918920E+02  . 23045 200E+00 
5  6  6  .90107529E+02  .32909000E+00 

666  .92504389E+02  .33788200E+00  < 

ITERATI0N  NUMBER  44 

266  .40192710E+02  .36735000E-01 
366  .660481 72E+02  . 10673900E+00 
466  ,81 743366E+02-. 17555400E+00 

5  6  j  .900857/6E+02-.21753000E-01 

6  6  6  .92657030E+02  . 15264100E+00 

ITERATION  NUMBER  45 

266  . 402027 26E+0 2  . 10016000E-01 
366  .6601 8333E+02-. 29839000E-01 
466  ,81 726789E+02-. 16577000E-01 
566  .90038601 E+02-. 47 175000E-01 
666  .92623301 E+02-.33729000E-01 

ITERATION  NUMBER  46 

266  .40206900E+02  .41 740000E-02 
366  .66031618E+02  . 1 3285000E-01 
4  6  6  .81 677748E+02-, 49041 OOOE-OI 

566  .89980809E+02-.57792000E-01 
666  .92571394E+02-.51907000E-01 

ITERATION  NUMBER  47 

266  .40216700E+02  .98000000E-02 
366  .660201 83E+02-. 1 ,435000E-0t 
466  .81643890E+02-.33858000E-01 

566  .89934384E+02-.46425000E-01  1 

666  .92519805E-K12-.51589000E-01 

ITERATION  NUMBER  48  ' 

2  6  6  .402187 1 4E+02  .20140000E-02  . 1 7061 224E+00 

366  .65997780E+02-. 22403000E-01  .66200262E+00 


160. 


I 


» 


t 


ITERATION  NUMBER  49 

2  6  6  ,402171  37E+02-.  1 5770000E-02 

3  6  6  .65992337F.+02-.54430000E-02 

4  6b  .81535184E+02-.58734000E-01 

5  6  6  .89807301 E+02-. 6691 0000E-01 

<>6  6  .924041 93E+02-.62883000E-01 

ITERATION  NUMBER  50 

2  6  6  .40214339E+02-.27980000E-02 

3  6  6  .65984920E+02-.2741 7000E-01 

488  .81505052E+02-.30132000E-01 
5  6  8  .897690/2E+02-.38229000E-01 

8  8  8  .92352080E+02-.52073000E-01 

ITERATION  NUMBER  51 

2  .8  6  .402053/8E+02-.8961 0000E-02 

368  .65949978E+02-.14942000E-01 
488  .81500347E+02-.47050000E-02 
588  .89751201E+02-.17871000E-01 
8  6  8  .92324220E+02-.27860000E-01 

ITERATION  NUMBER  52 

268  .40199481 E+02-.591 70000E-02 
388  .65945630E+02-.43480000E-02 
468  .81494320E+02-.80270000E-02 

5  6  6  .89736533E+02-. 1 4668000E-01 

6  6  6  .92308033E+02-. 1 6l 87000E-01 

ITERATION  NUMBER  53 

268  .401 96567E+02-.28940000E-02 
366  .65941 553E+02-.40770000E-02 

4  6  6  .81488275F.+02-.60450000E-02 

566  .8972671 2E+02-.98210000E-02 
666  .92295861 E+02-.1 21 72000E-01 

ITERATION  NUMBER  54 

466  .81432630E+02-.55645000E-01  .90190239E+00 

ITERATION  NUMBER  55 

266  .401899 l 4E+02-.23360000E-02 
366  .65925560E+02-, 1  841000E-01 
466  .81468953E+02  .36323000E-01 
566  ,89694835E-K)2-.56250060E-02 
666  .922621 28E+02-. 1 5831000E-01 


161. 


GA/Fhys/63-e 


ITERATION  NUMRER  56 


2  6  6  .40188478E+02-.14360000E-02 

3  6  6  .6592H388E+02  .28280000E-02 

46b  .81463657E+02-.52960000E-02 
966  .89693365E+02-.14700000E-02 
6  6  6  .92258601 E+02-.35270000E-02 

ITERATION  NUMBER  57 

266  .40187951 E+02-.52700000E-03 
366  .65925081 E+02-.33070000E-02 
466  .81462624E+02-. 10330000E-02 

5  6  6  .89690499E+02-.28660000E-02 

6  6  6  .92256453E+02-. 21 480000E-02 

ITERATION  NUMBER  58 

2  6  6  .40186625E+02-. 1 3260000E-02 

366  .65923581 E+02-.15000000E-02 
466  .81459925E+02-.26990000E-02 
566  .896881 34E +02-. 23650000E-02 
666  .92253860E+02-.25930000E-02 


FINAL  RESULTS  OF  ITERATION 


2  2  2  .1  338654E+02 

422  .18499369E+02 
622  .19876607E+02 
262  .1 987665 2E+02 
4  6  2  .36385764F.+02 

6  6  2  .401 8662 2F.+02 

226  .19876605E+02 
426  .36385683E+02 
6  2  6  .40186572F.+02 

266  .40186625E+02 
466  .81459925E+02 
6  66  .92253860F+02 

OIFFUS  IVITY  1.0000  CONDUCTIVITY  1.0000 

OELTA  TIME  .2000  SOURCE  STRENGTH  2000.0000 

XGR 1 0  . 1 0IK)  YGRI D-  .1000  ZGRID-  .1000 

MAX  Tl  ME  .2000ST0P  WHEN  ABSDIFF-  .005001 S  ZERO  OR  LESS 

OUTPUT  CONTROLS  ARE  (RITE  2  JR  I TE  4 
KRITE  4  KITER  6  JITER  6 
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Vita 

Hobert  Theodore  roppe  was  born  on  1?  February  1950  in 
Covington,  Kentucky,  the  son  of  Huth  A,  xoppe  ana  the  late 
Raymond  H.  roppe.  After  completing  his  work  in  1946  at 
Newport  nigh  school,  Newport,  Kentucky,  he  did  undergraduate 
work  at  the  University  of  Cincinnati,  Cincinnati,  Ohio  and 
Morehead  state  College,  Morehead,  Kentucky  until  January 
1901  when  he  enlisted  in  the  U.S.  Air  Force.  He  was  dis¬ 
charged  in  June  of  tfct  same  year  to  enter  the  U.  a.  Naval 
Academy  at  Annapolis,  Maryland.  Capt.  Poppe  graduated  in 
June  1955  with  a  degree  of  Bachelor  of  Science  and  was  com¬ 
missioned  a  second  Lieutenant  in  the  U.  6.  Air  Force.  After 
graduation,  he  attended  the  Ground  Electronics  school  at 
Kessler  Air  Foroe  Base,  Mississippi.  During  the  period  from 
1956  till  1961,  he  was  assigned  to  maintenance  activities 
in  ground  radar  and  classified  work  in  the  field  of  elec¬ 
tronics.  In  August  of  1961  he  entered  the  two  year  graduate 
program  in  Astronautics  at  the  Air  Foroe  Institute  of  Tech¬ 
nology  «e si dent  school  at  wright-Patterson  A.F.6. ,  Ohio. 


